Comparative study on image data model for integrating multiple resolution image/raster data
Taizou Yamamoto and Ryosuke Shibasaki
Institute of Industrial Science, University of Tokyo
7-22-1, Roppongi, Minato-Ku Tokyo 106, Japan
Tel: (81)-3-3402-6231 Fax: (81)-3-3479-2762
E-mail:yamamoto@shunji.iis.u-tokyo.ac.jp
Abstract
As a wide variety of satellite images with different resolutions become available and various types of geographic information is being accumulated, it is rather common among users to integrate those data for their application s. To conduct consistent and efficient processing of a wide range of geographic information, it is essential to develop data models for geographic information together with efficient operation method. However, we have no sufficient sets of data models and operations for handling image data as geographic or spatial data. This insufficiency leads to the generation of artifacts or the degradation of image data quality in overlying image data or transferring it to different coordinate systems by resampling. This study reveals possible degradation in the overlay and resampling operations, and makes comparison of data models and operations in term of computational load and its degradation. The results implies a possible future image data model.
1.Introduction
As a great variety of image data become available due to the launch of many remote sensing satellites, techniques of extracting necessary information by integrating these data are increasingly important. And as the development of GIS data and remote sensing image data are becoming very common. GIS is now an indispensable system as a platform of integrating information. For GIS, accurate and consistent processing are achieved by representing geographic information using data models and data handling.
So, what kinds of models should be used for imaged data? Image data is intuitively handled as raster data because in practice it is a two dimensional array of pixels which correspond with a certain unit area on the ground surface. But image data may not be raster data exactly all the time since the shape of each pixel neither be square nor rectangular due to viewing angle effects and topographic distortions. This means that we need to devise a new model to represent foot print more accurately, and a set of operations to handle them efficiently.
The objective of this study is to propose and make comparison of image data models together with operations consistent with those models. The performance of those models is evaluated in terms of computational complexity and size of the models. The shape of foot prints may be different, appending on the types of sensors. We paid attention here only to optical sensors whose foot prints are either square or rectangular.
2.Image data models
2.1 Necessity of image data model
Conventionally, pixel position is represented by its central points, while the shape of the pixel is only ambiguously represented by resolution. The shape of the pixel on the ground (foot print) is not represented explicitly. This kind of representation of image data is not sufficient for mapping foot print. We need to represent the foot print more clearly and explicitly. The ambiguity of the foot print
representation cause several serious problems in "resampling" process. Resampling is an operation of projecting an image data onto grid cells defined on a given two dimensional coordinate system (Fig.1). Usually, for the implementation of the resampling, centre position of each grid cell is transferred to the image plane, and pixel value at the transferred location is given with interpolation. In many cases, since the size of the grid cells are set to be almost the same with that of the pixel, each pixel of the image data is mapped to a single grid cell.

Fig.1 Resampling
Fig.2 summarizes the possible errors, which may occur in the resampling process. Error percentage denotes the ratio of unoverlapped portion to the original footprint. Fig.3 shows that the error percentage are strongly affected by the relative size of grid cells, rotation angles between the image coordinate systems (line and pixel directions) and the grid-cell is as large as the pixel, maximum error could be 90%. This kind of errors also affect severely to the quality of product, when we overlay high resolution image data on low resolution image data which are resampled to the same coordinate system of the high resolution image data. This problem comes from the same source a the shape of footprints of pixels are not properly represented and no sufficient algorithms for operating those data models are developed. By applying smaller grid cells, of course, this problem can be solved with a conventional resampling algorithm. In those cases, however, computational load of the resampling would be proportional to the numbers of grid cells, which means the load may become quite heavy if "high quality resampling" is required. As will see later, a data model which represent explicitly the footprint boundary can mitigate the increase of computational load rather effectively.

Fig.2 Degradation of Image data in Resampling.

Fig.3 Average Percentage of the Unoverlapping area by the rotation angle* and grid size.
See* Fig.2