Le portail francophone de la géomatique

Bonjour a tous,

Pourquoi 60m??? Vous pouvez vous contenter de ne degrader que l'image SPOT
en la reechantillonnat a 30m, non?
Sinon, il existe diverses techniques basees sur les pyramides d'images
ondelettes, pyramide morphologique ...
Un filtrage leger (filtre centre ou autre), suivi d'un echantillonnage
devrait permettre de faire ce que vous souhaitez. Apres, pour des
traitements plus complexes, vous pouvez aller voir du cote de la fusion,
mais les methodes ci dessus sont peu couteuses en temsp de calcul et
fournissent des resultats apreciables.

Cordialement.

Bonjour,

Dans le cadre d'un projet sur la Chine nous avions sous la main une
image Landsat multi spectral 30 m et IRS 5m que nous avons fusionnes
et qui a donne d'excellent resultat. Un article a ete ecrit sur le
sujet , en anglais. Je vous copie le resume ci-dessous.
Si vous desirez l'article complet, vous pouvez le loader depuis mon
site web (sous le lien CV).

Meilleures salutations

Olivier

Photo-Textured Digital Terrain Models as a Basis for Regional and
Local Planning

Li Zhang1, Jacques Feiner2, Olivier Louy2, Armin Gruen1, Willy A. Schmid2

1Institute of Geodesy & Photogrammetry, ETH-Hoenggerberg, CH-8093
Zurich, Switzerland, {zhangl,agruen}@geod.baug.ethz.ch
2Institute for Local, Regional and National Planning,
ETH-Hoenggerberg, CH-8093 Zurich, Switzerland,
{feiner,louy,schmid}@orl.arch.ethz.ch

KEY WORDS: Regional Planning, DTM, Ortho-rectification, Image
Fusion, 3D Visualisation

SUMMARY: Many cities in the second and third world experience
currently strong urbanization and are expanding fast. These
development dynamics are so strong, that the traditional maps as a
basis for planning are often out of date or are- in certain cases
just not available. However, to monitor, control and guide these
development trends, up-to date maps are needed as a basis for urban
and regional planning.
We will show in our contribution how a new form of digital map - a
photo-textured Digital Terrain Model (DTM) with integrated
information for planning measures - can be produced efficiently and
used in the planning process: For the study area of named Greater
Kunming Area we produced a DTM. Then we geocoded and orthorectified
LANDSAT 7 ETM MS and IRS-1C PAN satellite images. In a special
process of image fusion we combined the information of high spectral
information (MS - Multi-spectral with 30 m pixelsize) and high metric
resolution (PAN - Panchromatic, 5 m pixelsize) to a joint image
representation. Planning information is used as a third layer of
data. The combined datasets can be interpreted as a photorealistic 3D
map, which can be accessed and represented by modern visualization
tools. Finally we show how such a 3D digital model can be used for
illustrating planning approaches.

----------------------------------------------------------------
Olivier Louy, M.Sc.
Researcher / GIS Manager

Netzwerk Stadt und Landschaft (NSL)
Institut fur Raum- und Landschaftsentwiclung (IRL-Institute)
Institute for Territorial Development and Landscape
Department of Civil, Environmental and Geomatics Engineering
Swiss Federal Institute of Technology Zurich
ETH Hoenggerberg HIL H28.1
CH - 8093 Zurich
Schweiz / Switzerland

Tel Office : +41 1 633 29 99
Fax :      +41 1 633 11 02
Mobile:    +41 79 378 17 19
e-mail: louy@texte-a-enlever.nsl.ethz.ch
-----------------------------------------------------------------
Web Site: http://www.louy.ch/olivier

Ca depend evidemment des besoins, et de la finesse de modelisation
geometrique pour superposer les images.
Sans faire de traitement d'image complexe, on peut faire une analyse
visuelle en superposant les deux images et en jouant sur la transparence et
les radiometries de chaque image. Des astrophysiciens utilisent cette
visualisation en transparence d'images de mesures physiques ayant des
mailles de resolutions tres differentes pour completer les calculs.

Francois Aumonier

