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













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?









share























  • What are you specifying as -s_srs? It's not defined in your from info.

    – Michael Stimson
    1 min ago
















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?









share























  • What are you specifying as -s_srs? It's not defined in your from info.

    – Michael Stimson
    1 min ago














0












0








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?









share














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





share












share










share



share










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



















  • 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










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
});


}
});














draft saved

draft discarded


















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
















draft saved

draft discarded




















































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.




draft saved


draft discarded














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





















































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







Popular posts from this blog

Щит и меч (фильм) Содержание Названия серий | Сюжет |...

is 'sed' thread safeWhat should someone know about using Python scripts in the shell?Nexenta bash script uses...

Meter-Bus Содержание Параметры шины | Стандартизация |...