Hi all, I'm working on a batch script
to produce warped sst images of a specific region from avhrr files, the
calibration part works fine but I'm having problems with the warping
part, I need to convert to a mercator projection but am unsure of the
parameters that need to be entered (leading to my script failing) this
is the part that is failing, where am I going wrong?
Secondly is it possible to warp from input geometry, would that be a
better way of going about the final warping and georeferancing, if it is
how should I approach it.
Thanks
Neil Ackermann
;Warp the calibrated data to Sample Mercator
;
; Set the keywords
; Query file to get the size
envi_file_query, fid, ns=ns, nl=nl, nb=nb
;set the dimensions to be warped
dims = [-1l, 0, ns-1, 0, nl-1]
;set which band to be warped
pos = [0]
out_name = tmp3_name
units = envi_proj_units_translate('Meters')
; proj = envi_proj_create(/geographic, unit=units)
; Parameters are a, b, lat0, lon0, x0, y0
Params = [6378160.0, 6356774.7, 0.000000, 18.000000, 500000., 10000000.]
datum = 'cape'
name = 'West Coast Mercator'
proj = envi_proj_create(type=20, name=name, datum=datum, params=params)
;pixel size of 1100m by 1100m
ps = [1100,1100]
;
; Call the doit
;
envi_doit, 'envi_avhrr_warp_doit', fid=cl_fid, $
avhrr_fid=fid, dims=dims, pos=pos, $
out_name=out_name, method=8, zero_edge, $
background=0, grid=[51,500], proj=proj, $
pixel_size=ps
; ,r_fid=r_fid
|