Vol 101
Vol 101
Journal of Geomatics
(A publication of the Indian Society of Geomatics)
Editorial Board
Chief Editor: Prof. Ajai
(Address for Correspondence: Space Applications Centre, ISRO, Ahmedabad - 380 015)
Phone: +91-79-26913305 (O), 91-02717-235441 (R), Email: [email protected]
Associate Editors:
R. P. Dubey Ahmedabad, Email: [email protected]
Markand P. Oza Ahmedabad, Email: [email protected]
Members:
V. Balaji Hyderabad, Email: [email protected]
Mahesh Chandra New Delhi, Email: [email protected]
A.R. Dasgupta Ahmedabad, Email: [email protected]
P.K. Garg Dehradun, Email: [email protected]
A.K. Gosain New Delhi, Email: [email protected]
Ashok Kaushal Pune, Email: [email protected]
I.V. Murali Krishna Hyderabad, Email: [email protected]
V.N. Patkar Mumbai, Email: [email protected]
S.M. Ramasamy Tiruchirapalli, Email: [email protected]
Aniruddha Roy New Delhi, Email: [email protected]
P.S. Roy Hyderabad, Email: [email protected]
Milap Chand Sharma New Delhi, Email: [email protected]
P. Venkatachalam Mumbai, Email: [email protected]
Advisory Board
Paul J. Curran Vice-Chancellor, Bournemouth University, Poole, UK
V. Jayaraman Bengaluru, India
R. Krishnan Thiruvananthpuram, India
Sugata Mitra NIIT GIS Ltd, New Delhi, India
P. Nag Varanasi, India
M.P. Narayanan President, CSDMS, NOIDA, U.P., India
R.R. Navalgund ISRO H.Q., Bengaluru India
Y.S. Rajan ISRO H.Q., Bengaluru, India
R. Siva Kumar New Delhi, India
Josef Strobl Interfaculty Dept. of Geoinformatics, University of Salzburg, Austria
ii
Ex-Officio (Immediate Past President) Shailesh Nayak, MoES, New Delhi – 110003
Journal of Geomatics
(A publication of the Indian Society of Geomatics)
Vol. 10 No. 1 April 2016
Research articles
1
1 Flood hazard zoning using analytic hierarchy process: A case study for
Pampa river basin, Kerala, India
N.A. Mayaja and C.V. Srinivasa
5 Unfolding the time relationship of structural events through Landsat data: A case study 24
from Khandia formation, Champaner group, Gujarat
M.A. Limaye
6 Runoff estimation from a tributary of lower Tapi basin using SCS-CN method integrated 29
with remote sensing and GIS data
Sudhakar B. Sharma, Anupam K. Singh and Ajay S. Rajawat
7 Change detection on the Volta river due to the construction of the Kpong dam using 36
remote sensing techniques
E.M. Osei Jnr, A.S. Amoah, B.E.K. Prah, F.K. Anyah and K.A. Addo-Gyan
8 An objective method for detecting night time fog using MODIS data over northern India 40
Sasmita Chaurasia and B.S. Gohil
9 Morphometric changes of the Varuna river basin, Varanasi district, Uttar Pradesh 48
K. Prakash, S. Singh and U.K. Shukla
12 Evaluation of predictive ability of support vector machines and naive Bayes trees methods 71
for spatial prediction of landslides in Uttarakhand state (India) using GIS
Binh Thai Pham, DieuTien Bui, Indra Prakash and M.B. Dholakia
iv
14 Formosat-2 with Landsat-8 temporal - multispectral data for wheat crop identification 89
using Hypertangent Kernel based Possibilistic classifier
Rohit Nandan, Ankit Kamboj, Anil Kumar, A. Senthil Kumar and K. Venkata Reddy
Published biannually by the Secretary, Indian Society of Geomatics on behalf of the Society
Copyright Indian Society of Geomatics
ISG Website: www.isgindia.org
1
Abstract: The river Pampa is the third largest river of Kerala with a catchment area of 2235 km2. This river almost every
year causes substantial damage to human life, properties and the cropland during monsoon. In this study an attempt is
made to classify the regions in the river basin in order of risk and severity due to floods. The severity of flood hazard in
these locations varies due to various geospatial factors. The hazard due to flood in any particular location and its impact
can be assessed in relative terms by using an analytical approach as applied to a set of geospatial factors ranging from
qualitative to quantitative type. This paper evolves appropriate risk indices for the entire Pampa river basin and classifies
them according to the severity of flood risk using a popular Analytic Hierarchy Process (AHP). The study brought out
that two regions in the river basin fall under very high flood risk category whereas four villages come under high risk
category. It was revealed that highly populated and urbanized regions located in the downstream of this river basin are
more vulnerable to flood hazard.
Key Words: Flood Hazard Index (FHI), Analytic Hierarchy Process (AHP), Pampa river basin, Kerala floods
6. Conclusions
Figure 3: Flood risk map of Pampa river basin
Ranking the villages in the flood plain is of utmost
5. Results and discussions importance in flood management planning. In this
study, 52 villages in the Pampa river basin are classified
In this study, a set of composite flood hazard indices for according to their Flood Hazard using Analytic
the Pampa river basin has been worked out by adopting Hierarchy Process. The indices have been derived from
the AHP methodology. The indices have been derived a variety of parameters (factors) ranging from geospatial
from primary decision factors viz: population density, data to population density, rainfall, land use, type of
annual average rainfall, land use, type of soil, slope, soil, slope, geomorphic factors, quality of road and
geomorphic features, quality of roads and elevation. The elevation. The flood prone areas of Pampa river basin
Panchayaths/ municipalities of Pampa river basin have have been classified into four categories viz. low,
been accordingly classified into low, medium, high and medium, high and very high. The analysis reveals that
very high risk categories based on histogram human activities, which result in increased population
distribution of FHI. density, land use land cover changes etc. make the
region more vulnerable to flood hazards. Hence a
The analysis revealed two regions - Thiruvanvandoor, comprehensive basin planning, considering the above
and Veeyapuram - covering about 7.5 km2 of the basin factors only will be effective in mitigating flood hazard
as areas prone to ‘very high’ levels of flood risk. Four
regions, covering about 35 km2 come under the ‘high References
risk’ category of flood where as another 30 regions (395
km2) fall under the ‘medium risk’ category. Remaining Mayaja, N.A. and C.V. Srinivasa (2012). A GIS based
portions of the basin are relatively under ‘low risk’. It flood mapping study for Pampa river basin, Kerala.
is observed that though the regions falling under very Research dissertation of Visweswaraya Technological
high and high level of flood hazard constitute only 0.3% University, Karnataka.
and 2 % respectively of the basin area, these are densely
populated and highly urbanised regions (with density of Saaty, T.L. and J. Alexander (1989). Conflict
population more than 3000 persons / km2), located at the resolution: The analytic hierarchy process. Praeger,
downstream of the river. Further, the land use pattern of New York.
these regions reveals high level of built up area and they
have good network of paved highways. Thus, it can be Saaty, T.L. (1994). Fundamentals of decision making
inferred that the prime reasons of flood hazard are high and priority theory with the analytic hierarchy process.
rate of urbanization and human interventions in this ISBN 0-9620317-6-3.
region. The extensive road networks recently developed
in the river basin also testify this finding. Owing to the Saaty, T.L. and L.G. Vargas (2002). Models, methods,
same reasons the flood occurrence at these regions concepts & applications of the analytic hierarchy
causes more damages to both humans as well as process. Book. ISBN 0-7923-7267-0, Kluwer
infrastructure. Academic.
Amongst the eight primary decision factors influential Venkata Bapalu, G. and R. Sinha (2014). GIS in flood
in causing flood hazard, the most prominent hazard mapping: A case study of Kosi river basin,India.
anthropogenic factor identified is population density. download from internet, http://GISdevelopment.net
This calls for urgent need of an effective urban planning
5
Appendix-1
(FHI computation: Decision Hierarchy at levels 1, 2 and 3)
Appendix 2
(The comparison matrix for level- 2)
Appendix 3
(FHI computation)
The FHI for each location was determined by aggregating RIWs at each level of the hierarchy. FHI was calculated by
multiplying the RIWs of level 3-decision factor by the associated RIWs of the level 2 factors at each level and summing
the values of all grouped elements. The level 2 and level 3 Relative Importance Weight matrices computed are shown
below:
Appendix 4
(The consistency ratios at levels 2 and 3)
ISG Website
(http://www/isgindia.org)
The web site of Indian Society of Geomatics contains all pertinent information about
ISG and its activities. The latest announcements can be found on homepage itself.
“About ISG” link gives information about the constitution of ISG and its role in
Geomatics, both the technology and its applications in the Indian context. The site also
furnishes information about the members in different categories, like – Patron
Members, Sustaining Members, Life Members and Annual Members. One can
download Membership form from this section or through the Downloads link. The
website also has full information about the Executive Council Meetings’ Agenda of
past and present along with Executive Agenda and Minutes. The details of local
Chapters’ office bearers are also provided. The Annual General-body Meeting (AGM)
Agenda, their minutes for particular year can also be seen in the “AGM” section. The
list of Events organized by the society can be found through the “Events” link.
Website related queries, suggestions and feedback to improve the website can be sent to
the webmaster at e-mail: [email protected] or [email protected]
7
Abstract: Dimensionality reduction (DR) techniques help in reducing the volume of the hyperspectral data with
minimum loss of information. The three most commonly used DR techniques – PCA, MNF and ICA have their own
advantages and limitations in transforming the redundant hyperspectral data to non-redundant data thus aiding in an
improved feature extraction. In the present work, an attempt was made to analyze the performance of hybrid
dimensionality reduction method which uses a combination of three non – linear DR techniques for extracting the
concrete materials from the CHRIS hyperspectral data. SAM and SID classifiers were used for classifying six different
surface materials (concrete, paved and unpaved) in the study area along with four different vegetation types. Analysis
has shown that the hybrid method gave satisfactory results for classifying the surface materials in CHRIS data. The SAM
classifier gave the best results with an accuracy improvement of 10% after adapting the hybrid method. The classification
accuracies have increased from 79.54% to 85.14% for SAM classification and 80.24% to 84.90% for SID classification.
Keywords: Principle component analysis(PCA), Minimum noise function (MNF), Independent component analysis
(ICA), Hybrid method, Classification
Table 3: Classification accuracies of various surface SAM and SID classifiers had equal performance in
materials considered in classification classification. A SAM angle of 0.15 was used for
classification of both the original and transformed images.
Spectral Angle Spectral The overall average accuracies (Fig. 5) have increased
Mapper Informed from 79.54% to 85.14% for SAM classification and
Divergence 80.24% to 84.90% for SID classification.
Class Original Stacked Original Stacked
name CHRIS Vector CHRIS Vector 86.00% 85.14% 84.90%
New 78.56% 85.25% 79.56% 87.67% 84.00%
Concrete 82.00% 80.24%
79.54%
Tar 85.45% 87.63% 84.68% 88.65% 80.00%
White 88.56% 91.58% 89.25% 91.44% 78.00%
Paint 76.00%
Cemented 72.66% 85.86% 73.58% 80.25%
area
Sandy 74.56% 84.78% 81.21% 82.56%
surface
Old 80.45% 87.45% 81.56% 87.24%
Concrete Figure 5: Overall classification accuracies
Scrubs 81.25% 83.15% 79.54% 81.45%
Thomy 75.45% 77.63% 74.68% 78.65% The spatial resolution of the hyperspectral image plays a
trees major role in the image classification, apart from its high
Bushy 79.21% 83.47% 80.12% 84.65% spectral resolution. In the present work, the spatial
trees resolution of 17m of the CHRIS hyperspectral image is
Croplands 79.25% 84.56% 78.21% 86.46% still coarse to misclassify the pixel into a different class
Average 79.54% 85.14% 80.24% 84.90% (mixed pixel effect).
accuracy
Accuracy assessment was performed over the four Also it is accepted that the transformed image could not
classified images - SAM classified original and completely restore the spectral resolution of the original
transformed stacked images, SID classified original and CHRIS image after using the DR techniques. Hence, a
transformed stacked images. The accuracies for each class compromise in the overall accuracy rate of the classifiers
in all the four classified maps are shown in table 3. Both is accepted.
11
Figure 6: A - Improvements in the extraction of road paths up on classification; B - Improvements in the extraction
of urban and concrete buildings in a township; C- Extraction of curved path; D - Extraction of building roof coated
with artificial paints; E - Google Earth images of figures – A, B, C, D respectively
The authors are thankful to Dr Sateesh Reddy, the Tutorial, E.Z. (2010). ENVI user guide. Colorado Springs,
outstanding scientist and Director, Dr BV Rao, CO: ITT.
Technology Director, Dr Arindan Biswas, Scientist from
IIRS and Rama Sarma, DOMS, RCI, DRDO for their Wang, Y., G. Wu and L. Ding (2014). Plant Species
interest and for providing the necessary infrastructure and Identification Based on Independent Component Analysis
extending all support whenever required. for Hyperspectral Data. Journal of Software 9.6: 1532-
1537.
References
Yusuf, B.L. and Y. He (2011). Application of
Galal, A. and H. Hasan (2012). Learning Flexible hyperspectral imaging sensor to differentiate between the
Hyperspectral Features. International Journal of Remote moisture and reflectance of healthy and infected tobacco
Sensing Applications. IJRSA Vol.2 Iss. 2, 44-48 leaves. African Journal of Agricultural Research, 6(29),
6267-6280.
Green, A.A., M. Berman, P. Switzer and M.D. Craig
(1988). A transformation for ordering multispectral data in
12
Abstract: Precise Point Positioning (PPP) has been proved by many researchers in the last decade as a cost-effective
alternative for Differential GPS (DGPS) with an estimated precision sufficient for many applications. PPP
implementation needs state-of-art software to correct GNSS observations from different types of errors. Several online
PPP services have been developed recently by government agencies, universities, industries and individuals. The PPP
software centre managed by University of New Brunswick (UNB), Canada is offering the user four online PPP services
(CSRS-PPP, GAPS, APPS and magicGNSS) (UNB-PPP, 2015). This research presents an accuracy assessment
evaluation study for those services by processing 3h 52min. dual frequency-static GPS observations which were divided
into 10 sessions with different observation duration (10min, 20min, 30min, 45min., 1 hr, 1.5 hr, 2 hr, 2.5 hr, 3 hr and 3hr
52min.)
Precise point positioning (PPP) is an enhanced single The basic observable for PPP is the ionosphere-delay-
point positioning technique for code or phase free (up to first order effects) pseudo-range and carrier-
measurements using precise orbits and clocks instead of phase observations. The simplified observation
broadcast data. PPP became viable with the existence of equations are given as:
the extremely precise ephemerides and clock
corrections, offered by different organizations.To
compensate for ionospheric effects (the largest source
of error for GPS observations), dual frequency and
measurements are used for an ionosphere free
combination (Rizos et al., 2012).
Reference ITRF2008
system
Coordinate LLH/XYZ
format
Satellite orbit IGS IGS JPL IGS
and clock
ephemeris
Satellite phase IGS ANTEX
centre offsets
Receiver phase IGS ANTEX
centre offsets
Tropospheric
model
Dry model Davis (GPT) UNB-VMF1 Standard A batch least squares
(CMC) formula algorithm that minimizes
Wet model Hopfield model Gradient Chen Standard(0.10m) measurement residuals
(GPT) and Herring Gradient model solving for orbits, satellite
model and station clock offsets,
phase ambiguities and
station tropospheric zenith
delays.
APPS: Automatic Precise Point Service; CSRS: Canadian Spatial Reference System; GAPS: GPS Analysis and
Positioning Software; PPP: precise point positioning; IGS: International GNSS Service; JPL: Jet Propulsion Laboratory;
GPT: global pressure and temperature data; GMF: global mapping function; FES: Finite Element Solution;
(CMC):Canadian Meteorological Centre ;(GPS)Global Positioning System
15
Table 4: The (estimated coordinates – reference coordinates) from the four PPP-services for the 10
observation batches
10 min. observation batch
Estimated coordinate – Reference coordinate Latitude (m) Longitude (m) Height (m)
CSRS-PPP 0.042 0.051 0.071
GAPS 0.029 0.021 0.431
magicGNSS No data No data No data
APPS No data No data No data
20 min. observation batch
CSRS-PPP 0.033 -0.093 -0.255
GAPS 0.049 -0.099 0.031
magicGNSS 0.366 0.375 0.654
APPS 0.087 -0.275 -0.121
30 min. observation batch
CSRS-PPP 0.003 -0.069 -0.144
GAPS 0.044 -0.090 0.001
magicGNSS 0.177 0.126 1.592
APPS 0.075 -0.232 0.034
45 min. observation batch
CSRS-PPP -0.003 -0.012 0.002
GAPS 0.017 -0.013 0.153
magicGNSS 0.108 0.108 1.927
APPS 0.037 -0.085 0.041
1 hr. observation batch
CSRS-PPP 0 0.024 0.028
GAPS 0.014 0.021 0.211
magicGNSS -0.015 0.084 2.267
APPS 0.007 0.042 0.164
1.5 hr. observation batch
CSRS-PPP 0.009 0 0.008
GAPS No data No data No data
magicGNSS 0.039 -0.036 1.879
APPS 0.001 -0.008 0.252
2 hr. observation batch
CSRS-PPP 0.012 -0.003 -0.009
GAPS 0.023 -0.010 0.227
magicGNSS 0.024 -0.012 2.030
APPS 0.001 -0.011 0.251
2.5 hr. observation batch
CSRS-PPP 0.006 -0.003 0.013
GAPS No data No data No data
magicGNSS -0.003 -0.024 2.081
APPS 0.006 -0.011 0.268
3 hr. observation batch
CSRS-PPP 0.003 0 0.011
GAPS 0.014 -0.008 0.264
magicGNSS 0.009 -0.003 2.082
APPS 0.004 -0.012 0.245
3 hr. 52min observation batch
CSRS-PPP 0 0 0
GAPS No data No data No data
magicGNSS 0.009 -0.009 2.038
APPS 0.004 -0.006 0.229
Figure 1: The (estimated coordinates – reference coordinates) in meters from the four PPP-services for the 10
observation batches
The quality of observation data was checked using
4. Test Study UNAVCO’s translation, editing, and quality check
utility (TEQC) software (TEQC, 2012). This set of
To compare the accuracy of these four online PPP observations was divided using TEQC software into 10
services, an observation set of 3h 52 min. dual batches with different observation duration (10min,
frequency-static GPS observations was collected (GPS 20min, 30min, 45min., 1 hr, 1.5 hr, 2 hr, 2.5 hr, 3 hr and
day 17191) with Topcon GR-3 dual frequency receiver 3hr 52min.). The PPP coordinates for each session were
(Topcon GR-3, 2013) using 15 sec observation interval estimated using the four PPP services and compared
and 10o cut-off elevation angle. with the reference coordinates for the tested station. The
reference coordinates is obtained from processing the (3
hr 52 min.) batch using CSRS-PPP service.
18
It can be concluded from this research that different PPP Rizos, C., V. Janssen, C. Roberts and T. Grinter (2012).
services provide different solutions according to Precise point positioning: is the era of differential GNSS
processing strategy. The user should compare his positioning drawing to an end?. FIG Working Week
solutions from number of services and not depend on 2012. Rome, Italy, 6–10 May 2012.
one service. The difference between PPP estimates from
different services could reach 2 meters for height TEQC (2015). TEQC-UNAVCO tutorial.
coordinate and at a decimeter level for latitude and http://facility.unavco.org/software/teqc/doc/UNAVCO
longitude coordinates. The study recommends the _Teqc_ Tutorial.pdf. (Accessed 8/8/2015).
CSRS and APPS services to be used by PPP-users.
Topcon GR-3 (2015). Topcon Positioning Systems
References (GR-3). http://www.topcon.co.jp/en/positioning/
products/pdf/gr-3_e.pdf.(Accessed 8/8/2015).
Collins, P. (2008). Isolating and estimating
undifferenced GPS integer ambiguities. In: Proceedings UNB-PPP (2015). The precise point positioning service.
of ION NTM 2008, San Diego, 28–30 January 2008, pp University of New Brunswick (UNB)
720–732 http://www2.unb.ca/gge/Resources/PPP/. (Accessed
8/8/2015.
19
Abstract: Land is a vital asset for any federal government. Valuation of land is a very important activity as it has high
impact on revenue generation for any country. Present valuation technique is time consuming and cannot quantify the
spatial importance of the land for decision making related to real estate. GIS provides a technical platform for
management of geographic data and inherent location information to support application of spatial statistical and location
econometric tools. Spatial database helps to take decisions for their projects and act as a base data for further functions
like tax calculation, land purchase etc. The main objective of this study is to develop a Real Estate Information System
(REIS) for the valuation purpose especially for the buying, selling and taxation of land properties. This computerized
standalone application, as a repository of land value data is capable to providing easy access to the user/customer about
the land information. The technology used for development is open source and hence offers easy modification and
customization at user end also. The current version is developed for windows environment and can be used and installed
in any windows based system with very less effort. The system can be utilized by the users interested in real estate market
for the land procurement process. The system can also assist government agency for taxation purpose.
Keywords: Geographical information system (GIS), Real Estate Information System (REIS), Open source, Land
valuation
Table-2: Tools and Technologies used in the study help of the scanned and georeferenced wardmap. After
S. No. Component Specification digitization, an attribute table was prepared for the
Operating Windows (32 Bit wardmap layer with the attributes ward_no,
1 system &64 Bit) shape_length, shape_area and ward_name.
2 GIS package Map window
Programming
3 language VB.net
4 Database File based(.shp)
.
Map window active
5 Libraries control(.ocx)
(a) (b)
The data available for ward boundaries and road
network was available in image form. These data are
converted in .shp files by digitization. After base map,
attribute coding was done which includes ward name,
ward no., area for wards and road width, road length,
road name for roads. This work was done in open source
GIS and final base map was prepared. The available data
(c) (d)
was cadastral map which contains Khasra no. area, P. H.
No. and village name. This map was further improvised
Figure 2: (a) Digitized ward map; (b) attribute table
by adding unit circle land rates applicable for various
of ward map; (c) digitized roads map; and (d)
areas as per the guidelines for circle rates issued by
attribute table of road map
administration of that area. Multiplicative factors are
worked out for each of the parameter which affect the
5.2 Preparation of ward wise land value map
land unit rate and hence added to the base map file. The
Ward map so prepared is used to assign the individual
parameters and multiplicative factors as listed below:
cadastral units to different ward the land values as per
the Govt. guidelines and other attributes are coded in the
Landuse
attribute table with the help of Field Calculator tool.
Land Comme Indus Agricul Reside Field names in the value map attribute table are
use rcial trial tural ntial Khasra_no, type, ward_name, area, land_value and
MF 1.50 1.50 1.00 1.00 PH_HALKA.
Depth from road
Depth from Less than 20 More than 20 5.3 Finalizing parameters for valuation and
road Meter meter valuation of sample properties
MF 1.00 0.75 The parameters affecting the value of any property are
Type of road taken from the “Valuation Guidelines of Immovable
Road National State Other Properties” issued by the District Collector, Bhopal.
type Highway Highway Roads These guidelines describe various parameters such as
MF 2.00 1.50 1.20 depth from road, roadtype, land use etc, which are inputs
for working out valuation. Every parameter has been
5. Results and discussion given a weight and these weights are filled across the
polygons in the attribute table of cadastral maps.
The development of Real Estate Information System i.e.
REIS using geospatial techniques has been carried out
in the following stages-
Creation of geospatial database
Valuation of properties
Development of software solution.
Figure 2 shows the ward map and road map with the
22
Figure 4: A sample land value map Figure 6: Graphical User Interface (GUI) of REIS
5.6. Discussion
The application (REIS) is developed for realizing an
analytical system to study the real estate market by the
user. It provides results based on the input and also can
be used for finding costliest or cheapest properties in
any particular area or in the city. Any kind of analysis
can be done like comparison of values in same area or
comparison of same extent of land in different wards.
Figure 5: Valuation of properties This study has provided the GIS based system to work
out the valuation of open land parcels but further
5.4 GUI designing and coding of basic structure of parameters and multiplication factors can be added for
REIS built up parcels owned by individuals so that entire
To fulfil the third objective, a software application has valuation will become more scientific and customised.
been developed using open source technology that Creation of spatial database is very useful for valuation
facilitates prospective users to search a property purpose and it can be recommended to prepare spatial
according to their needs. In the present study a simple databases for the government departments related to the
graphical user interface (GUI) is created wherein land. It can help them to take decisions for their projects
dynamic map are displayed and it allows to perform and can act as a base data for the further functions like
basic GIS operation and user-based queries (Figure 6). Tax calculation, etc. by adding appropriate attributes.
The GUI have been developed using visual basic System integration or a development of centralized
programming language. system for all the departments related to the land records
is important so that benefit could be obtained from
The query window is the main part of the application business as well as technical point of view. REIS can
which has been customized to meet the needs of the also be linked to other databases of different
user. The query window is linked with the database to departments so that user can access the other types of
fetch the requested results and to display it on the map data. This would save the resources, time and effort and
window. The user interface contains basic display can help in the fast decision making. This system can be
components like zoom-in, zoom-out, pan etc. It also improved by adding the market rate, which can help user
provides the facility of spatial query and it provides the a better understanding of the real estate market.
attribute data of any land parcel where the user has
clicked. On click by the user on any part of the map, a
new window pops up and the attribute information is
displayed in that window.
23
http://www.microsoft.com/visualstudio/ Accessed on
22.06.2013
http://msdn.microsoft.com/en-us/library/hh127540.
Accessed on 22.06.2013
Abstract: Imprints of multistructural events recorded within the rocks are visualised through satellite data. The events
of superposed folding and shearing at Koba-Rustampura area belonging to Khandia Formation of Champaner group has
been studied. Classification of such multiple events becomes simple with the help of their respective trends. These
structural events have been delineated by using visual image interpretation techniques to study the spatial pattern and
textures on the Landsat image. By deciphering axial traces and directions of displacement, one can build the chronology
of the structural events revealing the deformational history.
Applications of remote sensing in geosciences are well The true colour composite of Landsat image of 2016
established and widely accepted for identifying and with 30m spatial resolution has been used to identify
correlating structures, regionally. Identification of different structural events present within the study area.
structures on regional scale have been a necessity to Deformational events from later to former have been
understand overall deformational pattern of any terrain. interpreted by studying the spatial pattern on the image
Such attempts includes mapping of large scale features, through visual image interpretation techniques. The
extracting lineament patterns, identifying regional fold information collected during ground truth carried out in
trends, quantifying the fault directions, etc and its 2016, have been used to build the time relationship and
correlation up to plate dynamics (Nama, 2004; Kenea, to understand overall deformational history pertaining
1997; Heddi et al., 1999; Semere and Ghebread, 2006; to the Koba-Rustampura area.
Marghany et al., 2009; Maged and Mazlan, 2010;
Stefouli and Osmaston, 1986; Shuichi, 2002; Stamouslis 3. Regional geological setup
and Rogers, 2003; Yamaguchi and Naito, 2003; Rowan
and Mars, 2003; Gomez et al., 2005; Harding and The Koba-Rustampura area belongs to a part of
Berghoff, 2000; Misra et al., 2014; Joshi et al., 2014). Champaner group, which is well known for its low-
grade meta-sedimentary sequence. The group consists
In order to establish the correlation, in terms of regional of lithological entities such as quartzite, phyllite, meta-
structures, it is prerequisite to appreciate the continuity conglomerate, schist, impure dolomitic limestone and
from meso to micro scale. The present work reports a intermixed variety of granites and gneisses (Gupta et al.,
study on the time relation of structural events in Koba- 1992, 1995; Joshi et al., 2014). Geographically
Rustampura area which is situated 24 Km east of Champaner group is surrounded on three sides (i.e.
Vadodara district, Gujarat (falls under latitude and north, east and south) by younger plutonic intrusive
(Godhra granite) and one side (i.e. west) by Deccan trap
longitude 22o16’38.61” – 22o21’21.18” N and
rocks. Geologically the Champaner group represent an
73o28’25.63” – 73o38’51.38” E respectively). The
example of inlier due to the presence of younger rocks
study area belongs to Khandia formation of Champaner
neighbouring from all sides (Gupta et al., 1997).
group, Aravalli Supergroup having meso-proterozoic
age. The region has experienced polyphase
Structurally, rocks of Champaner group display two
deformational history and are characterised by
significant trends of axial traces. D1 phase of
lithological entities such as phyllite, quartzite, meta- deformation has resulted F1 folds of E-W trend where as
conglomerate. Based on the Landsat image of 2016, D2 phase of deformation has resulted F2 folds of N-S
acquired from the google earth portal, trend. The proximity of F2 folds decreases from eastern
characterisation of rocks holding different structural end to western end of Champaner group (Jambusariya
trends have been attempted. The same have been and Merh, 1967; Gopinath et al., 1977; Merh, 1995;
delineated in a chronological order and supplemented by Shah et al., 1984).
field as well as micro-structural studies.
Figure 3: Field Photograph showing: (a) Closely spaced factures in quartzites. Loc. Koba village; (b) Breccia
containing angular fragments of quartz embedded in ferruginous matrix. Loc. Rustampura village; (c) Elongated
clasts of quatz in meta-conglomerate, ball pen signifies the stretching direction. Loc. Koba village; (d) Group II
mica fish in meta-conglomerate (10XCN); (e) Quartz fish in meta-conglomerate (10XCN); (f) Development of S-C
fabric of mica grains in meta-conglomerate (4XCN); (g) Breccia containing medium to coarse grained angular
quartz clasts in ferruginous matrix (4XCN); (h) ‘V’ pull apart microstructure with domino like arrangement of
quartz clasts, arrows indicated shear direction (4XCN)
27
These fold events have undergone post deformational island, Washington state. Proceedings of the American
shearing along the weak planes. Directions of these Society of Photogrammetry and Remote Sensing,
shears includes top to NW, top to NE and down to SE Annual Conference, Washington, DC, 9 pp.
displacement. Field evidences, such as brittle fracturing
in quartzites and elongation of clasts in meta- Heddi, M., D.J. Eastaff and J Petch (1999).
conglomerate, gives sustainable sense of shear. In Relationships between tectonic and geomorphological
addition to field evidences, supportive microstructural linear features in the Guadix-Basa basin, Southern
evidences are also envisaged. These include dominant Spain. Earth Surface Process and Land-forms, 24,
S-C fabric of oriented mica grains, quartz fish, group 2 pp.931-942.
mica fish and ‘V’ pull apart mechanism with ‘domino’
like microstructure. Jambusaria, B.B. and S.S. Merh (1967). Deformed
greywacke conglomerates of Jaban near Sivrajpur,
Based on the above facts the time relationship of Panchmahals district, Gujarat. Ind. Minerals. 8, pp. 6-
deformational events is established. Chronologically it 10.
can be represented from older to younger as: 1. E-W
trending folds; 2. N-S trending open folds on limbs of Joshi, Aditya and M.A. Limaye (2014). Evidence of
earlier folds; 3. Shearing. syn-deformational granitoid emplacement within
champaner group, Gujarat. Journal of, The Maharaja
Acknowledgement Sayajirao University of Baroda, 49(1), pp. 45-54.
The author is thankful to Prof. L.S. Chamyal, Head, Joshi, Aditya, M.A. Limaye and B.S. Deota (2014).
Dept of Geology, M.S. University, Baroda for providing Structural footprints extraction from rocks of Lunavada
necessary facilities. The author is sincerely thankful to region,Gujarat through IRS LISS III. Journal of,
Mr. Aditya Joshi for academic discussion. Author is also Geomatica, 8(2), pp. 170-173.
grateful to the Forest Department of Vadodara and
Panchmahal District, Gujarat for granting the Kenea, N.H. (1997). Improved geological mapping
permission to carry out fieldwork in the Sanctuary area. using Landsat TM data, Southern Red Sea Hills, Sudan:
Author will never forget the help rendered by Mr. H.J. PC and his decorrelation stretching. International
Bhatt during the fieldwork. Journal of Remote Sensing, 18, pp. 1233-1244.
Gupta, S.N., Y.K. Arora, R.K. Mathur, B.P. Iqbaluddin, Rowan, L.C. and J.C. Mars (2003). Lithologic mapping
T.N. Sahai and S.B. Sharma (1997). The Precambrian in the Mountain Pass, California, area using Advanced
geology of the Aravalli region, Southern Rajasthan and Spaceborne Emission and Reflection Radiometer
NE Gujarat, India. Mem. Geological Survey of India (ASTER) data. Remote Sensing of Environment, 82,
123, pp. 262. pp.350-366.
Harding, D.J. and G.S. Berghoff (2000). Fault scarp Semere, S. and W. Ghebread (2006). Lineament
detection beneath dense vegetation cover: airborne characterization and their tectonic significance using
Lidar mapping of the Seattle fault zone, Bainbridge Landsat TM data and field studies in the central
28
ISG Newsletter
Indian Society of Geomatics (ISG) brings out a newsletter which is very popular because
of its content on geomatics. The newsletter has featured special themes like desertification,
mountain ecosystem, watershed development, climate change etc.
The forth coming issue of ISG Newsletter will feature popular geomatics articles of current
interest.
ISG invites articles of general interest on current topics related to geomatics. The articles
may be sent to:
E-mail: [email protected]
Phone: 02717-235434
29
Abstract: The purpose of this research paper is to identify watersheds with high flood potential based on their
characteristics for formation of surface runoff. The SCS-CN method relies on remote sensing and GIS data for obtaining
watershed characteristics. A 30m raster grid size digital elevation model (DEM) has been generated from field survey
using Global Positioning System (GPS) of 3m accuracy integrating with Survey of India topographical maps of 1:50,000
scale having 10m contour interval. The undisturbed soil samples from field have been collected and laboratory analysis
was carried out using modified proctor compaction test as per ASTM D1557 and sieve analysis as per ASTM C136. This
has helped in establishing hydrological soil map while land use map has been prepared using Landsat 7 ETM+ image
band 2, 3, 4 (30m) merged with PAN band 8 (15m) for classification. The supervised maximum likelihood classification
approach has been employed for preparation of land use map for Varekhadi catchment having 442 km 2 of geographical
coverage. The major land use categories classified on 10 Nov 2001 Landsat 7 ETM+ image were agriculture (32%),
forest (29%), wasteland (20%), fallow land (14%), built-up (4%) and water bodies (2%). The hydrological soil groups
generated in GIS environment have identified two soil groups viz. group B and group C that exist under study area. The
Varekhadi catchment has been delineated into five watersheds viz. Amli, Zankhwaw, Visdaliya, Godsambha and Wareli
delineated using DEM and stream network. The SCS-CN model was applied for estimating of daily run-off for each sub-
watershed. The results obtained on the flood potential analysis shows that Wareli watershed has highest flood potential
while the Amli watershed lowest. It should be noted that highest value of flood potential belongs to lowest part of
watershed, where high population density is found. This analysis reflects an increased vulnerability and risks to floods
and inundations for Wareli watershed. Stream gauge data has been used for result validation with a common event of
2010 and it shows good agreement with the model. The flood potential analysis within the lower Tapi basin tributary
suggests that the SCS-CN method with hydrological parameters derived using remote sensing and GIS data can be applied
to predict run-off in poorly gauged watersheds.
Keywords: SCS-CN method, Remote sensing, GIS, Landsat 7 ETM+, Runoff, Ungauged catchment
P = Q + F + Ia ………………….. (1)
Ia 0.2S
When P S and generally , hence
equation 4 can be written as;
Qp = (P-Ia)2/(P-Ia+S)
(b)
5. Conclusion
Remote sensing and GIS data are of great use for surface
runoff estimation when conventional methods of runoff
estimation are inadequate. Both the techniques have
been used for generating model input for determination
(c) of physical characteristics of watershed such as land
use, hydrology soil group and CN number. Gaussian
maximum likelihood classifier has been used for
classification of land use and shows good field
acceptability. It has been integrated with SCS-CN
method for identification of watershed, estimation of
flood potential for a part of lower Tapi basin.
Acknowledgement
Email: [email protected]
Abstract: Construction of dams brings about changes and this is the situation on the Volta River with the construction
of the Kpong dam. The aim of this study is to detect these changes and the trend of change after some years of construction
of the dam. Topographic map of the lake (1974), Landsat Thematic Mapper (TM) image (1990), Landsat Enhanced
Thematic image (2000) and Landsat (UTM) image (2010) were used. The images were georeferenced and classified into
six (6) classes namely, closed forest vegetation, open forest vegetation, dense shrub, grass/herbaceous cover, built up/bare
surfaces and water body for the purpose of our study, using supervised classification. Maps were produced to show the
changes in the land cover features. Careful observations of the produced maps showed that the most dominant change
happened after the construction of the dam. wherebythe river overflew its boundary submerging some islands and
communities along the river.This caused the Islands to reduce from 340.607ha in 1974 to 319.959ha in 1990 whilst the
lake increased from 2266.398ha in 1974 to 4007.07ha in 1990. After that the water increase has been gradual and all
other changes have been as a result of anthropogenic activities.
Keywords: Volta river, Kpong dam, Supervised classification, Land cover changes, Landsat images
2.2 Methodology
Landsat Thematic Mapper (TM) image of 1990,
Landsat Enhanced Thematic Mapper (ETM) 2000 and
Landsat Thematic Mapper (UTM) 2010 of the study
area on the Volta tiver were obtained. These images
were dry season images captured in the months of
March and April. The images underwent radiometric
and geometric corrections. Figure 2: Map showing the boundaries of the river
at 1974 and 1990 after construction of the dam
The topographic map of 1974 and the 1990 image were
used to produce another map showing the boundaries of
the lake for the two different years. The Landsat images
were classified into six (6) different land covers namely
closed forest, Open forest, dense shrub, grass, built
up/bare lands and water body using maximum
likelihood method of supervised classification.
Classifications for the various images were guided by
observations on the ground, local residents and other
techniques such as PCA and NDVI. Statistics were
generated for these classes on each image to know the
area covered by each land cover. Each classified image
was superimposed with the data of the boundary of the
river.
The accuracy of the 1990 and 2000 images were not Ledec, G. and J.D. Quintero (2003). Good dams and bad
verified but that of 2010, 100 points were picked with dams: Environmental criteria for site selection of
the GPS and used to validate the classification. This hydroelectric projects. Sustainable Development
gave an overall accuracy of 84% and a Kappa of Working Paper16, November 2003
80.25%.
Moshen, A. (1999). Environmental land use change
Acknowledgement detection and assessment using with multi – temporal
satellite imagery. Zanjan University.
Sincere thanks to the Survey department of the Centre
for Remote Sensing and Geographic Information Narumalani, S., Y. Zhou and J.R. Jensen (1997).
System (CERSGIS) and Environmental department of Application of remote sensing and geographic
the Volta River Authority. information systems to the delineation and analysis of
riparian buffer zones. Aquatic Botany 58: 393-409.
References
Prabaharan, S., K. Srinivase Raju, C. Lakshumanan and
Agbenyo, P.K. (2009). Effects of the Akosombo and M. Ramalingam (2010). Remote sensing and GIS
Kpong power schemes on six selected mafi- applications on change detection study in coastal zone
communities in the Volta region. MSc Thesis: The using multi temporal satellite data. International Journal
Institute for Development Studies, Universsity of Cape of Geomatic and Geosciences, Volume 1, No. 2, 2010.
coast, Ghana
Singh, A. (1989). Digital change detection techniques
Amankwaa, O.P.J. (2002). Ghana: A human using remotely sensed data. International Journal of
geography for secondary schools. Takoradi, Ghana: St. Remote Sensing. Vol. 10, No. 6, p. 989-1003.
Francis Press Ltd.
Tortajada, C. (2001). Environmental sustainability of
Campbell, J.B. (1996). Introduction to remote sensing. water projects. Doctoral thesis: Division of Hydraulic
Second edition, The Guilford Press, USA. Engineering, Department of Civil and Environmental
Engineering, Royal Institute of Technology,
Girmay, Y. (2006). Assessing the environmental Stockholm.
impacts of a hydropower project: The case of
40
Abstract: An objective method for detecting night time fog based on bi-spectral difference of 3.9µm and 11µm channel
brightness temperature has been developed. The thresholds used in this are dynamically derived, based on the data and
have been tested with MODIS AQUA/TERRA data for December 2012-January 2013 and December 2013-January 2014
winter season over northern India. The generated fog maps for December 2012-January 2013 have been validated
qualitatively with the fog map generated by IMD (India Meteorological Department) using MODIS data. Quantitative
validation has been carried out for January 2014, against visibility data at five locations of Northern India. Fog was
detected with 70% success rate using this objective method over this region for the validation period.
Keywords: Fog, MODIS, Bi-spectral, Brightness Temperature Difference (BTD), Indo Gangetic Plains (IGP)
3. Methodology Cloud
(a) (b)
Normalized Histogram of 10.8 µm BT Normalized Histogram of 3.9 µm BT
1 1
Normalized Histogram
Normalized Histogram
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0
0
220 240 260 280 300
220 240 260 280 300
BT 10.8 μm BT 3.9 µm
Figure 2: Histogram of brightness temperature in two channels (a) 10.8 µm and (b) 3.9 µm
In order to separate out the fog pixels only, the scatter is also observed that for these values of BT the BTD is
plot between TIR brightness temperature and BTD having large positive value and appear as a separate
values of the image is studied. Figure 3 shows the scatter cluster. It is a remarkable feature (Marked as ellipse in
plot between TIR BT and BTD for a typical case having figure 3) which separates out the foggy pixels from the
fog, clear pixels, snow and cloud. From the scatter plot underlying land pixels. The information in this cluster is
it is seen that for brightness temperature less than 270° used to find out the minimum threshold value for BTD
K, the brightness temperature difference values are and the pixels having BTD value more than the
mostly negative or takes a very small positive value. minimum threshold value is classified as fog. The
This corresponds to the high level clouds and snow. ΔTBmin is determined automatically from the histogram,
of the BTD image. In a typical fog image, the histogram
of BTD image shows one primary maximum, a
secondary maximum and a minimum in between (figure
4a). In such a condition, BTD corresponding to the
BTD (BT 10.8 -BT 3.9 ) 0K
However, the MODIS data which has been used for this
analysis is in tile form and for each acquisition the frame
of the image changes and thus the proportion of
different features in the image also changes. Therefore,
for some cases in the histogram of BTD image of
BT 10.8 0K MODIS data, only single maximum is obtained (figure
Figure 3: Scatter plot of brightness temperature of 4b). In such cases the threshold for BTD (ΔTBmin) is set
10.8µm and brightness temperature difference (The to a point where the slope of the clear sky pixel drops
cluster marked within the ellipse indicates the considerably and ΔTB more than ΔTBminis classified as
presence of pixels with fog) fog. By applying this methodology fog region has been
identified with MODIS nighttime data from 01
However, for land, low clouds and fog regions where December 2012 to 31 January 2013 as well as for 01
the brightness temperature is more than 280°K, the December 2013 to 31 January 2014 and compared with
brightness temperature difference is positive. Due to MODIS fog image of IMD (India Meteorological
temperature inversion the brightness temperature of fog Department, www.imd.gov.in), whenever available as
is either similar or greater than that of the underlying well as with the early morning MODIS RGB data of
land. Thus the pixels with TIR brightness temperature next day. However, a quantitative validation of fog map
greater than 280 and less than 290° K, will correspond generated for January 2014 has been carried out with
to land pixels over laid by fog or low cloud pixels during visibility data at five different locations i.e. Delhi,
night (as observed from the figure 3). From the figure it Lucknow, Jaipur, Amritsar and Varanasi situated in
northern India.
43
(a) (b)
Histogram of Difference Image Histogram of Difference Image
1 1
Normalized Histogram
Normalized Histogram
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 0
-15 -10 -5 0 5 10 15 -10 -5 0 5 10
Figure 4: Histogram of brightness temperature difference image (a) BTD image showing double peak (b) BTD
image showing single peak
4. Results and discussion shallow, dense and very dense depending on the change
in brightness temperature difference. The minimum
The evolution of fog over Indo-Gangetic plain in difference in brightness temperature ΔTBmin which is
December 2012 and January 2013, using MODIS obtained to detect fog varies from 3° to 8° K for the
Terra/AQUA data (Different colour shows the BTD period of analysis. For each image the threshold is found
values) is shown in figure 5. In order to show the to be different. This change in threshold is attributed to
evolution of fog to the maximum extend over the IG the change in fog optical depth and droplet size
plains, both the fog map generated using MODIS (Chaurasia et al., 2011). BTD value of 3° K corresponds
AQUA/TERRA has been shown in this figure. Care has to very thin fog. However, for BTD values between 4 to
been taken to incorporate maximum fog period to show 5° K corresponds to moderately dense fog. BTD value
the evolution. Therefore, few days with in December of 6 and >6 °K represents very thick fog. This gradation
2012 and January 2013 is missing in the figure 5. As is empirically made based on our earlier study
observed from the fog map generated by the above (Chaurasia et al., 2011) as well as from the visibility
method, the 2012-2013 fog periods started around early data.
December. On December 09, 2012, fog was initiated in
the eastern states like Bihar and gradually it extended The generated fog map has also been compared with the
spatially and covered Bihar and part of Uttar Pradesh by next day morning RGB image and fog map generated
December 15, 2012. The spatial extent continued to by IMD using MODIS data. The qualitative inter-
increase both from east to west as well as in the southern comparison of the two maps is shown in figures 6 and 7
direction starting from the foot hill of Himalaya. During indicating good agreement with each other. Bi-spectral
December 22-25, 2012 it covered the entire Indo- differencing technique has also been used in the
Gangetic plain starting from Bihar at the east to Punjab generation of fog map by IMD using MODIS data, and
at the west. This is the time when the wheat cultivation a pixel is classified as fog/low stratus when the BTD is
in the four major wheat growing states Bihar, UP, greater than 2.50 K (Product catalog, www.imd.gov.in ).
Haryana and Punjab are in full swing and the fields are This threshold is fixed for all images. In order to make
irrigated, which adds to the local moisture availability the two maps comparable instead of showing the
apart from the prevailing meteorological condition and variation of BTD over fog, only fog and no- fog regions
temperature inversion in this region which is conducive are generated. The blue colour in both the images shows
for fog formation. Similar situation continued till fog. It is observed from the figure 7 that, for December
January 07, 2013. 25, 2012, there is a systematic data loss in IMD fog map.
It is also observed that it has also picked up some of the
On January 10, 2013 part of north-western region of high cloud part in the southern part of the image, which
Uttar Pradesh was covered with fog which extended in is filtered by new developed dynamical thresholding
width till January 14, 2013. On January 19, 2013 there algorithm. Similarly, the cloud captured by IMD fog
was heavy fog condition over Punjab and parts of map on January 10, 2013 is also not captured by our
Haryana. On January 21, 2013 again Uttar Pradesh algorithm.
region was engulfed with fog which continued till
February 2013 (results shown till end of January). The A quantitative validation of the generated fog map has
gradient in the fog map shows the category of fog i.e. been carried against the in-situ visibility data over five
44
15 Dec 2012, 1700UTC 19 Dec 2012, 1640UTC 21 Dec 2012, 1625UTC 22 Dec 2012, 1710UTC
23 Dec 2012, 1750UTC 24 Dec 2012, 1655UTC 25 Dec 2012, 2020UTC 29 Dec 2012, 1715UTC
31 Dec 2012, 1715UTC 01 Jan 2013, 1745UTC 02 Jan 2013, 1650UTC 14 Jan 2013, 1715UTC
23 Jan 2013, 1710UTC 24 Jan 2013, 1615UTC 26 Jan 2013, 1600UTC 27 Jan 2013, 1645UTC
Figure 5: Evolution of fog over Indo-Gangetic plain in December 2012 and January 2013, using MODIS Terra/AQUA data
(Different colour shows the BTD values)
45
(a) (b)
Figure 6: Inter-comparison of (a) generated fog map using MODIS Terra/AQUA of night time (Different colour shows the
BTDvalue) with (b) next day morning MODIS RGB image (recoloured)
(b)
13 Dec 2012, 1955UTC 15 Dec 2012, 1657UTC 20 Dec 2012, 1720UTC 25 Dec 2012, 2015UTC 10 Jan 2013, 2015UTC
Figure 7: Inter-comparison of generated fog map using (a) MODIS AQUA data with that of (b) IMD fog map using
MODIS data (recoloured)
46
Total no. of Observations 140 Ellord, G.P. (1995). Advances in the detection and
analysis of fog at night using GOES multi spectral
No. of non-concurrent data 25 infrared imagery. Weather Forecasting,10, 606-619.
Number of Hits 81
Eyre, J.R., J.L. Brownscombe and R.J. Allam (1984).
Number of Miss 16 Detection of fog at night using Advanced Very High
Number of false detection 18 Resolution Radiometer. Meteorological Magazine, 113,
266–271.
Percent of Detection (POD) 70%
Percent of False Detection 15% Greenwarld, T.J. and S.A. Christopher (2000). The
GOES I-M imagers: New tools for studying
Percent of Miss 12.8% microphysical properties of boundary layer stratiform
clouds. Bulletin of American Meteorological Society,
5. Conclusions 81(11), 2607-2619.
This technique yields a very good probability of night Gultepe, I., S.G. Cober, G. Pearson, J. A. Milbrandt, B.
time fog detection with acceptable false alarm Hansen, G.A. Isaac, S. Platnick, P. Taylor, M.Gordon
conditions. It has been observed that with the proposed and J.P., Oakley (2009). The Fog Remote Sensing and
dynamical thresholding, for 70% of the cases the Modelling (FRAM) field project and preliminary report.
generated fog map successfully detects fog using Bulletins of American Meteorological Society, 90(3),
satellite IR imagery during night time. In future, it may 341-359.
be required to develop algorithm to minimize mis-
classification between low clouds and fog. Hunt, G.E. (1973). Radiative properties of terrestrial
clouds at visible and infra-red thermal window
Acknowledgements wavelengths. Quarterly Journal of Royal
Meteorological Society, 99, 346-369.
The authors wish to thank the Director, Space
Applications Centre, for his keen interest in this work. Lee, T.F., F.J. Turk and K. Richardson (1997). Stratus
We are also thankful to www.ladsweb.nascom.nasa.gov and fog products using GOES-8-9 3.9µm data. Weather
for MODIS data and India Meteorological Department Forecasting, 2, 664-677.
for providing the visibility data.
Levy, W.F., L.A. Remer, S. Mattoo, E. Vermote and
References Y.J. Kaufman (2007). Second-generation operational
algorithm: Retrieval of aerosol properties over land
Bendix, J. (2002). A satellite-based climatology of fog from inversion of moderate resolution imaging
and low-level stratus in Germany and adjacent areas. specroradiometer spectral reflectance. Journal of
Atmospheric Research, 64, 3-18. Geophysical research, 112: D13211,
DOI:10.1029/2006JD007811.
Bendix, J. and M. Bachmann (1991). A method for
detection of fog using AVHRR-imagery of NOAA Ram Kirpa, M., M. Sarin, A. K.Sudheer and R.
satellites suitable for operational purposes. Rangaranjan (2012). Carbonaceous and secondary
MeteorologischeRundschau, 43, 178 (in German), inorganic aerosols during wintertime fog and haze over
1991. urban cities in Indo-Gangetic plain. Aerosol and Air
Quality Research, 12, 359-370.
Cermak, J. and J. Bendix (2007). Dynamical nighttime
fog/low stratus detection based on Meteosat SEVIRI Schueler, C.F. and W.F. Barnes (1998). Next generation
data – a feasibility study. Pure and Applied Geophysics, MODIS for polar operational environmental satellites.
64, 1179-1192. Journal of Atmospheric and Ocean Technology, 5, 430-
439.
Cermak, J. and J. Bendix (2008). A novel approach to
fog/low stratus detection using Meteosat 8 data. Singh, B., H.N. Dutta and A. Kaushik (2004).
Atmospheric research, 87, 279-292. Ecological restoration by fog water. Presented in XIII
National Space Science Symposium (NSSS) Koteyam.
Chaurasia, S., V. Sathiyamoorthy, B. Paul Shukla, B.
Simon, P.C. Joshi and P.K. Pal (2011). Night time fog Tiwari, S., S. Payra M. Mohan, S. Verma and D.S.
detection using MODIS data over Northern India. Bhisht (2011). Visibility degradation during foggy
Meteorological Applications, 8(4), 483–494.
47
Abstract: The Varuna is one of the important interfluves rivers joining the Ganga river near Varanasi city. The
morphometric analysis of the Varuna river basin has been carried out. It covers an area of about 3622 km2 of the Ganga
plain. The dendritic drainage pattern diagnosed in the area exhibits homogeneous permeable substratum and gentle
gradient. Spatio temporal changes (Land Use and Land Cover i.e. LULC) of the Varuna river basin, in Varanasi district,
using Landsat multispectral imageries spanning 38 years (1972, 1988, 2002 and 2010) are also studied. The LULC
patterns illustrate that in early 80’s the basin is largely covered by the salt affected waste land. Later on the salt affected
waste land area is reduced and area of agricultural land and built-up area has increased. The study underlines the necessity
of a scientific approach for the sustainable river basin management, with the help of the hydrological conditions, recent
climatic anomalies and geological control of the basin.
Figure 1: Location map of the Varuna basin (a) Map showing subdivision of Ganga plain in Uttar Pradesh
(modified after Singh, 1996); (b) Landsat imagery (MSS Path 153 and Row 42 of 1972) representing Varuna river
basin
Figure 2: Varuna river basin (a) DEM (digital elevation model) from SRTM data; (b) Drainage pattern and 4 sub
basins
Figure 3: (A) Two discreet southward shifting of the confluence of the Varuna river near Varanasi, traced out in
the field with hand GPS; (B) Traces (in doted lines) of the old Varuna River near Sarnath (after Jayaswal et al.);
(C) Present confluence of the Varuna river; (D) Old confluence of the Varuna River; (E) Abandoned channel of
the Varuna River; and (F)Low lying area showing abandoned channel of the Varuna river
Dendritic drainage pattern is the most common and The variations among the stream length ratio (Rl) values
widespread pattern found in the study area (Fig.2b). The are directly related to the geomorphology including
dendritic patterns evolved in the area closely resemble topography and lithology and hence it governs the
to the area having homogeneous bed materials (mainly erosional stage of the watershed and discharge
Gangetic alluvial) with a very gentle regional slope. The (Sreedevi et al., 2004). The Rl values in case of Varuna
perimeter of the Varuna watershed is 482.07 km (Table river basin varies between 0.66 and 0.75 (Table 1) and
1). The values of the bifurcation ratio (Rb) generally set that implies accomplishment of geomorphic maturity.
in between 2.0 and 5.0 for the drainage network. Such The Rho coefficient is used for determining the storage
network generally develops in consistent lithologies and capacity of drainage network (Horton, 1945). The
also signifies the minimum structural control over it. values of Rho coefficient of the Varuna river and its sub-
When the values are higher than 10.0, it indicates that watersheds vary from 0.14 to 0.26 (Table 1). The higher
structural control has played a dominant role on values of Rho coefficient of SW2 and SW3 (Fig.2b) are
drainage network development (Mekel, 1970). Low Rb indicators of higher hydrologic storage during floods
values indicate elongated shape of the basin (Morisawa, and decreased effects of erosion during elevated
1985). The Varuna river belongs to this category. discharge.
51
Table 1: Mathematical formula adopted for the quantitative measurement and estimated values of the
morphometric parameter
The GIS techniques have been efficiently utilized for Khan, A.A., P.C. Nawani and M.C. Srivastava (1988).
the assessment of the drainage characteristics of the Geomorphological evalution of the area around
Varuna river watershed and to comprehend the Varanasi, Uttar Pradesh, with aid of aerial photography
significance of morphometric studies and shifting of the and Landsat imageries. Geological Survey of India,
confluence. The low bifurcation ratio and elongation Rec. 113(8), pp. 31-39.
ratio of the Varuna river signify the elongated shape of
the basin. The elongated shapes of sub- watersheds, Kondolf, G.M. and P. Hervè (2003). Tools in fluvial
with the larger basin area, are due to head-ward erosion geomorphology. John Wiley & Sons Ltd., England.
leading to expansion of channels. The study area is
demonstrating very low drainage density because of Kotoky, P., D. Bezbaruah, J. Baruah and J.N. Sarma
permeable subsurface and dense vegetation. The low (2005). Nature of bank erosion along the Brahmaputra
drainage density in the study area leads to very coarse river channel, Assam, India. Current Science 88 (4), pp.
drainage texture. The low values of circulatory ratio are 634-640.
characterized by attenuated flood- discharge periods
and youth stages of watershed development. Thus the Mekel, J.F.M. (1970). The use of aerial photographs in
incorporation of morphometric analysis together with geological mapping. ITC text book of photo-
predictable watershed assessment methods would have interpretations 8, pp. 1-169.
a better understanding the status of land form, drainage
54
Schumm, S.A., J.F. Dumont and J.M. Holbrook (2000). Swamee, P.K., B. Parkash, J.V. Thomas and S. Singh
Active tectonics and alluvial rivers. Cambridge (2003). Changes in channel pattern of river Ganga
University Press, Cambridge. between Mustafabad and Rajmahal, Gangetic plains
since 18th century. International Journal of Sediment
Shukla, U.K., I.B. Singh, M. Sharma and S. Sharma Research 18(3), pp. 219-231.
(2001). A model of alluvial megafan sedimentation:
Ganga Megafan. Sedimentary Geology 144, pp. 243– Thomas, J., S. Joseph and K.P. Thrivikramaji (2010).
262. Morphometric aspects of a small tropical mountain river
system, the southern Western Ghats, India. International
Shukla, U.K. and N.J. Raju (2008). Migration of the Journal of Digital Earth 3(2), pp. 135-156.
Ganga river and its implications on hydrogeological
potential of Varanasi area, U.P. India. Journal of Earth Turner, M.G., W.H. Romme, R.H. Gardner, R.V.O’
System Science 117, pp. 489–498. Neill and T.K. Kratz (1993). A revised concept of
landscape equilibrium: Disturbance and stability on
Shukla, U.K., P. Srivastava and I.B. Singh (2012). scaled landscapes. Land Ecology 8(3), pp. 213- 227.
Migration of the Ganga river and development of cliffs
in the Varanasi region, India during the late Quaternary: Wells, N.A. and J.A. Don (1987). A reconnaissance of
Role of active tectonics. Geomorphology 171–172, pp. sedimentation on Kosi alluvial fan of India. Special
101–113. publication of SEPM, Tulsa 39, pp. 51-61.
Abstract: Malaria is a major global health problem. About 3.5 billion people are at risk of infection of malaria worldwide,
with environmental factors contributing for about 70–90% of the disease-risk. Over one million cases of malaria are
reported each year, out of which more than 80% is from the sub-Saharan Africa. Ethiopia is a predominantly malaria-
prone country as about 75% of the landscape of the country is favorable for breeding of the malaria vector. Geographical
Information System has emerged as a spatial technology, which integrates a wide range of datasets available from
different sources including remote sensing and Global Positioning System. During the present study, a malaria-risk map
of Mecha district of Ethiopia was prepared by establishing the relationship of various climatic and non-climatic factors
related to the disease using regression analysis. Weighted overlay technique of multi-criteria evaluation was used to
develop the malaria-risk map. Temperature, rainfall, altitude, distance from streams, distance from swamps and ponds,
population density, health facilities and land-use/land-cover patterns were used to prepare malaria-risk areas. Malaria
hazard, elements of risk and vulnerability layer were overlaid, and further verified by ground truth and village-wise
reports of malaria cases to produce the final malaria-risk map. Four categories of malaria-risk ranging from very high to
low were derived. Most of the study area (99.01%) was found to belong to high and moderate malaria-risk. It is suggested
that effective identification and mapping of malaria-risk levels can be made using geospatial tools, to contribute for the
prevention and control of this disease.
Keywords: GIS, Landsat, Malaria, Remote sensing, Regression analysis, Weighted overlay
Vulnerability map
Malaria-hazard map Element at risk map
AHP
Weighted overlay
Malaria-Risk map
Image processing starting from image pre-processing environmental conditions for malaria transmission
(geometric and radiometric correction), layer stacking based on climatic and non-climatic factors. All factor
(band 2–7), image enhancement and image parameters compatible to hazard analysis were
classification to the final accuracy assessment were generated before weighted overlay. Hazard map was
done in ERDAS Imagine software. To correct the produced by computing and reclassifying the five
satellite image radiometrically, Operational Land parameters viz. meteorological (rainfall and
Imager (OLI) band data were converted to top of temperature) data, distance from ponds and swamps,
atmosphere (TOA) to a planetary reflectance using altitude and distance from streams layers. Each of the
reflectance rescaling coefficients provided in the hazard parameters was ranked according to the
product metadata file. The following equation was used importance for mosquito breeding and transmission.
to convert DN values to TOA reflectance: The process of weighting each factor was performed in
IDRISI software. After assigning weight, the hazard
𝜌𝜆′ = 𝑀𝜌𝑄𝑐𝑎𝑙 + 𝐴𝜌 (3) map was computed by overlaying the five selected
factors.
The ρλ' does not contain the sun angle correction and
hence the image was again converted to TOA Vulnerability map was generated from distance from
reflectance with a correction for the sun angle as health facility map and population density map. The two
follows: layers were overlaid with 54% weight to population
density map and 46% to health facility map. The weight
𝜌𝜆′ 𝜌𝜆′ was given by consulting health experts, who have
𝜌𝜆 = = (4)
𝑐𝑜𝑠(𝜃𝑠𝑧) 𝑠𝑖𝑛(𝜃Ѕ𝐸) advanced knowledge about malaria based on the
regression result coefficients and available information.
After preparing the factors and validating the model and The element at the risk-map was computed by
image processing, three malaria-risk layers were reclassifying the land-use/land-cover pattern of the
generated using their factors. The malaria-risk map was study area. The land-use/land-cover types were ranked
developed by combining the suitability of
58
R=V×H×E (5)
3. Results
Table 1: Malaria cases and rainfall and temperature in Mecha district (2002–2012)
Months July Aug. Sep. Oct. Nov. Dec. Jan. Feb. Mar. Apr. May June
Malaria case 39164 2770 3638 4760 5705 4487 3313 3162 2909 3070 7188 7683
Rainfall 418.9 388.3 238.2 100.3 35.4 10.4 2.9 2.06 29.7 60.9 175.6 358.5
Temperature 20.07 20.5 20.78 19.72 15.52 17.45 18.94 19.89 22.5 23.78 22.76 20.69
3.3 Regression analysis for model validation and adjusted R2 for this model were 657.72329, 0.8052
As shown in Table 2, rainfall, temperature and and 0.7566, respectively. Multiple R-squared and
population density have strong positive relationship and adjusted R-squared were both statistics derived from the
altitude and slope have strong negative relationship with regression equation to quantify model performance. In
malaria incidence. These are the main factors for this model, R2 was 80.5299%. Hence, 80.53% variation
malaria prevalence in the study area relative to other in the dependent variable (malaria cases) could be
factors analysed (Tables 2 and 3). The AIC, multiple R2 explained by the model.
59
Bautista, C.T., A.S. Chan, J.R. Ryan, C. Calampa, M.H. Martens, W.J., L.W. Niessen, J. Rotmans, T.H.Jetten
Ropers, A.W. Hightower and A.J. Magill (2006). and A.J. McMichael (1995). Potential impact of global
Epidemiology and spatial analysis of malaria in the climate change on malaria risk. Environmental Health
Northern Peruvian Amazon. American Journal of Perspective 103:458−64.
Tropical Medicine and Hygiene 75:1216−22.
Meron, M. (2010). Web GIS in decision support to
Carlos, A.V., A.F. Julián, M.C. Zulma, R. Patricia, C.L. control malaria: Case study in TiroAfetaWoreda,
Myriam and D. Sofia (2010). Correlation between Oromia Region, Ethiopia, MSc. Thesis, Addis Ababa
malaria incidence and prevalence of soil-transmitted University, Addis Ababa.
helminthes in Colombia: An ecologic evaluation.
Biomédica 30:501–508. Oliveira, E.C., E.S. Santos, P. Zeilhofer, R. Souza-
Santos and M. Atanaka-Santos (2013). Geographic
Ceccato, P., S. Connor and M.C. Thomson (2005). information systems and logistic regression for high-
Application of geographical information systems and resolution malaria risk mapping in a rural settlement of
remote sensing technologies for assessing and the southern Brazilian Amazon. Malaria Journal,
monitoring malaria risk. Parassitologia 47: 81−96. doi:10.1186/1475-2875-12-420.
Dongus, S., D. Nyika, K. Kannady, D. Mtasiwa, H. Omukunda, E., A. Githeko, M.F. Ndong, E.
Mshinda and U. Fillinger (2007). Participatory mapping Mushinzimana and G. Yanr (2013). Effect of swamp
of target areas to enable operational larval source cultivation on distribution of anopheline larval habitats
management to suppress malaria vector mosquitoes in in Western Kenya. Journal of Vector Borne Disease 49:
Dar es Salaam, Tanzania. International Journal of 61–71.
Health Geographics 6:37. doi:10.1186/1476-072X-6-
37. Palaniyandi, M. (2012). The role of remote sensing and
GIS for spatial prediction of vector-borne diseases
Ehlkes, L., C. Anne, B. Kreuels, R. Krumkamp and O. transmission. Journal of Vector Borne Disease 49: 197–
Adjei (2014). Geographically weighted regression of 204.
land cover determinants of Plasmodium falciparum
transmission in the Ashanti Region of Ghana, Patz, J.A., A.K. Githeko, J.P. McCarty, S. Hussein, U.
International Journal of Health Geographics 13:13−35. Confalonieri and N. Wet (2003). Climate change and
infectious diseases. Climate Change and Human Health:
Erin, M., S. Jennifer, G. Katya, Y. Amrish, B. Teun, O. Risks and Responses103–137.
Wycliffe, K. Simon, D. Chris, A. Thomas, J.C. Smith
and C. Nakul (2014). The cost effectiveness of malaria Richard, Y. and I. Poccard (1998). A statistical study of
control interventions in the highlands of Western NDVI sensitivity to seasonal and inter- annual rainfall
Kenya. PLoS One.; 9: e107700. variations in Southern Africa. International Journal of
Remote Sensing 19:2907–2920.
Hay, S.I., J.A. Omumbo, M.H. Craig and R.W. Snow
(2000). Earth observation, geographic information Saxena, R., B.N. Nagpal, A. Srivastava, S.K. Gupta and
system and Plasmodium falciparum malaria in sub- A.P. Dash (2009). Application of spatial technology in
Saharan Africa. Advances in Parasitology 47:173–215. malaria research & control: Some new insights. Indian
Journal of Medical Research 130:125−132.
Kathleen, W. (2002). A review of control methods for
African malaria vectors. Environmental Health Project, Sithiprasasna, R., W.J. Lee, D.M. Ugsang and K.J.
Activity Report 108, U.S. Agency for International Linthicum (2005). Identification and characterization of
larval and adult anopheline mosquito habitats in the
Development, Washington.
Republic of Korea: Potential use of remotely sensed
data to estimate mosquito distributions. International
Kleinschmidt, I., B.L. Sharp, G.P.Y. Clarke, B. Curtis
Journal of Health Geographics 4:1−11.
and C. Fraser (2001). Use of generalised linear mixed
models in the spatial analysis of small area malaria
Solomon, N., K.V. Suryabhagavan and E. Gadisa
incidence rates in KwaZulu Natal, South Africa.
(2015). Visceral Leishmaniasis (Kala-Azar) risk
American Journal of Epidemiology 153:1213−21.
mapping using geo-spatial tools: A case study in Kafta
Humera district, North Western Ethiopia, International
Kleinschmidt, I., J. Omumbo, O. Briët, N. Van de
Journal of Geoinformatics 11:21−26.
Giesen, N. Sogoba and N.K. Mensah (2001). An
empirical malaria distribution map for West Africa.
Stephen, M. (2006). Landscape determinants and
Trop Med Int Health 6:779–86.
remote sensing of mosquito larval habitats in the high
64
Abstract: Cultivation of citrus orchards is characterized by small field size, sparse distribution with different age groups
and coexisting spectrally similar crops. Thus, classification of citrus crop using satellite data is quite challenging. The
present study was carried out to evaluate object based classification techniques for mapping citrus orchards located at
Indi hobli of Bijapur district, Karnataka state. Two date LISS-IV and single date Cartosat-1 data were used for
classification. Spectral signature of young orchards less than 5 to 6 years and poorly managed orchards were mixed with
pomegranate, sugarcane and grape orchards. Two approaches of segmentation techniques namely, threshold & clump
and lambda schedule, were tested. The results revealed that the single date satellite data showed classification accuracy
(around 75%) using both threshold and clump and lambda schedule segmentation approach. Inclusion of second date data
along with vegetation indices significantly improved the mapping accuracy (around 85%) by eliminating short duration
crops from evergreen citrus orchards. The study explored the potential use of high resolution data for inventory of citrus
orchards and the methodology could be refined for operational application using textural and contextual information.
Keywords: Object based image analysis, Full lambda schedule Segmentation, Threshold and clumping, Citrus, Brovey
transform
where,
Oi is the area of region i.
Oj is the area of region j.
Figure: 1a) LISS IV May Figure: 1b) LISS IV ui is the average value in region i.
November uj is the average value in region j.
║ui-uj║is the Euclidean distance between the
Classification was done classification for the data spectral values of regions i and j.
obtained in the month of May (figure 1a) and also using length(∂(Oi,Oj)) is the length of the common boundary
the stacked data of May and November. Summer crops of Oi and Oj (Robinson, 2002).
get harvested in the month of May and monsoon crops
get harvested in the month of November (figure 1b). While not necessarily derived from the Pixel Probability
Therefore, these two months data was used to exclude Layer, the raster object segments will have the zonal
seasonal crops as much as possible in the images. mean pixel probabilities as attributes. Output from the
Probability Pixels to Raster Objects Operator contains
LISS IV sensors have three bands with a spatial pixels that are grouped as raster objects which have
resolution of 5.8 m. CARTOSAT-1 panchromatic has associated probability attributes (Erdas, 2013).
one broad band with very good spatial resolution—2.5
m. Combining these two images to yield a three band 3.2 Threshold and clumping
data set with 2.5 m resolution provides the best This operator performs a threshold operation on a pixel
characteristics of both sensors (NRSC, 2004). probability layer which keeps only pixels with a
probability greater than or equal to the threshold value
Brovey Transform was used to merge the two data sets. specified. It converts the pixels to binary values (0’s and
The Brovey Transform was developed to visually 1’s), then performs a contiguity operation (clump) on
increase contrast in the low and high ends of an image’s the binary values of 1(Imagine Objective, 2010).
histogram (i.e., to provide contrast in shadows, water
and high reflectance areas such as urban features). Now the objects formed were edited from the above two
Consequently, the Brovey Transform should not be used methods by applying filters over the objects to refine the
if preserving the original scene radiometry is important. results. Focal filter was used in the beginning to make
However, it is good for producing RGB images with a the object edges sharp and remove very narrow strips
higher degree of contrast at the tails of the image from the objects formed. This filter replaces the most
histogram. repeating value in the window of specific size around
the pixel. The optimum value for this was found to be 3.
Our objective was to perform a binary class The next filter used was probability filter which
classification for the identification of Citrus orchards in removes the objects of lower probability and gives us
the given area. Object based analysis of the image was more reliable objects after filtering. Later size filter was
carried out to identify the citrus orchards and mark all used to filter out very small objects. Finally, re-clump
the other regions as the background. filter was used on objects to form bigger objects. Raster
to vector conversion was made using polygon tracing
3. Methodology and work flow and later two more filter are applied on the vector
objects. The smoothening filter which eliminates the
In ERDAS imagine, Objective Workstation was used. sharp edges of objects and the island filter takes out
The first step was pixel level binary classification of the island like gaps in the formed objects were used in post
image and assigns probability for each and every pixel processing operation. The work flow is pictorially
depending on how close it is to Citrus class. For this step represented in a flowchart in the Appendix-1.
sixtyfour training samples were given from both Citrus
and Non Citrus equally covering all the variations. The 4. Results and discussions
next step i.e., objects formation step which is an
important step in the analysis. Two different methods The main aim of the entire process was to classify the
were followed to create objects in this step. One was mature citrus orchards. It can be seen how different
segmentation using full lambda schedule and the second parameter sets (Table 1) had performed in citrus
was threshold and Clump. dominant and non-citrus dominant regions from images
shown (Figure 2a, 2b).
67
Table 1: Object creation methods, parameters used and respective area statistics and accuracy obtained in
identifying citrus orchards
Single Date/ Object Creation No. Spectral Texture Size Shape Probability Estimated Mapping
Multi Date Method Threshold area (ha) Accuracy
(%)
Set1 0.8 1 0.1 0.05 NA 2490 61.81
Segmentation- Set2 0.6 1 0.15 0.05 NA 1146 63.59
Lambda Schedule
Single Date Set3 0.6 1 0.15 0.05 NA 1140 63.27
Threshold/ Clump Set4 NA NA NA NA 0.95 1888 95.24
Set5 NA NA NA NA 0.9 1152 63.91
Set6 NA NA NA NA 0.95 1557 86.40
Segmentation- Set7 0.5 1 0.02 0.02 NA 843 46.79
Lambda Schedule Set8 0.6 1 0.15 0.05 NA 1402 77.79
Multi Date
Set9 0.8 1 0.15 0.05 NA 1411 78.30
Threshold/ Clump Set10 NA NA NA NA 0.9 1460 81.03
Multi Date Segmentation- Set11 0.7 1 0.15 0.05 NA 1440 79.90
+ NDVI Lambda Schedule
4.1 Lambda schedule segmentation (single date) increasing the probability filter value to 0.6 and
In raster object operation (ROO) process the probability reducing the spectral value to 0.6, which gave very less
filter values were tested iteratively from 0.45 to 0.8. The area estimate 1146 hectares as compared with 1800
optimized value for probability filter was found to be hectares of reported area. However, misclassification of
0.6 for lambda schedule and 0.7 for the Threshold and non-citrus regions was reduced.
Clump process. For Lambda Schedule it was observed
that by giving 0.7 as probability filter value area under 4.2 Threshold and clump segmentation (single date)
citrus drastically changed resulting in under In threshold and clump the threshold probability in ROC
classification and decreasing it to 0.5, resulted in was initially set to be 0.95 (table 2) and probability
misclassifying non-citrus regions. Hence the optimum value of 0.6 in ROO, which resulted in smaller object
value for this parameter was found to be 0.6. In size and misclassification of non-citrus segments. Even
threshold and clump method the value was increased to though the area estimated looked very precise (95.24 per
0.7 because in this step, only high probable pixels were cent) spatially there was significant misclassification.
taken into account in the raster object creation (ROC) Hence, the probability filter value in ROO step was
step itself which precedes ROO step in segmentation. increased to 0.8 which reduced the misclassification. It
Hence, higher values were needed to be kept as cut-off also had smaller segments of objects existing adjacently
for probability filter. which could join to form a bigger object if threshold in
the object formation step is reduced. To achieve this,
The set 1 estimated the area as 2490 hectares and overall there was a need to decrease probability value in ROC.
accuracy was 73.13 per cent. But it was observed that a Hence, two changes were made i.e., the threshold
lot of non-citrus regions were classified as citrus probability was set to 0.9 and probability filter in ROO
orchards. Hence, changes were made in probability was set to 0.7.
filter in ROO step and spectral value in ROC step by
68
` Figure 2a Figure 2b
Figure 2: From top left to bottom right in each of the figures show the classification results for 11 parameter
sets from (a) citrus dominant region; and (b) non-citrus dominant region
Table 2: Table of confusion matrices for selected sets using point based assessment method
Lambda Schedule (Set 1) Threshold and Clump(Set 5)
Citrus Non citrus Total Citrus Non citrus Total
Classified classified
Citrus 31 15 46 Citrus 27 11 38
Non citrus 3 18 21 Non Citrus 7 22 29
Total samples 34 33 Total samples 34 33
taken taken
Accuracy 91.17 54.55% 73.13% Accuracy 79.41 66.67% 73.13%
% %
Lambda Schedule (Set 2) Lambda Schedule (Set 9)
Citrus Non citrus Total Citrus Non citrus Total
classified classified
Citrus 28 10 38 Citrus 29 9 38
Non citrus 6 23 29 Non Citrus 5 24 29
Total samples 34 33 Total samples 34 33
taken taken
Accuracy 82.35 69.70% 76.12% Accuracy 85.29 72.72% 79.10%
% %
Threshold and Clump (Set 10) Lambda Schedule (Set 11)
Citrus Non citrus Total Citrus Non citrus Total
classified classified
Citrus 30 12 42 Citrus 31 6 37
Non citrus 4 21 25 Non Citrus 4 26 30
Total samples 34 33 Total samples 34 33
taken taken
Accuracy 88.23 63.63% 76.12% Accuracy 91.14 78.90% 85.07%
% %
69
APPENDIX – 1
1
Department of Civil Engineering, Gujarat Technological University, Ahmedabad - 382424, India
2
Department of Geotechnical Engineering, University of Transport Technology, Hanoi, Vietnam
3
Geographic Information System Group, Department of Business Administration and Computer Science, Telemark
University College, HallvardEikasPlass 1, N-3800 BøiTelemark, Norway
4
Bhaskarcharya Institute for Space Applications and Geo-Informatics (BISAG), Gandhinagar, India
Email: [email protected] or [email protected]
Abstract: The main objective of this study is to apply and evaluate the predictive capability of the Support Vector
Machines (SVM) and Naïve Bayes Trees (NBT) methods for spatial prediction of landslides in a part of Uttarakhand
state (India). SVM is one of the most efficient machine learning methods that has been applied widely in landslide
prediction whereas NBT has not been applied for landslide problems. In these models, a total of 430 historical landslide
locations have been first identified to construct landslide inventory map. Landslide locations have been split randomly
into two parts to generate training dataset (70% landslide locations) and testing dataset (30% landslide locations).
Secondly, landslide affecting factors such as slope angle, slope aspect, elevation, plan curvature, lithology, soil, land
cover, distance to roads, distance to rivers, distance to lineaments, and rainfall have been selected to assess the spatial
relationship with landslide occurrences. The predictive capability of these factors has been evaluated using the Gain Ratio
method. Using training dataset, the SVM and NBT models have been constructed to assess the susceptibility of landslide
occurrences. Finally, the performance of the SVM and NBT models have been validated and compared using receiver
operating characteristic curve technique and statistical index-based evaluations. The results show that both the SVM and
NBT models perform well for spatial prediction of landslides. Out of these, the SVM model (AUC = 0.881) outperforms
the NBT model (AUC = 0.832). Overall, SVM and NBT indicate promising methods which could be used for spatial
prediction of landslides in landslide prone areas. Moreover, the results obtained from this study could be helpful for
planning and decision making in landslide hazard management.
Table 1: Landslide affecting factors and their class utilized in this study
(a) (b)
(a) (b)
(a) (b)
Figure 6: Landslide Susceptibility Maps (LSM) using (a) the SVM model and (b) the NBT model
77
Acknowledgement
References
Guha-Sapir, D., P. Hoyois and R. Below (2014). Annual Pham, B.T., D. Tien Bui, MB. Dholakia, I. Prakash and
disaster statistical review 2013: The numbers and H.V. Pham (2016a). A comparative study of least square
trends. Brussels: CRED. support vector machines and multiclass alternating
decision trees for spatial prediction of rainfall-induced
Guzzetti, F., P. Reichenbach, M. Cardinali, M. Galli and landslides in a tropical cyclones area. Geotechnical and
F. Ardizzone (2005). Probabilistic landslide hazard Geological Engineering 34:1-18.
assessment at the basin scale. Geomorphology, 72(1–4),
pp 272-299. Pham, B.T., D. Tien Bui, I. Prakash and M.B. Dholakia
(2016b). Rotation forest fuzzy rule-based classifier
Hassanien, A.E., M.A. Salama, J. Platos and V. Snasel ensemble for spatial prediction of landslides using GIS.
(2015). Rough local transfer function for cardiac Natural Hazards:1-31.
disorders detection using heart sounds. Logic Journal of
IGP. Pradhan, B. (2013). A comparative study on the
predictive ability of the decision tree, support vector
Kantardzic, M. (2011). Data mining: Concepts, models, machine and neuro-fuzzy models in landslide
methods and algorithms. John Wiley & Sons. susceptibility mapping using GIS. Computers &
Geosciences, 51, pp 350-365.
Karegowda, A.G., A. Manjunath and M. Jayaram
(2010). Comparative study of attribute selection using Quinlan, J.R. (1986). Induction of decision trees.
gain ratio and correlation based feature selection. Machine Learning, 1(1), 81-106.
International Journal of Information Technology and
Knowledge Management, 2(2), pp 271-277. Salama, M.A., O.S. Soliman, I. Maglogiannis, A.E.
Hassanien and A.A. Fahmy (2013). Rough set-based
Kavzoglu, T., E.K. Sahin and I. Colkesen (2014). identification of heart valve diseases using heart sounds.
Landslide susceptibility mapping using GIS-based Rough Sets and Intelligent Systems-Professor Zdzisław
multi-criteria decision analysis, support vector Pawlak in Memoriam. Springer, pp 475-491.
machines and logistic regression. Landslides, 11(3), pp
425-439. Sarkar, S., D.P. Kanungo and G.S. Mehrotra (1995).
Landslide hazard zonation: A case study in Garhwal
Kohavi, R., 1996. Scaling Up the Accuracy of Naive- Himalaya, India. Mountain Research and Development,
Bayes Classifiers: A Decision-Tree Hybrid, KDD. 15(4), pp 301-309.
Citeseer, pp. 202-207.
Sengupta, A., S. Gupta and K. Anbarasu (2010).
Lucà, F., M. Conforti and G. Robustelli (2011). Rainfall thresholds for the initiation of landslide at
Comparison of GIS-based gullying susceptibility Lanta Khola in north Sikkim, India. Natural hazards,
mapping using bivariate and multivariate statistics. 52(1), pp 31-42.
Northern Calabria, South Italy, Geomorphology,
134(3–4), pp 297-308. Song, Y., J. Gong, S. Gao, D. Wang, T. Cui, Y. Li and
B. Wei (2012). Susceptibility assessment of earthquake-
NCEP (2014). Global weather data for SWAT. induced landslides using Bayesian network: A case
http://globalweather.tamu.edu/home. study in Beichuan, China. Computers & Geosciences,
42(0), pp 189-199.
Onagh, M., V. Kumra and P.K. Rai (2012). Landslide
susceptibility mapping in a part of Uttarkashi district Soni, J., U. Ansari, D. Sharma and S. Soni (2011).
(India) by multiple linear regression method. Predictive data mining for medical diagnosis: An
International Journal of Geology, Earth and overview of heart disease prediction. International
Environmental Sciences. Journal of Computer Applications, 17(8), pp 43-48.
Pham, B.T., D. Tien Bui, P. Indra and M.B. Dholakia Tien Bui, D., B. Pradhan, O. Lofman and I. Revhaug
(2015a). Landslide susceptibility assessment at a part of (2012). Landslide susceptibility assessment in Vietnam
Uttarakhand Himalaya, India using GIS – based using support vector machines, decision tree and naive
statistical approach of frequency ratio method. Bayes models. Mathematical Problems in Engineering.
International Journal of Engineering Research &
Technology 4(11), pp 338-344. Tien Bui, D., T.A. Tuan, H. Klempe, B. Pradhan and I.
Revhaug (2015). Spatial prediction models for shallow
Pham, B.T., D. Tien Bui, H.R. Pourghasemi, P. Indra landslide hazards: A comparative assessment of the
and M.B. Dholakia (2015b). Landslide susceptibility efficacy of support vector machines, artificial neural
assesssment in the Uttarakhand area (India) using GIS: networks, kernel logistic regression and logistic model
A comparison study of prediction capability of naïve tree. Landslides, pp 1-18.
79
Abstract: Satellite based remote sensing (RS) data at different spatial and temporal scales can provide crop sown area
estimates needed by decision makers. Due to conflict of spatial resolution versus temporal frequency of data collection,
getting early and accurate crop sown area at large scale is very difficult. This paper presents a methodology of early
estimation of crop sown area at large scale by making use of high temporal coarse spatial resolution data and low temporal
fine spatial resolution data. It also uses previous years’ data for extracting a-priori knowledge of crop sowing area. Early
crop area estimate was made for Gujarat state (India) for 2011-12 rabi season. Multi-date MODIS (MODerate resolution
Imaging Spectroradiometer) data and two-date Resourcesat-2 AWiFS (Advance Wide Field Sensor) data upto mid-
December were used for crop sown area early estimates. Multi-date MODIS data for previous five years provided a-priori
information on crop presence / absence over the previous five crop seasons. While ISODATA was used for classifying
multi-date MODIS and AWiFS data; hierarchical decision tree approach was used for integrating multi-source
information. Incorporating two date AWiFS data and a-priori information with multi-date MODIS data increased crop
sown area early estimates accuracy significantly.
(a) (c)
Sig. Ext. ISODATA
Mask
GT
Spectral
Pattern
(b) (d) Ref. Matching Spectral
Sign. Clusters
Rationale –
Matching of Temporal Spectral
Profiles of Clusters With
Figure 1: Examples of Harmonic ANalysis of Time Reference Temporal Spectral
Series (HANTS) applications: (a) scaled NDVI Profiles Prepared from GT
profile of a rabi crop with thin cloudy atmosphere, Classified Image Crop Layer
(b) scaled NDVI profile of another rabi crop with Figure 2: Conceptual flow of ISODATA clustering
missing data due to thick cloudy atmosphere; (c) and spectral pattern matching showing major data
after applying HANTS over profiles in (a); and (d) analysis components
after applying HANTS over profiles shown in (b)
Multi-date full-season MODIS NDVI time series data
3.3 ISODATA Classification from 2006-07 to 2010-11 were subjected to data
The Iterative Self-Organizing Data Analysis Technique screening and extrapolation before classification. The
(ISODATA) was used as an unsupervised classifier by HANTS corrected data were then classified using
recognizing multi-temporal patterns in the dataset (Ball ISODATA clustering. The temporal spectral profiles of
and Hall, 1965). It is an iterative and heuristic procedure ISODATA clusters were matched with the reference
which assigns first an arbitrary initial cluster vector temporal spectral profiles and the clusters were labelled
based on user’s input. In the second step, it classifies to rabi-crop and other than rabi classes. Thus classified
each pixel of the data to the closest cluster in spectral images were obtained for 5 rabi seasons i.e. 2006-07 to
domain. Merging and splitting of clusters is done, if 2010-11.
conditions are met. Clusters are merged if either the
number of pixels in a cluster is less than a certain Two sets of multi-date AWiFS data, as described
threshold or if the centres of two clusters are closer than earlier, were processed to prepare classified images. Set
a certain threshold. Similarly, clusters are split into two 1 was classified by subjecting 7 dates of NDVI stack to
different clusters if the cluster standard deviation ISODATA clustering and then labelling of clusters to
exceeds a predefined/user defined value and the number rabi-crop and other classes. In case of Set 2, firstly
of pixels is twice the threshold for the minimum number December 11 AWiFS data was classified by ISODATA
of pixels in a cluster. In the third step the new cluster in vegetation (crops, forest, shrubs, plantations etc.) and
mean vectors are calculated based on all the pixels in other classes and ΔNDVI (i.e. NDVIDec11 – NDVINov11)
that cluster. The second and third steps are repeated until image was used for classifying vegetation class in 4 sub-
the "change" between two consecutive iterations is classes based on NDVI gradient.
small. ISODATA has been found very effective at
identifying spectral clusters in data. It is especially very 3.4 Integration by decision tree approach
useful while analysing a new data as we don’t need to Decision tree approach is a commonly used method of
know much about the data beforehand. Care has to be information extraction in data mining. The objective is
taken that the data is structured well otherwise to create a model that estimates the value of a target
ISODATA may take long time if data is largely variable based on several input variables available from
unstructured. Conceptual flow of the procedure single or multiple data sources. A decision tree is a flow-
followed in ISODATA clustering and spectral matching chart-like structure, where each internal (non-leaf) node
is shown in Fig. 2. denotes a test on an attribute, each branch represents the
outcome of a test, and each leaf (or terminal) node holds
Before classifying 9-dates MODIS NDVI time-series a class label. The topmost node in a tree is the root node
data of 2011-12, temporal patterns of reference crops (Quinlan, 1993).
were obtained from 9-dates NDVI values at the crop
locations collected during the in-season field surveys. In this study three input images namely (i) a-priori crop
Stacked 9-dates NDVI data was subjected to ISODATA history image from 5-years multi-date MODIS data, (ii)
and clusters were obtained for 2011-12 season. The 9-dates MODIS derived classified image, (iii) 2-dates
temporal NDVI patterns of the ISODATA clusters were
84
Vermote, E.F., S.Y. Kotchenova and J.P. Ray (2011). Wen Jun, Su Zhongbo, and Ma Yaoming (2004).
MODIS surface reflectance user’s guide (ver 1.3). Reconstruction of a Cloud-free Vegetation Index Time
MODIS land surface reflectance science computing Series for the Tibetan Plateau. Mountain Research and
facility. http://modis- Development 24(4):348-353, doi:
sr.ltdri.org/guide/MOD09_UserGuide_v1_3.pdf (last http://dx.doi.org/10.1659/0276-4741(2004).
accessed April 1, 2016).
Xiangming Xiao, S. Boles, S. Frolking, C. Li, J.Y.
Vyas, S., R. Nigam, N.K. Patel and S. Panigrahy (2013). Babu, W. Salas and B. Moore III (2006). Mapping
Extracting regional pattern of wheat sowing dates using paddy rice agriculture in South and Southeast Asia
multispectral and high temporal observations from using multi-temporal MODIS images. Remote Sensing
Indian geostationary satellite. J Indian Soc Remote of Environment, 100: 95 – 113.
Sens, 41(4):855–864, DOI 10.1007/s12524-013-0266-
3. Xu, Yongming and Yan, Shen (2013). Reconstruction
of the land surface temperature time series using
Wardlow, B.D., S.L. Egbert and J.H. Kastens (2007). harmonic analysis. Computers & Geosciences 61:126–
Analysis of time-series MODIS 250 m vegetation index 132.
89
Abstract: Agriculture plays major role in India’s economy, and provides undoubtedly the largest livelihood with its
allied sectors. Crop type identification serves in number of applications such as crop yield forecasting, collecting crop
production statistics, facilitating crop rotation records, mapping soil productivity, identification of factors influencing
crop stress, assessment of crop damage and monitoring farming activity. To identify particular crop type in a single date
imagery is a challenging task. However, Classification facilitates the multi-temporal images by taking into account
changes in reflectance as a function of plant phenology. This research work deals with Possibilistic c-means classifier
with Hypertangent kernel for wheat (Triticumaestivum) identification in Haridwar, Uttarakhand, India. The vegetation
index outputs of Formosat-2 and Landsat-8 (Operational Land Imager) sensors were arranged in chronological order of
their date and prepared three temporal datasets which cover whole phenological cycle of wheat. It was evaluated that for
2.7, 2.5, and 2.5 values of weighted constant (m), images of 4 (Formosat-2: 04 Dec 2014, 30 Jan 2015, 21 Feb 2015;
Landsat-8: 16 Mar 2015), 5 (Formosat-2: 04 Dec 2014, 30 Jan 2015, 21 Feb 2015; Landsat-8: 16 Mar 2015, 01 Apr
2015), and 6 (Formosat-2: 04 Dec 2014, 30 Jan 2015, 21 Feb 2015, 09 Apr 2015; Landsat-8: 16 Mar 2015, 01 Apr 2015)
date combination respectively represent the nicely separated wheat crop from other vegetation and were easily
differentiated between early harvested and late harvested wheat crops. This study demonstrates that 5 date combination
was sufficient to discriminate late harvested wheat crops and 6 date combination was sufficient to discriminate early
harvested wheat crops.
Keywords: Wheat identification, Phenology, Soft classification, Possibilistic c-means (PCM), Kernel, Weighted constant
The major objective of this research is to propose a soft where, 𝑑𝑖𝑗 represents the distance between the pixels
classifier algorithm which has the capability to extract value i and mean of the class j.
wheat crop with better accuracy dealing with non-
linearity within the classes. Specific objectives include: PCM is robust to handle mixed pixels but it fails to
(1) To implement Hypertangent kernel based correctly classify the pixels when the classes are linearly
Possibilistic classifier for specific crop identification in non-separable (Wu, 2006). Kernel methods provide a
bi-sensor multi-spectral data, and (2) To evaluate compatible and reliable framework for developing
number of temporal images and optimized value of nonlinear technique of classification and have useful
weighted constant that may be best suited for wheat crop properties when dealing with low number of training
identification. data, presence of heterogeneous land cover and different
noise sources in the data.
2. Temporal vegetation index and classification
approach In this paper a robust supervised fuzzy classification
technique, Kernel based Possibilistic c- Means
The Normalized difference vegetation index (NDVI) algorithm (KPCM) has been presented. Its basic idea is
layers of different date were prepared and stacked in to transform the low dimensional input data into a
chronological order and prepared three temporal higher dimensional feature space via a kernel method.
datasets of vegetation index. The NDVI was proposed After the implementation of kernels, the classes become
by Kriegler et al. (1969) and the mathematical linearly separable and PCM is performed on the feature
expression is given by equation (1). space.
91
Reference data for wheat was identified on the imagery Supervised Soft Classification using
from the GPS data collected during a field visit on 16 Hypertangent Kernel based Possibilistic c-
Means Classifier
March 2015. During field survey the early sowing and
late sowing wheat crop samples were collected based on
the information provided by concern farmers. Further,
Evaluation for optimized weighted
this reference data was used as a training and testing constant and best temporal date
data for classification and validation. combination
4. Methodology
Figure 2: Methodology adopted
Four temporal data of Formosat-2 from 04 Dec 2014 to
09 April 2015 and two temporal data of Landsat-8 (OLI) 255
of 16 March 2015 and 01 April 2015, sensors have been 205
NDVI (on 8 bit scale)
0.1
0.05
0
1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3
0.2
0.1
0
1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 3
Weighted Constant (m)
Difference between Early Harvested and Late
Harvested classes in Early Harvested test sites
Difference between Late Harvested and Early
Harvested classes in Late Havested test sites
Figure 5: Classified images extracted by using Hyper
Tangent Kernel based PCM classifier, (a) Class1-
(c) Early harvested wheat crop (Six date; optimized m=
0.3
Difference in membership
0.1 6. Conclusion
value
Jensen, J.R. (2009). Remote sensing of the environment. Wardlow, B.D., S.L. Egbert and J.H. Kastens (2007).
2nd Edition. Pearson Education, pp. 382-399. Analysis of time-series MODIS 250 m vegetation index
data for crop classification in the US Central Great
Kaur, P., P. Gupta and P. Sharma (2012). Review and Plains. Remote Sensing of Environment, 108(3), pp.
comparison of kernel based fuzzy image segmentation 290-310.
techniques. International Journal of Intelligent Systems
and Applications (IJISA), 4(7), pp. 50. Wu, X.H. (2006). A possibilistic c-means clustering
algorithm based on kernel methods. In
Kriegler, F.J., W.A. Malila, R.F. Nalepka and W. Communications, Circuits and Systems Proceedings,
Richardson (1969). Preprocessing transformations and 2006 International Conference on (Vol. 3, pp. 2062-
their effects on multispectral recognition. In Remote 2066). IEEE.
Sensing of Environment, 4(1), pp. 97.
Zadeh, L.A. (1965). Fuzzy sets. Information and
control, 8(3), pp. 338-353.
95
-------------------------------------------------
In addition to students and professionals in the field who will get a broad account of the functioning of
space based camera systems and the nuances in the design, development and deployment of them, those
in policy making and technical management in space agencies across the globe will also find the book
useful to understand the path taken by India to achieve pre‐eminence in this field.
“I hope the book will interest a broad range of readers both within and outside the country.
Those who were part of the journey will feel a sense of satisfaction and proud of what they
could achieve and the younger readers will be inspired and encouraged to be part of this
excitement. The book should interest all those who want to know how India has achieved
preeminence in space based remote sensing.”
Author: Dr. George Joseph was director of the Space Applications Centre (ISRO), Ahmedabad. Under his
overall guidance the development of electro‐optical sensors started in ISRO.
96
Abstract: Due to unemployment, the crime is increasing. Criminal acts like murder, rape, kidnapping, home breaking,
theft, robbery are prohibited but exists all around the world. India is one major country where crime is increasing. The
present study deals with the crime mapping of Ajmer city with GIS approach. The study is having objectives to generate
crime maps to identify the crime pattern of Ajmer city. It is based on time series analysis to identify crime direction and
hotspots of crimes. It also analyses type of hotspot, proximity of crimes to police stations, displacement of crime across
time, crime rate of each ward and the socio-economic characteristics of city. This analysis provides insight to police with
a view to decreasing crime rate.
Key Words: Crime rate, GIS, Hotspot, Land use/land cover, Crime map
Jaishankar et al. (2004) showed that use of GIS provides Ajmer is the district of Rajasthan state. Ajmer city is the
a convenient tool for crime pattern analysis due to its head quarter of the district. Ajmer city has a population
geographic referencing capabilities. It provides valuable of 542321. Ajmer city is a religious and tourist place. It
information concerning property of crimes including is surrounded by NH-8, NH-14, NH-79 and NH-89.
data on the social and physical characteristics of areas Ajmer city is surrounded by the Aravali hills. Ajmer city
that contribute to localized criminal activity. stretches from 26˚23’ North to 26˚23’ North and 74˚36’
Thangavelu et al. (2013) discussed about the importance East to 74˚40’ east. It has nine Police Stations (PS) and
of GIS, as it can be used as a tool to identify factors one Mahila PS. All nine PS boundaries of Ajmer city
contributing to crime and thus allow police department have been selected for study. Figure 1 shows the
to proactively respond to the situations before they location map of study area.
become problematic. Crime analysis mapping is a
valuable problem solving tool because it can lead to the
identification of new problems facing law enforcement,
lend a visual perspective to an analysis, assist in the
development of an effective response, aid in the
formation of partnerships by providing a common point
of reference and assist evaluation procedures (Velasco
and Boba, 2005).
(a) (b)
(c)
(d)
Figure 5: Crime map (2009 – 2014) of four Police Stations (PS) viz., (a) Aadarsh Nagar, (b) Alwar gate, (c ) Ganj
and (d) Christianganj
100
Abstract: In the present study, a comparison of accuracy and precision between two formulae for assessing Helmert's
orthometric correction was carried out. Two test levelling lines were used for comparison. For the two lines, the actual
levelling and gravity data were used as reference (real world) data, whereas respective UCPH2002 derived gravity values
and SRTM30-interpolated elevation data were executed as erroneous input information. Practically, one formula proved
to be less sensitive to the propagation of biases and random errors of the input data, thus yielding a remarkably accurate
and precise orthometric correction. So, it is recommended to use this formula, especially when using interpolated end
benchmarks elevations in the computation of the height correction. Furthermore, the determination of realistic
uncertainties of the height corrections could help model them as observed parameters during the adjustment of levelling
loops.
gB OC g
the mean gravity relevant to B B
1 , (4f)
0
g B
g B 0.0424H B , H B
g A & g B the observed gravity at A and B, respectively. g A g A 0.0424 H A , (4g)
g B g B 0.0424 H B . (4h)
Equivalently, but alternatively expressed, the OC could
be formulated as follows (Hwang and Hsiao, 2003) So, Eq. (3) may be re-written as follows
k g g g
n i
k k
H
OC i 1 n i A 1 H A . (2) OC g i i 1n i A g A
g B
g 0 i 1 0 0
i 1 B i 1
2 2
1 k
Accordingly, Eq. (6) may be formulated as follows
2 g i ni H A g A 2 g B gg A 1 2H A .
k
n i g H k
B g i 1
OC g i i 1 n i A g A
B
(12)
i 1 gB i 1 gB gB
1 k g Also, using the law of error propagation, it follows that
2
g i n i H A g A g B A 1 H A .
g B i 1 gB 2 g A 2 g A (0.0424)2 2 H A , (13a)
(8)
The impact of random errors on OC follows from the 3. Input data and input errors for the test levelling
application of the law of variance-covariance lines
propagation to Eqs. (1) and (2). So, neglecting the error
covariances and applying the error propagation The current study considers two levelling lines in Egypt
principle to Eq. (1), the OC variance, 2OC , is with observed elevations and gravity data. The first line
expressed by (e.g. Ghilani and Wolf, 2006) (Line I) runs along the Nile valley, while the second
(Line II) is located in the Western desert. Line I and line
2 2 2 II exhibit relatively moderate and mountainous terrain
k k
2OC OC 2 g i OC 2ni OC 2 g A roughness, respectively. Figure 1 depicts two post maps
i 1 g i i 1 ni
g A for the gravity points along both lines. While the
available terrestrial gravity and elevation data are used
2 2 2 as real world reference values in the current
OC OC OC
gB
2 HA
2 HB ,
2 investigation, the input "erroneous" gravity and height
g B H A H B information (input into the error analysis algorithms)
(9) are derived from the UCPH2002 geo-potential model
and the SRTM30 terrain model, respectively. First, geo-
where g i , n i , g A , g B , H A and H B
potential model derived surface gravity was derived as
the sum of the UCPH2002-synthesized surface gravity
stand for the error standard deviations relevant to the disturbances relative to WGS-84 and the respective
input quantities, which reflects their precisions. The normal gravity values (Forsberg and Tscherning, 2008;
partial derivatives in Eq. (9) are the same derived in Eq. Filmer and Featherstone, 2011; Filmer et al., 2013;
(4). So, Eq. (9) can be written as follows Hassouna, 2013). On the other hand, the SRTM30 grid
was used to interpolate respective elevation data at the
2
k n 2 k g
2 points of lines I and II, via the B-spline interpolation
2OC i 2 g i i 1 2ni H A 2 g A algorithm (Cimmery, 2010).
i 1 0 i 1 0 0
Table 1 lists the different features of the two levelling
2 2 2 lines under consideration. Table 2 shows the statistics
HB
2 g B g A 1 2H A g B 1 2H B . of the differences among the derived "erroneous"
0 0
0
gravity and benchmark elevation values; and the
(10) respective reference values over the investigated lines.
g A
& g B were appropriately computed from
Figure 1: Post maps for Lines (I) and (II) Eqs. (4g & 4h) and (13), respectively, taking into
account whether or not the gravity or benchmark
elevations are kept errorless.
Table 1: Different features of the two levelling
profiles In the above sense, Tables 3a and 3b summarize the
accuracies of OC, as computed from Eq. (5) and Eq. (8).
No. of Mean Also, Tables 4a and 4b show the respective items, but
H Length regarding the precisions of OC determination obtained
Line gravity elevation
(m ) (km ) via Eqs. (10) and (12).
sections (m )
I 93 -24.68 34.85 280.4
5. Discussion and concluding remarks
II 78 -391.82 268.18 369.2
Tables 3a and 3b show that in general, the accuracies of
The observational noises of the observed (reference) OC determination by the two formulae deteriorate as
gravity data for the test lines are in the order of 1 mgal. the terrain roughness increases. Table 3a implies that
Also, the error budget for the reference elevations of the Eqs. (1) and (2) yield equally accurate OC values, if one
start and end benchmarks could be as worse as a few deals with highly accurate (i.e. nearly errorless)
centimeters. However, such small uncertainties might benchmark elevations. Table 3a shows also that in such
not represent those of any eventual interpolated gravity case, geo-potential models derived gravity could safely
or elevation values, which are greatly affected by the be used to assess quite accurate OC values in case of
local gravity and terrain signal roughness. So, it was moderate terrain roughness, as that represented by line
decided to use the spatial standard deviations of the I. This is easy to conclude from the respective two zero
106
Table 2: Statistics of the discrepancies among the erroneous and reference quantities for the two lines
Max .
g
Min .
Mean
g g
Line Item Unit
g
i g i
g i
A B
Table 3a: Comparison among the OC accuracies Tables 4a and 4b. So, Eq. (2) is more precise than Eq.
(gravity accuracy emphasized) (1), when interpolated benchmark elevations are used.
gi gA gB OC OC So, generally, for a high quality OC determination, it is
Line ( mm ) ( mm ) recommended to use Eq. (2). Such results should be
( mgal ) ( mgal ) ( mgal ) understood to hold to the investigated interpolated
(Eq. 5) (Eq. 8) benchmark elevations and synthesized gravity values
I Section- 11.50 20.78 0.0 0.0 along levelling lines. Besides being more economic,
wise such levelling data sources may be the only available
II varying 10.35 17.95 -4.6 -4.6
tool for assessing OC. The use of the SRTM30 and the
Table 3b: Comparison among the OC accuracies UCPH2002 low resolution models in the current study
(elevation accuracy emphasized) was for the sake of handling an instance worst case. A
HA HB OC OC future work may test other high resolution DTMs and
Line ( mm ) ( mm ) geo-potential models for representing erroneous data.
(m ) (m )
(Eq. 5) (Eq. 8)
A further related criterion could be the assessment of the
I -3.34 -2.37 2.0 0.6
values and qualities of OCs for height differences along
II 14.94 -8.50 -40.2 -4.7 profiles, which are extracted from DEMs. Such terrain
models may be derived from satellite imagery or
Table 4a: Comparison among the OC precisions photogrammetry. In this respect, the published
(gravity precision emphasized) uncertainty of elevations would often be of the order of
gi g g OC OC that discussed in the current study. Accordingly, Eq. (2)
Line A B
( mm ) ( mm ) is expected to show better error propagation
( mgal ) ( mgal ) ( mgal ) characteristics.
(Eq. 10) (Eq. 12)
I 8.33 0.5 0.5 Finally, it could seem rational and innovative to model
II 13.17 8.3 8.3 the OCs as biases over levelling circuits in vertical
networks. In this sense, OCs over levelling lines might
Table 4b: Comparison among the OC precisions be treated as observed weighted parameters during the
(elevation precision emphasized) adjustment of levelling nets (e.g. Ghilani and Wolf,
2006). Of course, the weights of such parameters should
H and H B
OC OC lean on the respective error standard deviations, and
hence, the importance of OC precision may arise. This
A
Line ( mm ) ( mm )
(m ) in turn could optimally help pick out the random loop
(Eq.10) (Eq. 12) closures parts during the adjustment.
I 1.51 3.1 0.3
II 32.90 78.7 11.3 Acknowledgment
This award has been instituted to recognize outstanding and conspicuously important contribution in promoting geomatics
technology and applications at the country level. The contributions should have made major impact on the use of this
technology for national development.
The award shall consist of Rs. 50,000/- in cash, a medal and citation.
Eligibility
Any citizen of India, engaged in activities related to geomatics technology and its applications is eligible for this award. The
prize is awarded on the basis of work primarily done in India.
The age limit for awardees is 45 years or above as on June 30 of the year of award.
Selection
A duly constituted Award Committee will evaluate all nominations received. The committee shall consist of eminent experts
in the field of geo-spatial technology, to be identified by the Executive Council, ISG. The committee shall forward selected
name/s to ISG – EC for approval and announcement. Apart from those persons, whose nominations have been received, the
Committee may consider any person or persons who, in their opinion, have made outstanding contributions to development
of geo-spatial technology and applications.
The award can be withheld in any year if, in the opinion of the committee, no candidate is found suitable in that particular
year.
The award shall be presented during the Annual Convention of ISG. Local Hospitality shall be taken care by ISG & Air fare
(low cost) may be reimbursed if awardees request for it.
The nominations can be proposed by Head of a major research institute/ centre; Vice-Chancellor of a university; Secretary
of Government Scientific Departments; President of a National Academy, President, Indian Society of Geomatics / Indian
Society of Remote Sensing / Indian National Cartographic Association / ISG fellow or two life members of the society with
more than 10 year old membership.
A candidate once nominated would be considered for a total period of two years. Nomination should be sent in the prescribed
format to Secretary, ISG.
The last date for receiving nominations shall be September, 31 or otherwise extended.
vi
Proposer:
Signature
Name
Address
Phone/ Fax
E-mail
Life Membership No. (in case of ISG Member):
Signature
Name
Address
Phone/ Fax
E-mail
Life Membership No. (in case of ISG Member):
3. Photogrammetry
4. Digital Cartography
Eligibility
Any citizen of India engaged in scientific work in any of the above-mentioned areas of research is eligible for the award.
The awards are to be given for the work largely carried out in India.
First award will be given for original contribution in the field of Geomatics supported by publications in a refereed
journal of repute.
The contribution for the first award should have been accepted by peers through citation of the work.
Work based on the applications of existing technologies will not be considered for the first award.
The work should have made impact on the overall development of Geomatics.
Nominations should be sent in the prescribed format, completed in all aspects to the Secretary, Indian Society of Geomatics,
Space Applications Centre Campus, Ahmedabad 380 015 by August 31, 2016.
Selection Process
An expert committee, consisting of at least three members, constituted by the Executive Council of the Indian Society of
Geomatics, will scrutinize the nominations and recommend the awardees’ names to the Executive Council. The Council will
decide on the award based on the recommendations.
ix
The Applications in the above format (five copies) should be submitted (by Registered Post or Speed Post) to
*ISG is, however, not bound to accept these names and can refer the nomination to other experts/peers
x
P-1 Director, Space Applications Centre (ISRO), Jodhpur Tekra Satellite Road, Ahmedabad - 380 015
P-2 Settlement Commissioner, The Settlement Commissioner & Director of Land Records-Gujarat, Block No. 13, Floor 2, Old Sachivalay, Sector-10,
Gandhinagar – 382 010
P-3 Commissioner, Mumbai Metropolitan Region Development Authority, Bandra-Kurla Complex, Bandra East, Mumbai - 400 051
P-4 Commissioner, land Records & Settlements Office, MP, Gwalior - 474 007
P-5 Director General, Centre for Development of Advanced Computing (C-DAC), Pune University Campus, Ganesh Khind, Pune - 411 007
P-6 Chairman, Indian Space Research Organization (ISRO), ISRO H.Q., Antariksha Bhavan, New BEL Road,Bengaluru 560 231
P-7 Director General, Forest Survey of India, Kaulagarh Road, P.O. I.P.E., Dehra Dun – 248 195
P-8 Commissioner, Vadodara Municipal Corporation, M.S. University, Vadodara - 390 002
P-9 Director, Centre for Environmental Planning and Technology (CEPT), Navarangpura, Ahmedabad - 380 009
P-10 Managing Director, ESRI INDIA, NIIT GIS Ltd., 8, Balaji Estate, Sudarshan Munjal Marg, Kalkaji, New Delhi - 110 019
P-11 Director, Gujarat Water Supply and Sewerage Board (GWSSB), Jalseva Bhavan, Sector – 10A, Gandhinagar - 382 010
P-12 Director, National Atlas & Thematic Mapping Organization (NATMO), Salt Lake, Kolkata - 700 064
P-13 Director of Operations, GIS Services, Genesys International Corporation Ltd., 73-A, SDF-III, SEEPZ, Andheri (E), Mumbai - 400 096
P-14 Managing Director, Speck Systems Limited, B-49, Electronics Complex, Kushiaguda, Hyderabad - 500 062
P-15 Director, Institute of Remote Sensing (IRS), Anna University, Sardar Patel Road, Chennai - 600 025
P-16 Managing Director, Tri-Geo Image Systems Ltd., 813 Nagarjuna Hills, PunjaGutta, Hyderabad - 500 082
P-17 Managing Director, Scanpoint Graphics Ltd., B/h Town Hall, Ashram Road, Ahmedabad - 380 006
P-18 Secretary General, Institute for Sustainable Development Research Studies (ISDRS), 7, Manav Ashram Colony, Goplapura Mod, Tonk
Road, Jaipur - 302 018
P-19 Commandant, Defense institute for GeoSpatial Information & Training (DIGIT), Nr. Army HQs Camp, Rao Tula Ram Marg, Cantt., New
Delhi - 110 010
P-20 Vice President, New Rolta India Ltd., Rolta Bhavan, 22nd Street, MIDC-Marol, Andheri East, Mumbai - 400 093
P-21 Director, National Remote Sensing Centre (NRSC), Deptt. of Space, Govt. of India, Balanagar, Hyderabad - 500 037
P-22 Managing Director, ERDAS India Ltd., Plot No. 7, Type-I, IE Kukatpalli, Hyderabad - 500 072
P-23 Senior Manager, Larsen & Toubro Limited, Library and Documentation Centre ECC Constr. Gp., P.B. No. 979, Mount Poonamallee
Road, Manapakkam, Chennai - 600 089.
P-24 Director, North Eastern Space Applications Centre (NE-SAC), Department of Space, Umiam, Meghalaya 793 103
P-25 Progamme Coordinator, GSDG, Centre for Development of Advanced Computing (C-DAC), Pune University Campus, Pune –
411 007
P-26 Chief Executive, Jishnu Ocean Technologies, PL-6A, Bldg. No. 6/15, Sector – 1, Khanda Colony, New Panvel (W), Navi Mumbai – 410 206
P-27 Director General, A.P. State Remote Sensing Applications Centre (APSRAC), 8th Floor, “B” Block, Swarnajayanthi Complex,
Ameerpet, Hyderabad- 500 038
P-28 Director, Advanced Data Processing Res. Institute (ADRIN), 203, Akbar Road, Tarbund, Manovikas Nagar P.O., Secunderabad –
500 009
P-29 Managing Director, LEICA Geosystems Geospatial Imaging Pvt. (I) Ltd., 3, Enkay Square, 448a Udyog Vihar, Phase-5,
Gurgoan- 122 016
P-30 Director, Defense Terrain Research Limited (DTRL), Ministry of Defense, Govt. of India, Defense Research & Development
Organisation, Metacafe House, New Delhi – 110 054
P-31 Chairman, OGC India Forum, E/701, Gokul Residency, Thakur Village, Kandivali (E), Mumbai – 400 101
P-32 Managing Director, ML Infomap Pvt. Ltd., 124-A, Katwaria Sarai, New Delhi – 110 016
P-33 Director, Rolta India Limited, Rolta Tower, “A”, Rolta Technology Park, MIDC, Andheri (E), Mumbai – 400 093
P-34 Director, State Remote Sensing Applications Centre, Aizawl – 796 012, Mizoram
xi
The journal covers all aspects of Geomatics – geodata acquisition, Normally length of a published paper should be about 6-10 pages
pre-processing, processing, analysis and publishing. Broadly this in A4 size including figures. Use of illustrations in colour should
implies inclusion of areas like GIS, GPS, Photogrammetry, be restricted and resorted to only where it is absolutely necessary
Cartography, Remote Sensing, Surveying, Spatial Data and not for enhancing the look of the paper. If the number of
Infrastructure and Technology including hardware, software, colour illustrations exceeds five, authors’ institution may be
algorithms and model. It endeavours to provide an international asked to reimburse the extra cost involved, which at current rates
forum for rapid publication of developments in the field – both in is about Rs. 2500 per coloured figure/diagram/plate/illustration.
technology and applications.
Submission of Manuscript
To begin with the frequency of publication will be six- monthly.
However, depending on the response and interest, frequency of Submissions should be in electronic form via email. The
publication may be reviewed. manuscript may be sent by email to [email protected]. In
exceptional cases hard copy submission in camera ready form may
A manuscript for publication must be based on original research be allowed with the prior permission of the Chief Editor.
work done by the author(s). It should not have been published in Submission in any other form will be returned to the author. To
part or full in any type of publication nor should it be under speed up the review process, authors are advised to provide a list
consideration for publication in any periodical. Unsolicited review of three probable reviewers with their institutional address and e-
papers will not be published. mail IDs.
The Editorial Board or the Indian Society of Geomatics is not Guidelines for Citing References
responsible for the opinions expressed by the authors.
Names of all cited publications should be given in full.
Language No abbreviations should be used. Following procedure is
to be adopted.
The language of the Journal will be English (Indian). However,
manuscripts in English (US) and English (British) are also Journal Publications
acceptable from authors from countries located outside India.
Bahuguna, I.M. and A.V. Kulkarni (2005). Application of digital
Manuscript Format elevation model and orthoimages derived from IRS-1C Pan
stereo data in monitoring variations in glacial dimensions,
Each paper should have a title, name(s) of author(s), and affiliation Journal of the Indian Society of Remote Sensing, 33(1), 107- 112.
of each of the authors with complete mailing address, e-mail (to be referred to in the text as Bahuguna and Kulkarni (2005) or
address, an abstract, four to six keywords, and the text. The text if more than two sets of authors are to be referred to, as (Bahuguna
should include introduction/background, research method, results, and Kulkarni, 2005; Jain et al., 1994)) When more than two
discussion, followed by acknowledgements and references. The authors are to be referred to, use Jain et al. (1994). However, in
main text should be divided in sections. Section headings should References, all authors are to be mentioned.
be concise and numbered in sequence, using a decimal system for
subsections. Figures, images and their captions should be inserted Publication in a Book
at appropriate points of the text. Figures, images and tables should
fit in two column format of the journal. If absolutely necessary, Misra, V.N. (1984). Climate, a factor in the rise and fall of the
figures, images and tables can spread across the two columns. Indus Civilization – Evidence from Rajasthan and Beyond in
Figures and images, however, should not exceed half a page in Frontiers of the Indus Civilization (B.B. Lal and S.P. Gupta: Chief
height. A title should be provided for each Table, Image and Editors) Books and Books, New Delhi, pp. 461-489
Figure. All figures and images should be in 600 dpi resolution and
sized as per column/margin width. Authors must ensure that Papers Published in Seminar/ Symposium
diagrams/figures should not lose easy readability upon reduction Proceedings
to column size. The SI (metric) units and international quantities
should be used throughout the paper. In case measurements are Jain, A., A.R. Shirish, M. Das, K. Das, M.C. Porwal, and P.S. Roy
given in any other system, equivalent measurements in SI (metric) (1994). Remote Sensing and Geographic Information System – An
units should be indicated in brackets. approach for the assessment of biotic interference in the forest
ecosystem. Proceedings. 15th Asian Conference on Remote
Use MS Word with English (UK/US) or English (Indian) Sensing,Bangaluru, November 17-23, 1994, pp. 65-72.
dictionary. The page size should be A4 paper, with 2 cm margin
on all sides. Title, authors and affiliation should be centred. Books
Abstract should be justified across margins. The manuscript text
should be in two columns of 8.2 cm each with a gutter of 6mm Possehl, Gregory L. (1999). Indus Age: The beginnings. Oxford
between them. Use only Times New Roman fonts. Title should be and IBH Publishing Corporation, New Delhi.
12 points bold. Authors and affiliation should be 9 points. All
other text including headings should be 10 points. Heading
numbering scheme should be decimal e.g. 1, 1.1, 1.2.3, etc.
Headings should be in bold.
xii
Journal of Geomatics Vol 10 No. 1 April 2016
Reviewing
If the authors have used any copyright material in their
Each paper will be reviewed by three peers. Papers forwarded by manuscript, it is understood that they have obtained permission
members of the Editorial or Advisory Boards along with their from the owner of the copyright material and they should convey
comments would get processed faster and may be reviewed by two the same along with the manuscript to the Chief Editor.
referees only.
Certificate of Original Work
Sample format for Authors is available in downloadable form at
ISG website: www.isgindia.org/JOG/Sample_format.doc The authors will also provide a certificate that the paper is an
original work, not published or being considered for publication
Copyright elsewhere.
The copyright of the paper selected for publication will rest with In the event the certificate turns out to be false, the Journal shall
the Indian Society of Geomatics. Corresponding author shall be ban the author(s) from publishing in the Journal for a period of
required to sign a copyright assignment form, on behalf of all five years and inform the same to all other related publications.
authors, once the paper is selected for publication. Authors are,
however, at liberty to use this material elsewhere after obtaining Reprints
permission from the Indian Society of Geomatics.
Authors will be provided soft copy (PDF) of their paper, no hard
copy reprints will be provided.
xiii
Journal of Geomatics Vol 10 No. 1 April 2016
Journal of Geomatics
Advertisement Rates
1 Issue 4 Issues
Back Cover Page in colour Rs. 25,000 Rs. 80,000
Inside Cover Page in colour Rs. 20,000 Rs. 64,000
Full Page inside in colour Rs. 15,000 Rs. 48,000
Full Page inside in B/W Rs. 10,000 Rs. 32,000
Advertisement Details
Mechanical Details
Double Spread/Center Spread (42 x 29.7) cm
Full page bleed (21 x 29.7) cm
Full page non-bleed (19 x 27.7) cm
Art Requirements
Negatives: Art must be right reading, emulsion, down. Film must be supplied in one piece per color, each identified
by color. Camera-ready art is accepted for black & White adds; however, film is preferred. Electronic Files are also
accepted.
Electronic File Requirements: All material must be received before ad close dates.
Software: Adobe illustrator 9.0 (saved as EPS). Adobe Photoshop CS (saved as EPS or TIFF). Please convert higher
versions down. If you can only supply an IBM format, the file must be in viewable EPS or TIFF format with fonts
embedded as that format.
Colour Ads: Colour separations must be provided, right reading, emulsion down. Please note that files using RGB or
Pantone colours (PMS) must be converted to CMYK before we receive files.
xiv
The Secretary
Indian Society of Geomatics
39, Basant Bahar - II
Bopal, Ahmedabad – 380058
Sir,
I want to become a Member/ Life Member/ Sustaining Member/ Patron Member/ Foreign Member/ Student Member of
the Indian Society of Geomatics, Ahmedabad for the year ______. Membership fee of Rs. ______/- is being sent to you
by Cash/DD/Multicity Cheque. (In case of DD/ Cheque No. _________, drawn on Bank
________________________________.
Date:
Place:
Signature
1. Name:___________________________________________________________________________________________
2. Address:_________________________________________________________________________________________
_____________________________________________________________________________________ PIN: __________
Phone:________________ Fax: _______________ Email: ____________________________________________________
Signature of Proposer
MEMBERSHIP SUBSCRIPTION
MEMBERSHIP GUIDELINES
1. Subscription for Life Membership is also accepted in two equal installments payable within a duration of
three months, if so desired by the applicant. In such a case, please specify that payment will be in
installments and also the probable date for the second installment (within three months of the first
installment).
2. A Member of the Society should countersign application of membership as proposer.
3. Subscription in DD or Multicity Cheque should be made out in the name of ‘INDIAN SOCIETY OF
GEOMATICS’ and payable at Ahmedabad.
4. Financial year of the Society is from April 1 to March 31.
5. For further details, contact Secretary, Indian Society of Geomatics at the address given above.
6. ISG has chapters already established at the following places. Ahmedabad, Ajmer, Bhagalpur, Bhopal,
Chennai, Dehradun, Delhi, Hyderabad, Mangalore, Mumbai, Mysore, New Delhi, Pune, Srinagar,
Tiruchirappalli, Vadodara and Visakhapatnam. Applicants for membership have the option to contact
Secretary/Chairman of the local chapter for enrolment. Details can be found at the website of the Society:
www.isgindia.org
7. Journal of the Society will be sent only to Patron Members, Sustaining Members and Life Members.