I am very much so a beginner.. but here goes. I have two maps that should line up with each other, but their projections are different. I've tried to use the gdalwarp utility to match them up, but it doesn't seem to be working. Here's some info
Sus_10m.tif Driver: GTiff/GeoTIFF Files: ../dems/Sus_10m.tif Size is 2252, 2115 Coordinate System is: PROJCS["REGVEN / UTM zone 18N", GEOGCS["REGVEN", DATUM["Red_Geodesica_Venezolana", SPHEROID["GRS 1980",6378137,298.2572221010002, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6189"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433], AUTHORITY["EPSG","4189"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-75], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AUTHORITY["EPSG","2201"]] Origin = (238238.500000000000000,4514375.500000000000000) Pixel Size = (10.000444049733570,-10.001891252955083) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates:
Upper Left ( 238238.500, 4514375.500) ( 78d 6' 0.03"W, 40d44'19.34"N)
Lower Left ( 238238.500, 4493221.500) ( 78d 5'28.35"W, 40d32'54.29"N)
Upper Right ( 260759.500, 4514375.500) ( 77d50' 1.03"W, 40d44'44.02"N)
Lower Right ( 260759.500, 4493221.500) ( 77d49'32.07"W, 40d33'18.80"N)
Center ( 249499.000, 4503798.500) ( 77d57'45.37"W, 40d38'49.39"N)
Band 1 Block=2252x1 Type=Float64, ColorInterp=Gray NoData Value=nan
GRASS 6.4.4 (demolocation)> gdalinfo tmax_1980_11752.nc | head -n 50 Driver: netCDF/Network Common Data Format Files: tmax_1980_11752.nc Size is 212, 248 Coordinate System is: PROJCS["unnamed", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], TOWGS84[0,0,0,0,0,0,0], AUTHORITY["EPSG","6326"]], PRIMEM["Greenwich",0, AUTHORITY["EPSG","8901"]], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9108"]], AUTHORITY["EPSG","4326"]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["standard_parallel_1",25], PARAMETER["standard_parallel_2",60], PARAMETER["latitude_of_origin",42.5], PARAMETER["central_meridian",-100], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (1718000.000000000000000,218000.000000000000000) Pixel Size = (1000.000000000000000,-1000.000000000000000) Metadata: lambert_conformal_conic#false_easting=0 lambert_conformal_conic#false_northing=0 lambert_conformal_conic#grid_mapping_name=lambert_conformal_conic lambert_conformal_conic#latitude_of_projection_origin=42.5 lambert_conformal_conic#longitude_of_central_meridian=-100 lambert_conformal_conic#standard_parallel={25,60} NC_GLOBAL#citation=Please see http://daymet.ornl.gov/ for current Daymet data citation information NC_GLOBAL#Conventions=CF-1.4 NC_GLOBAL#references=Please see http://daymet.ornl.gov/ for current informatio n on Daymet references NC_GLOBAL#source=Daymet Software Version 2.0 NC_GLOBAL#start_year=1980 NC_GLOBAL#tileid=11752 NC_GLOBAL#Version_data=Daymet Data Version 2.1 NC_GLOBAL#Version_software=Daymet Software Version 2.0 NETCDF_DIM_EXTRA={time} NETCDF_DIM_time_DEF={365,6} NETCDF_DIM_time_VALUES={0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5,12.5 ,13.5,14.5,15.5,16.5,17.5,18.5,19.5,20.5,21.5,22.5,23.5,24.5,25.5,26.5,27.5,28.5 ,29.5,30.5,31.5,32.5,33.5,34.5,35.5,36.5,37.5,38.5,39.5,40.5,41.5,42.5,43.5,44.5 ,45.5,46.5,47.5,48.5,49.5,50.5,51.5,52.5,53.5,54.5,55.5,56.5,57.5,58.5,59.5,60.5
I used gdalwarp as such to try and convert:
GRASS 6.4.4 (demolocation)> gdalwarp -overwrite -s_srs EPSG:9001 -t_srs "+proj= lcc +lat_1=25 +lat_2=60 +lat_0=42.5 +lon_0=-100 +x_0=0 +y_0=0 +datum=WGS84 +unit s=m +no_defs" -tr 10 10 -r bilinear -multi -dstnodata 0 -of GTiff ../dems/Sus_10 m.tif Sus_10mLCC.tif Creating output file that is 2269P x 2147L. Processing input file ../dems/Sus_10m.tif. Using internal nodata values (e.g. nan) for image ../dems/Sus_10m.tif. 0...10...20...30...40...50...60...70...80...90...100 - done.
GRASS 6.4.4 (demolocation)>
the projection looks right, but the coordinates are off. I overlay these on top of a DEM for the entire north american continent and the Susquehanna shall dem ends up somewhere in Idaho.. that's not right.
Anyways, like I said, I'm a complete beginner. The proj4 string is correct, I verified it with my domain scientist. I would ask him to help with this, but he's very busy at the moment. thanks for any help, it's much appreciated.
Aucun commentaire:
Enregistrer un commentaire