proj geo
|
Contents: Description, Arguments, Usage, Examples, Sub-Functions, Related Functions, Source Supported Under Version: > 1.02
Description Project/re-project a Geo-Struct or Geo-Aware file. This function relies on GDAL > 1.8. Note: Currently there is an issue with re-projecting ISIS data that are in a SimpleCylindrical Projection.
Arguments and Return Values Arguments: Geo-Struct or filename Return Value: A reprojected davinci Geo-Struct
Usage Syntax: proj_geo(Geo-Struct, t_srs [,match=GeoStruct][,te=VAL][,resolution=FLOAT][,resampling=STRING][,ignore=FLOAT][,opts=STRING][,a_srs = STRING][,subdata = INT][,raw = BOOL]) or Syntax: proj_geo(filename, t_srs [,match=GeoStruct][,te=VAL][,resolution=FLOAT][,resampling=STRING][,ignore=FLOAT][,opts=STRING][,a_srs = STRING][,subdata = INT][,raw = BOOL]) Geo-Struct - davinci Geo-Struct filename - input full path to data t_srs - target projection argument (either a proj string or wkt) resolution - resolution of output in georeferenced units (Default=native) resampling - resampling type (near,bilinear,cubic,cubicspline,lanczos) ignore - ignore value for the data (Default="") opts - additional gdalwarp options (Default="") a_srs - source projection argument (either a proj string or wkt), assigned by default values in Geo-Struct te - extent xmin//ymin//xmax//ymax, assigned by default values in Geo-Struct match - match the projection/extent/resolution of the supplied match Geo-Struct, with overrides from other parameters raw - return the raw gdalinfo output in the structure subdata - write out only the subdataset specified. Currently we only support the writing of 1 subdataset.
Examples dv> a=read_geo($1)
struct, 14 elements
anc: struct, 11 elements
wkt: "PROJCS["SinusoidalMars",GEOGCS["GCS_Mars",DATUM["D_Mars",\
SPHEROID["Mars",3396190,0]],PRIMEM["Reference_Meridian",0],\
UNIT["degree",0.0174532925199433]],PROJECTION["Sinusoidal"],\
PARAMETER["longitude_of_center",315],\
PARAMETER["false_easting",0],\
PARAMETER["false_northing",0]]"
proj4: "'+proj=sinu +lon_0=315 +x_0=0 +y_0=0 +a=3396190 +b=3396190 +units=m +no_defs '"
origin: 1x2x1 array of double, bsq format [16 bytes]
pixel_size: 1x2x1 array of double, bsq format [16 bytes]
x: 5
y: -2983
lat: 1x5x1 array of float, bsq format [20 bytes]
lon: 1x5x1 array of float, bsq format [20 bytes]
lpo: 1x5x1 array of float, bsq format [20 bytes]
spo: 1x5x1 array of float, bsq format [20 bytes]
center: 2x1x1 array of float, bsq format [8 bytes]
ignore: -32768
projection: "Sinusoidal"
spheroid: struct, 3 elements
name: "Mars"
radius: "3396190"
iflattening: "0"
resolution: 100.0000000
x: 5
y: -2983
ul: 2x1x1 array of float, bsq format [8 bytes]
ur: 2x1x1 array of float, bsq format [8 bytes]
ll: 2x1x1 array of float, bsq format [8 bytes]
lr: 2x1x1 array of float, bsq format [8 bytes]
w: 736
h: 2966
dv> out=proj_geo(a,"+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=3396190 +b=3396190 +units=m +no_defs",v=1,ignore=-32768)
struct, 14 elements
anc: struct, 11 elements
wkt: "PROJCS["unnamed",GEOGCS["unnamedellipse",DATUM["unknown",\
SPHEROID["unnamed",3396190,0]],PRIMEM["Greenwich",0],\
UNIT["degree",0.0174532925199433]],PROJECTION["Equirectangular"],\
PARAMETER["latitude_of_origin",0],\
PARAMETER["central_meridian",0],\
PARAMETER["standard_parallel_1",0],\
PARAMETER["false_easting",0],\
PARAMETER["false_northing",0],\
UNIT["metre",1,AUTHORITY["EPSG","9001"]]]"
proj4: "'+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=3396190 +b=3396190 +units=m +no_defs '"
origin: 1x2x1 array of double, bsq format [16 bytes]
pixel_size: 1x2x1 array of double, bsq format [16 bytes]
x: -26644
y: -2980
lat: 1x5x1 array of float, bsq format [20 bytes]
lon: 1x5x1 array of float, bsq format [20 bytes]
lpo: 1x5x1 array of float, bsq format [20 bytes]
spo: 1x5x1 array of float, bsq format [20 bytes]
center: 2x1x1 array of float, bsq format [8 bytes]
ignore: -32768.00000
projection: "Equirectangular"
spheroid: struct, 3 elements
name: "unnamed"
radius: "3396190"
iflattening: "0"
resolution: 100.0912704
x: -26644
y: -2980
ul: 2x1x1 array of float, bsq format [8 bytes]
ur: 2x1x1 array of float, bsq format [8 bytes]
ll: 2x1x1 array of float, bsq format [8 bytes]
lr: 2x1x1 array of float, bsq format [8 bytes]
w: 747
h: 2963
data: 747x2963x10 array of float, bip format [88,534,440 bytes]
|
DavinciWiki Mini-Nav Bar Contents
Contact Developers
All other topics
Related Functions
Recent Library Changes Created On: 04-11-2013 |
