CGIAR DEM + gdalwarp — projection to Cartesian is not happeningHow can I fix badly specified geographic...
Plagiarism or not?
How to tell a function to use the default argument values?
How does a predictive coding aid in lossless compression?
GFCI outlets - can they be repaired? Are they really needed at the end of a circuit?
Can a virus destroy the BIOS of a modern computer?
What exploit Are these user agents trying to use?
Alternative to sending password over mail?
Is it acceptable for a professor to tell male students to not think that they are smarter than female students?
Do scales need to be in alphabetical order?
I would say: "You are another teacher", but she is a woman and I am a man
Reverse dictionary where values are lists
Arrow those variables!
Why is consensus so controversial in Britain?
Venezuelan girlfriend wants to travel the USA to be with me. What is the process?
CAST throwing error when run in stored procedure but not when run as raw query
In 'Revenger,' what does 'cove' come from?
Why are the 737's rear doors unusable in a water landing?
Short story with a alien planet, government officials must wear exploding medallions
Expand and Contract
Examples of smooth manifolds admitting inbetween one and a continuum of complex structures
Why was the shrinking from 8″ made only to 5.25″ and not smaller (4″ or less)?
Are there any examples of a variable being normally distributed that is *not* due to the Central Limit Theorem?
Size of subfigure fitting its content (tikzpicture)
One verb to replace 'be a member of' a club
CGIAR DEM + gdalwarp — projection to Cartesian is not happening
How 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
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
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago
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 4 mins ago
IlonpilaajaIlonpilaaja
1195
1195
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago
add a comment |
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago
add a comment |
0
active
oldest
votes
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
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
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
What are you specifying as -s_srs? It's not defined in your from info.
– Michael Stimson
1 min ago