CGIAR DEM + gdalwarp — projection to Cartesian is not happeningCGIAR DEM converted to XYZ consists of Nodata onlyHow can I fix badly specified geographic coordinate system with units in arc-seconds?how to georeference my png using gdalWedge of missing data reprojecting to stereographicgdalwarp malfunctioning for GEOS projectionClipping raster takes forever in QGISUnderstanding gdalwarp reprojectionSlope raster seems to come out distorted?Mapserver / QGis and DEM GeoTIFF different nodata values behaviourTransforming Lambert Conformal Conic to EPSG:3857UTM Projection of DEM-data with gdalwarp
How is it possible for user's password to be changed after storage was encrypted? (on OS X, Android)
Lied on resume at previous job
Doomsday-clock for my fantasy planet
Einstein metrics on spheres
Can I find out the caloric content of bread by dehydrating it?
How to manage monthly salary
Can I legally use front facing blue light in the UK?
Does it makes sense to buy a new cycle to learn riding?
Is this food a bread or a loaf?
What does 'script /dev/null' do?
Is Social Media Science Fiction?
How to deal with fear of taking dependencies
I see my dog run
What is the command to reset a PC without deleting any files
Is it wise to focus on putting odd beats on left when playing double bass drums?
extract characters between two commas?
Unbreakable Formation vs. Cry of the Carnarium
How to make payment on the internet without leaving a money trail?
When blogging recipes, how can I support both readers who want the narrative/journey and ones who want the printer-friendly recipe?
How can I add custom success page
"My colleague's body is amazing"
I’m planning on buying a laser printer but concerned about the life cycle of toner in the machine
Need help identifying/translating a plaque in Tangier, Morocco
Why doesn't a const reference extend the life of a temporary object passed via a function?
CGIAR DEM + gdalwarp — projection to Cartesian is not happening
CGIAR DEM converted to XYZ consists of Nodata onlyHow can I fix badly specified geographic coordinate system with units in arc-seconds?how to georeference my png using gdalWedge of missing data reprojecting to stereographicgdalwarp malfunctioning for GEOS projectionClipping raster takes forever in QGISUnderstanding gdalwarp reprojectionSlope raster seems to come out distorted?Mapserver / QGis and DEM GeoTIFF different nodata values behaviourTransforming Lambert Conformal Conic to EPSG:3857UTM Projection of DEM-data with gdalwarp
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I have downloaded a 30-by-30 tile from
http://srtm.csi.cgiar.org/srtmdata/
Then gdalinfo
on the obtained single *.tif
-file would yield
Driver: GTiff/GeoTIFF
Files: cut_n30e000.tif
cut_n30e000.tif.aux.xml
Size is 36024, 36023
Coordinate System is `'
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841)
Lower Left ( -0.0104166, 29.9904167)
Upper Right ( 30.0095841, 60.0095841)
Lower Right ( 30.0095841, 29.9904167)
Center ( 14.9995837, 45.0000004)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
Min=-444.000 Max=4783.000
Minimum=-444.000, Maximum=4783.000, Mean=348.628, StdDev=398.861
NoData Value=-32768
Metadata:
STATISTICS_MAXIMUM=4783
STATISTICS_MEAN=348.62818252965
STATISTICS_MINIMUM=-444
STATISTICS_STDDEV=398.86067773429
So far so good. So I would try to project this onto a 2D plane, as follows:gdalwarp -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m +no_defs' cut_n30e000.tif result_europe.tmerc.tif
.
After this operation, result_europe.tmerc.tif
is supposed to be represented
with respect to a Cartesian system of coordinates.
But I was surprised to find out that the lat/lon system persists:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36024, 36023
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841) ( 15d 0' 0.00"E, 44d 0' 1.94"N)
Lower Left ( -0.0104166, 29.9904167) ( 15d 0' 0.00"E, 44d 0' 0.97"N)
Upper Right ( 30.0095841, 60.0095841) ( 15d 0' 1.35"E, 44d 0' 1.94"N)
Lower Right ( 30.0095841, 29.9904167) ( 15d 0' 1.35"E, 44d 0' 0.97"N)
Center ( 14.9995837, 45.0000004) ( 15d 0' 0.67"E, 44d 0' 1.46"N)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
Why the original lat/lon bounds are still there? I.e. what was wrong with my gdalwarp
command above?
dem gdalwarp
add a comment |
I have downloaded a 30-by-30 tile from
http://srtm.csi.cgiar.org/srtmdata/
Then gdalinfo
on the obtained single *.tif
-file would yield
Driver: GTiff/GeoTIFF
Files: cut_n30e000.tif
cut_n30e000.tif.aux.xml
Size is 36024, 36023
Coordinate System is `'
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841)
Lower Left ( -0.0104166, 29.9904167)
Upper Right ( 30.0095841, 60.0095841)
Lower Right ( 30.0095841, 29.9904167)
Center ( 14.9995837, 45.0000004)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
Min=-444.000 Max=4783.000
Minimum=-444.000, Maximum=4783.000, Mean=348.628, StdDev=398.861
NoData Value=-32768
Metadata:
STATISTICS_MAXIMUM=4783
STATISTICS_MEAN=348.62818252965
STATISTICS_MINIMUM=-444
STATISTICS_STDDEV=398.86067773429
So far so good. So I would try to project this onto a 2D plane, as follows:gdalwarp -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m +no_defs' cut_n30e000.tif result_europe.tmerc.tif
.
After this operation, result_europe.tmerc.tif
is supposed to be represented
with respect to a Cartesian system of coordinates.
But I was surprised to find out that the lat/lon system persists:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36024, 36023
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841) ( 15d 0' 0.00"E, 44d 0' 1.94"N)
Lower Left ( -0.0104166, 29.9904167) ( 15d 0' 0.00"E, 44d 0' 0.97"N)
Upper Right ( 30.0095841, 60.0095841) ( 15d 0' 1.35"E, 44d 0' 1.94"N)
Lower Right ( 30.0095841, 29.9904167) ( 15d 0' 1.35"E, 44d 0' 0.97"N)
Center ( 14.9995837, 45.0000004) ( 15d 0' 0.67"E, 44d 0' 1.46"N)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
Why the original lat/lon bounds are still there? I.e. what was wrong with my gdalwarp
command above?
dem gdalwarp
1
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
1
yes that's cruddy, the files don't have the CRS inside - add-s_srs "+init=epsg:4326"
to your gdalwarp call
– mdsumner
Apr 4 at 5:57
add a comment |
I have downloaded a 30-by-30 tile from
http://srtm.csi.cgiar.org/srtmdata/
Then gdalinfo
on the obtained single *.tif
-file would yield
Driver: GTiff/GeoTIFF
Files: cut_n30e000.tif
cut_n30e000.tif.aux.xml
Size is 36024, 36023
Coordinate System is `'
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841)
Lower Left ( -0.0104166, 29.9904167)
Upper Right ( 30.0095841, 60.0095841)
Lower Right ( 30.0095841, 29.9904167)
Center ( 14.9995837, 45.0000004)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
Min=-444.000 Max=4783.000
Minimum=-444.000, Maximum=4783.000, Mean=348.628, StdDev=398.861
NoData Value=-32768
Metadata:
STATISTICS_MAXIMUM=4783
STATISTICS_MEAN=348.62818252965
STATISTICS_MINIMUM=-444
STATISTICS_STDDEV=398.86067773429
So far so good. So I would try to project this onto a 2D plane, as follows:gdalwarp -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m +no_defs' cut_n30e000.tif result_europe.tmerc.tif
.
After this operation, result_europe.tmerc.tif
is supposed to be represented
with respect to a Cartesian system of coordinates.
But I was surprised to find out that the lat/lon system persists:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36024, 36023
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841) ( 15d 0' 0.00"E, 44d 0' 1.94"N)
Lower Left ( -0.0104166, 29.9904167) ( 15d 0' 0.00"E, 44d 0' 0.97"N)
Upper Right ( 30.0095841, 60.0095841) ( 15d 0' 1.35"E, 44d 0' 1.94"N)
Lower Right ( 30.0095841, 29.9904167) ( 15d 0' 1.35"E, 44d 0' 0.97"N)
Center ( 14.9995837, 45.0000004) ( 15d 0' 0.67"E, 44d 0' 1.46"N)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
Why the original lat/lon bounds are still there? I.e. what was wrong with my gdalwarp
command above?
dem gdalwarp
I have downloaded a 30-by-30 tile from
http://srtm.csi.cgiar.org/srtmdata/
Then gdalinfo
on the obtained single *.tif
-file would yield
Driver: GTiff/GeoTIFF
Files: cut_n30e000.tif
cut_n30e000.tif.aux.xml
Size is 36024, 36023
Coordinate System is `'
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841)
Lower Left ( -0.0104166, 29.9904167)
Upper Right ( 30.0095841, 60.0095841)
Lower Right ( 30.0095841, 29.9904167)
Center ( 14.9995837, 45.0000004)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
Min=-444.000 Max=4783.000
Minimum=-444.000, Maximum=4783.000, Mean=348.628, StdDev=398.861
NoData Value=-32768
Metadata:
STATISTICS_MAXIMUM=4783
STATISTICS_MEAN=348.62818252965
STATISTICS_MINIMUM=-444
STATISTICS_STDDEV=398.86067773429
So far so good. So I would try to project this onto a 2D plane, as follows:gdalwarp -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m +no_defs' cut_n30e000.tif result_europe.tmerc.tif
.
After this operation, result_europe.tmerc.tif
is supposed to be represented
with respect to a Cartesian system of coordinates.
But I was surprised to find out that the lat/lon system persists:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36024, 36023
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-0.010416618510297,60.009584059976817)
Pixel Size = (0.000833333353512,-0.000833333353512)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -0.0104166, 60.0095841) ( 15d 0' 0.00"E, 44d 0' 1.94"N)
Lower Left ( -0.0104166, 29.9904167) ( 15d 0' 0.00"E, 44d 0' 0.97"N)
Upper Right ( 30.0095841, 60.0095841) ( 15d 0' 1.35"E, 44d 0' 1.94"N)
Lower Right ( 30.0095841, 29.9904167) ( 15d 0' 1.35"E, 44d 0' 0.97"N)
Center ( 14.9995837, 45.0000004) ( 15d 0' 0.67"E, 44d 0' 1.46"N)
Band 1 Block=36024x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
Why the original lat/lon bounds are still there? I.e. what was wrong with my gdalwarp
command above?
dem gdalwarp
dem gdalwarp
asked Apr 4 at 4:43
IlonpilaajaIlonpilaaja
1366
1366
1
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
1
yes that's cruddy, the files don't have the CRS inside - add-s_srs "+init=epsg:4326"
to your gdalwarp call
– mdsumner
Apr 4 at 5:57
add a comment |
1
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
1
yes that's cruddy, the files don't have the CRS inside - add-s_srs "+init=epsg:4326"
to your gdalwarp call
– mdsumner
Apr 4 at 5:57
1
1
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
1
1
yes that's cruddy, the files don't have the CRS inside - add
-s_srs "+init=epsg:4326"
to your gdalwarp call– mdsumner
Apr 4 at 5:57
yes that's cruddy, the files don't have the CRS inside - add
-s_srs "+init=epsg:4326"
to your gdalwarp call– mdsumner
Apr 4 at 5:57
add a comment |
1 Answer
1
active
oldest
votes
Just re-stating the answers provided by @MichaelStimson and @mdsummer:gdalwarp -s_srs '+init=epsg:4326' -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m' cut_n30e000.tif result_europe.tmerc.tif
will do the job.
After this, gdalinfo
gives the proper:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36696, 43220
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-1456758.150925980880857,1876581.142161503667012)
Pixel Size = (79.393568841911133,-79.393568841911133)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1456758.151, 1876581.142) ( 10d15'41.62"W, 58d20'50.81"N)
Lower Left (-1456758.151,-1554808.903) ( 0d 7' 0.41"E, 29d 8'47.69"N)
Upper Right ( 1456668.251, 1876581.142) ( 40d15'36.63"E, 58d20'51.87"N)
Lower Right ( 1456668.251,-1554808.903) ( 29d52'56.37"E, 29d 8'48.06"N)
Center ( -44.950, 160886.119) ( 14d59'57.93"E, 45d26'51.99"N)
Band 1 Block=36696x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
add a comment |
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317722%2fcgiar-dem-gdalwarp-projection-to-cartesian-is-not-happening%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
1 Answer
1
active
oldest
votes
1 Answer
1
active
oldest
votes
active
oldest
votes
active
oldest
votes
Just re-stating the answers provided by @MichaelStimson and @mdsummer:gdalwarp -s_srs '+init=epsg:4326' -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m' cut_n30e000.tif result_europe.tmerc.tif
will do the job.
After this, gdalinfo
gives the proper:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36696, 43220
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-1456758.150925980880857,1876581.142161503667012)
Pixel Size = (79.393568841911133,-79.393568841911133)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1456758.151, 1876581.142) ( 10d15'41.62"W, 58d20'50.81"N)
Lower Left (-1456758.151,-1554808.903) ( 0d 7' 0.41"E, 29d 8'47.69"N)
Upper Right ( 1456668.251, 1876581.142) ( 40d15'36.63"E, 58d20'51.87"N)
Lower Right ( 1456668.251,-1554808.903) ( 29d52'56.37"E, 29d 8'48.06"N)
Center ( -44.950, 160886.119) ( 14d59'57.93"E, 45d26'51.99"N)
Band 1 Block=36696x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
add a comment |
Just re-stating the answers provided by @MichaelStimson and @mdsummer:gdalwarp -s_srs '+init=epsg:4326' -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m' cut_n30e000.tif result_europe.tmerc.tif
will do the job.
After this, gdalinfo
gives the proper:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36696, 43220
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-1456758.150925980880857,1876581.142161503667012)
Pixel Size = (79.393568841911133,-79.393568841911133)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1456758.151, 1876581.142) ( 10d15'41.62"W, 58d20'50.81"N)
Lower Left (-1456758.151,-1554808.903) ( 0d 7' 0.41"E, 29d 8'47.69"N)
Upper Right ( 1456668.251, 1876581.142) ( 40d15'36.63"E, 58d20'51.87"N)
Lower Right ( 1456668.251,-1554808.903) ( 29d52'56.37"E, 29d 8'48.06"N)
Center ( -44.950, 160886.119) ( 14d59'57.93"E, 45d26'51.99"N)
Band 1 Block=36696x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
add a comment |
Just re-stating the answers provided by @MichaelStimson and @mdsummer:gdalwarp -s_srs '+init=epsg:4326' -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m' cut_n30e000.tif result_europe.tmerc.tif
will do the job.
After this, gdalinfo
gives the proper:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36696, 43220
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-1456758.150925980880857,1876581.142161503667012)
Pixel Size = (79.393568841911133,-79.393568841911133)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1456758.151, 1876581.142) ( 10d15'41.62"W, 58d20'50.81"N)
Lower Left (-1456758.151,-1554808.903) ( 0d 7' 0.41"E, 29d 8'47.69"N)
Upper Right ( 1456668.251, 1876581.142) ( 40d15'36.63"E, 58d20'51.87"N)
Lower Right ( 1456668.251,-1554808.903) ( 29d52'56.37"E, 29d 8'48.06"N)
Center ( -44.950, 160886.119) ( 14d59'57.93"E, 45d26'51.99"N)
Band 1 Block=36696x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
Just re-stating the answers provided by @MichaelStimson and @mdsummer:gdalwarp -s_srs '+init=epsg:4326' -t_srs '+proj=tmerc +lat_0=44.00 +lon_0=15.00 +units=m' cut_n30e000.tif result_europe.tmerc.tif
will do the job.
After this, gdalinfo
gives the proper:
Driver: GTiff/GeoTIFF
Files: result_europe.tmerc.tif
Size is 36696, 43220
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["unknown",
SPHEROID["WGS84",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",44],
PARAMETER["central_meridian",15],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-1456758.150925980880857,1876581.142161503667012)
Pixel Size = (79.393568841911133,-79.393568841911133)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1456758.151, 1876581.142) ( 10d15'41.62"W, 58d20'50.81"N)
Lower Left (-1456758.151,-1554808.903) ( 0d 7' 0.41"E, 29d 8'47.69"N)
Upper Right ( 1456668.251, 1876581.142) ( 40d15'36.63"E, 58d20'51.87"N)
Lower Right ( 1456668.251,-1554808.903) ( 29d52'56.37"E, 29d 8'48.06"N)
Center ( -44.950, 160886.119) ( 14d59'57.93"E, 45d26'51.99"N)
Band 1 Block=36696x1 Type=Int16, ColorInterp=Gray
NoData Value=-32768
answered Apr 4 at 11:26
IlonpilaajaIlonpilaaja
1366
1366
add a comment |
add a comment |
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317722%2fcgiar-dem-gdalwarp-projection-to-cartesian-is-not-happening%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
1
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
Apr 4 at 4:47
1
yes that's cruddy, the files don't have the CRS inside - add
-s_srs "+init=epsg:4326"
to your gdalwarp call– mdsumner
Apr 4 at 5:57