This code is an implementation of the Step filtering and the Polynomial filtering method proposed in 
providing an enhancement on partial curvilinear structures.
The method was applied on geophysical images in an effort to recognize the linear patterns
of subsurface architectural structures that exist in archaeological sites .
This method efficiently combines a rotation
and scale-invariant filter and a pixel-labeling method, providing
a robust enhancement and detection of mostly line structures in
2-D grayscale images, respectively. This is the code for the rotation and scale-invariant filter , .
You can find more details in
Usage: runFiltering.m : An example to run the function step_pol_filtering (Step filtering and Polynomial filtering)
data.mat : A synthetic magnetic image (see ) used for testing
dataR.mat: A real magnetic image (see ) that can be used for testing
It returns the filter output that provides an enhancement on partial curvilinear structures.
We will appreciate if you cite our paper in your work.
 C. Panagiotakis, E. Kokinou and A. Sarris, Curvilinear Structure Enhancement and Detection in Geophysical
Images, IEEE Trans. on Geoscience and Remote Sensing, vol. 49, no. 6, pp. 2040-2048, 2011.
 C. Panagiotakis and E. Kokinou, Linear Pattern Detection of Sea Faults via a Topology and Shape Optimization Method, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, vol. 8, no. 1, pp. 3-11, 2015.
The results are plotted in Grayscale