Je crois avoir trouver moi même la résolution à mon pb...
c'était juste que mon image SPOT avait 1 pixel de large de plus que mon image landsat...
merci qd même et bonne aprem
anaïs

Effectivement, elle n'a été que "callée" suivant la projection cartographique standard WGS84... elle n'a pas été orthorectifiée, ils m'ont juste demandé une altitude moyenne sur l'ensemble de la zone couverte. mais je travaille sur une zone littorale désertique très plane... donc est ce vraiment nécesaire de l'orthorectifiée ?

Si j'ai bien compris comment ca fonctionne, pas nécessaire. L'orthorectification est utile en cas de fort reliefs, mais dans votre cas, j'ai l'impression que ce n'est pas utile

Ah... C'est à dire ? Je m'y remet doucement Bernard, éclairez moi de vos lumières...
Question bêête : c'est quoi l'autocorrélation ? Wikipédia me crache plien de formules

Salut robin,

J'ai pas bien compris la citation...............

Voila : t'as deux images, différentes en dates , meme si ce sont theoriquement des dates simultanées.

t'as des differences de positionnements de pixel a pixel. C'est pas tout a fait le meme angle donc pas tout a fait la meme geometrie. Mais tu peux "marier ", "maching" en saxon, chaque pixel d'une image a un pixel de l'autre image par quelques calculs simples ou compliqués (voir la these de laurent Renouard), c'est ce que fait la photogrametrie. C'est ce qu'on appelle "autocorrellation" . Quand tu as tes points, homologues, tu fais une deformation par triangulation de Delaunay. L'image ainsi recalculée est "tres" proche, geometriquement de l'image de reference.

C'est ce que pratique tous les fournisseurs d'images, spot , Ikonos ou autre.  C'est la seule facon que je connaisse de superposer correctement deux images. Mais Emilie a pas encore bien compris.


Formule simple : la valeur absolue de la plus petite difference entre les deux images dans une fenetre 3x3 ou 7X7 donne le point homologue.  En mailles glissantes, evidemment; La taille de la fenetre et son utilisation fait le must des softs actuels.

Dans le cas de Francoise, c'est assez simple. Le problème de la resolution des images n'intervient vraiment que lorsqu'un meme objet est tres deformé. par exemple un batiment vertical et un autre oblique a une echelle qu'on va dire du 1/1000. La c'est plus coton.

Pas sur d'etre clair.

Cordialement

Bonjour,
Me revoici, me revoila...
Merci Robin et Lortic pour vos réponses...
Je pense effectivement devoir faire un calage de mes images... j'ai deux trois pixels de décalage entre mes Landsat (très bien callée entre elles : 1987 et 2000) et ma SPOT (2005).
J'ai par ailleurs deux images SPOT (2005-05 et 2005-12) qui ne sont pas callées du tout.
Lortic, tu parles de "Maching" et "mariage" entre les pixels...
Savez vous comment je dois procéder sous ENVI ??
merci

Dernière modification par Anaïs (Fri 21 December 2007 09:44)

ok... merci !
j'ai essayé. ça marche pas mal.
en fait j'ai beaucoup de mal avec l'aide car je lis très mal l'anglais... d'où toutes mes questions....
en fait j'ai deux images SPOT, une au 10m et l'autre au 20m...
Est ce qu'il vaut mieux faire le calage d'abord ou le rééchantillonage ?
Merci

dac ! merci.
le callage ne serait pas meilleur si je passe en 10m pour les deux images ?

Je ne sais pas en fait pour ta question précedente il faut comparer les RMS des transformations avec en 1>re-échantillonnage avant et en 2>re-échantillonnage après.

En fait tu veux sur-échantillonner l'image à 20m, mais cela ne changera pas les résultats du Map registration, car tu ne peux pas créer de l'information?
C'est pour cela que pour ton problème de calage, moi j'utilisera plutôt une image 10m dégradée à 20m? Et je comparerai les images à 20m.

Après, tu as raison il faut prendre comme image de référence toujours l'image à la plus haute résolution.
Je regarderai cela plus en détails un peu plus tard.
Bon courage.