Ciao Jose', two things: 1> it would be nice to prse such metadata, with a little bit of help from you ;-) 2> please dont' scale them down, use overviews!
Simone. On Wed, Oct 22, 2008 at 3:14 PM, zé <[EMAIL PROTECTED]> wrote: > > Olá Simone! > > This metadata is given by the Brazilian National Institute of Space > Researches (INPE - www.inpe.br) as attachments of images downloaded from its > remote sensoring program. The images tested comes from LANDSAT5, BAND4, and > are scaled into 10% of originals for speed improvement. > > Thanks for the quick answer! > Best regards. > > > Simone.Giannecchini wrote: >> >> Ciao Jose, >> I am going to do a quick check on your code. >> Iam curios about where this type of metadata comes from. >> >> Simone. >> On Wed, Oct 22, 2008 at 2:38 PM, zé <[EMAIL PROTECTED]> wrote: >>> >>> Hi! >>> >>> I'm trying to plot a same geo point into two differents images, but the >>> regions don't match (see the following images). The images are displaying >>> the same Earth's region in differents envelopes and time, referenced >>> using >>> DATUM SAD69. There're too a XML metadata about each image. >>> >>> http://www.nabble.com/file/p20109846/imagem3.png >>> >>> See that the point in first image is over on the river and it's out in >>> second image. >>> >>> Here, a snippet of the XMLs files: >>> >>> <datumName>SAD69</datumName> >>> <semiMajorAxis>6378160.000000</semiMajorAxis> >>> <flattening>0.003353</flattening> >>> <xShiftToWGS84>-66.870000</xShiftToWGS84> >>> <yShiftToWGS84>4.370000</yShiftToWGS84> >>> <zShiftToWGS84>-38.520000</zShiftToWGS84> >>> <projectionName>UTM</projectionName> >>> <originLongitude>-60</originLongitude> <=== This's the >>> only difference >>> <originLatitude>0</originLatitude> >>> <scaleFactor>0.999600</scaleFactor> >>> <falseEasting>500000.000000</falseEasting> >>> <falseNorthing>10000000.000000</falseNorthing> >>> >>> <datumName>SAD69</datumName> >>> <semiMajorAxis>6378160.000000</semiMajorAxis> >>> <flattening>0.003353</flattening> >>> <xShiftToWGS84>-66.870000</xShiftToWGS84> >>> <yShiftToWGS84>4.370000</yShiftToWGS84> >>> <zShiftToWGS84>-38.520000</zShiftToWGS84> >>> <projectionName>UTM</projectionName> >>> <originLongitude>-57</originLongitude> <=== This's the >>> only difference >>> <originLatitude>0</originLatitude> >>> <scaleFactor>0.999600</scaleFactor> >>> <falseEasting>500000.000000</falseEasting> >>> <falseNorthing>10000000.000000</falseNorthing> >>> >>> Here a snippet of my code (using GeoTools 2.4.4): >>> >>> ... >>> >>> //Read the source images (these images are scaled in 10% of >>> orignals >>> for speed improvement). >>> File path0 = new File("images/20000905_BAND4.tif"); >>> File path1 = new File("images/20080927_BAND4.tif"); >>> GridCoverage2D src0 = (GridCoverage2D) (new >>> GeoTiffReader(path0)).read(null); >>> GridCoverage2D src1 = (GridCoverage2D) (new >>> GeoTiffReader(path1)).read(null); >>> >>> //Calculate the affine transformations to map point between >>> images. >>> AffineTransform at0 = RendererUtilities.worldToScreenTransform( >>> new ReferencedEnvelope( >>> src0.getEnvelope2D(), >>> src0.getCoordinateReferenceSystem() >>> ), >>> new Rectangle( >>> src0.getRenderedImage().getWidth(), >>> src0.getRenderedImage().getHeight() >>> ) >>> ).createInverse(); >>> AffineTransform at1 = RendererUtilities.worldToScreenTransform( >>> new ReferencedEnvelope( >>> src1.getEnvelope2D(), >>> src1.getCoordinateReferenceSystem() >>> ), >>> new Rectangle( >>> src1.getRenderedImage().getWidth(), >>> src1.getRenderedImage().getHeight() >>> ) >>> ); >>> >>> /* Produce the same result as above... So can I use it instead? >>> AffineTransform at0 = (AffineTransform) >>> src0.getGridGeometry().getGridToCRS(); >>> AffineTransform at1 = (AffineTransform) >>> src1.getGridGeometry().getGridToCRS().inverse(); >>> */ >>> >>> //Map the first image's center into second image. >>> Point2D.Double p0 = new Point2D.Double( >>> src0.getRenderedImage().getWidth() / 2d, >>> src0.getRenderedImage().getHeight() / 2d); >>> Point2D.Double p1 = (Point2D.Double) >>> at1.transform( >>> at0.transform( >>> p0, null), >>> null); >>> >>> //Writes the results. The method pointMark creates a 31x31 >>> boundbox >>> //around the points once calculated. >>> GridCoverage2D dest0 = pointMark(src0, p0); >>> GridCoverage2D dest1 = pointMark(src1, p1); >>> (new GeoTiffWriter(new File("images/tmp0.tif"))).write(dest0, >>> null); >>> (new GeoTiffWriter(new File("images/tmp1.tif"))).write(dest1, >>> null); >>> >>> ... >>> >>> Is My code wrong or is it possible the images are incorrectly geo >>> referenced? How can I workaround on this? >>> >>> All files - project, images and metadatas - are zipped into a test case ( >>> http://www.nabble.com/file/p20109846/CoverageTest.zip CoverageTest.zip ). >>> Thanks for any help. >>> >>> José >>> -- >>> View this message in context: >>> http://www.nabble.com/Problem-with-images-world-coordinates-tp20109846p20109846.html >>> Sent from the geotools-gt2-users mailing list archive at Nabble.com. >>> >>> >>> ------------------------------------------------------------------------- >>> This SF.Net email is sponsored by the Moblin Your Move Developer's >>> challenge >>> Build the coolest Linux based applications with Moblin SDK & win great >>> prizes >>> Grand prize is a trip for two to an Open Source event anywhere in the >>> world >>> http://moblin-contest.org/redirect.php?banner_id=100&url=/ >>> _______________________________________________ >>> Geotools-gt2-users mailing list >>> [email protected] >>> https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users >>> >> >> >> >> -- >> ------------------------------------------------------- >> Eng. Simone Giannecchini >> GeoSolutions S.A.S. >> Owner - Software Engineer >> Via Carignoni 51 >> 55041 Camaiore (LU) >> Italy >> >> phone: +39 0584983027 >> fax: +39 0584983027 >> mob: +39 333 8128928 >> >> >> http://www.geo-solutions.it >> http://www.geo-solutions.it/simone.giannecchini >> http://www.linkedin.com/in/simonegiannecchini >> >> ------------------------------------------------------- >> >> ------------------------------------------------------------------------- >> This SF.Net email is sponsored by the Moblin Your Move Developer's >> challenge >> Build the coolest Linux based applications with Moblin SDK & win great >> prizes >> Grand prize is a trip for two to an Open Source event anywhere in the >> world >> http://moblin-contest.org/redirect.php?banner_id=100&url=/ >> _______________________________________________ >> Geotools-gt2-users mailing list >> [email protected] >> https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users >> >> > > -- > View this message in context: > http://www.nabble.com/Problem-with-images-world-coordinates-tp20109846p20110470.html > Sent from the geotools-gt2-users mailing list archive at Nabble.com. > > > ------------------------------------------------------------------------- > This SF.Net email is sponsored by the Moblin Your Move Developer's challenge > Build the coolest Linux based applications with Moblin SDK & win great prizes > Grand prize is a trip for two to an Open Source event anywhere in the world > http://moblin-contest.org/redirect.php?banner_id=100&url=/ > _______________________________________________ > Geotools-gt2-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users > -- ------------------------------------------------------- Eng. Simone Giannecchini GeoSolutions S.A.S. Owner - Software Engineer Via Carignoni 51 55041 Camaiore (LU) Italy phone: +39 0584983027 fax: +39 0584983027 mob: +39 333 8128928 http://www.geo-solutions.it http://www.geo-solutions.it/simone.giannecchini http://www.linkedin.com/in/simonegiannecchini ------------------------------------------------------- ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Geotools-gt2-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
