Olá Simone! "Therefore, if you map the center of the first image on the second one you won't get the center of that image marked." Ok! But the point into the second image wouldn't be over the river as the first one?
"I will try to suggest a better snippet tomorrow anyway." I'm really thanks for! :) "1> it would be nice to prse such metadata, with a little bit of help from you ;-)" I'm a newbie in GeoTools and even in geoscience, but I like to help anyway, with a bit community's help... ;) "please dont' scale them down, use overviews!" Ok, I'll search for, thanks for the tip. Best regards. José Simone.Giannecchini wrote: > > 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 > > -- View this message in context: http://www.nabble.com/Problem-with-images-world-coordinates-tp20109846p20110755.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
