begin

? ? ;y = cbinread("./N22E113.hgt",-1,"short")

? ? ;printVarSummary(y)

setfileoption ("bin", "ReadByteOrder", "BigEndian")

demdata = cbinread("./N22E113.hgt",(/3601,3601/),"short")

;print(demdata(:,0))

;print(max(demdata))

;print(min(demdata))

lat_s = 22.00

lat_n = 23.00

lon_w = 113.0

lon_e = 114.0

;nlat = 3601

;nlon = 3601?

;lat1d = fspan(lat_s,lat_n,nlat) ;---結(jié)果出現(xiàn)坐標(biāo)間距可變bug

;lon1d = fspan(lon_w,lon_e,nlon)

dxy = 1d0/3600? ;---雙精度---

fxy = fspan(0, 3600*dxy, 3601)?

lon1d? = fxy + lon_w

lat1d? = lat_n - fxy ;---y reverse

demdata!0 = "lat"

demdata!1 = "lon"

demdata&lat = lat1d

demdata&lon = lon1d

demdata&lat@units = "degrees_north"

demdata&lon@units = "degrees_east"

;---write---

system("/bin/rm -f PRDdem_o.nc")? ; remove any pre-existing file

ncdf = addfile("PRDdem_o.nc" ,"c")? ; open output netCDF file? ?

ncdf->z = demdata

end