Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Make FoV routine #449

Open
wants to merge 232 commits into
base: develop
Choose a base branch
from
Open

Make FoV routine #449

wants to merge 232 commits into from

Conversation

aburrell
Copy link
Contributor

@aburrell aburrell commented Aug 6, 2021

A preliminary version of the routine that automatically determines the backscatter return direction. This partially addresses #135.

Test results using a fitacf 3.0 file:

make_fov -cn a ~/Programs/Data/SuperDARN/Downloads/19960916.0400.00.han.fitacf3 > han.tst_fov

This creates an ascii file with 13727 lines and the following columns:

#STID DATE TIME BMNUM BMAZM CPID INTT_SC INTT_US NAVE FRANG RSEP RXRISE FREQ NOISE ATTEN CHANNEL NRANG RG GFLG FOVFLG FOV_PAST GRPFLG GRPNUM GRPID P_0 P_0_ERR V V_ERR W_L W_L_ERR P_L P_L_ERR PHI0 PHI0_ERR ELV ELV_LOW ELV_HIGH VH VH_ERR VH_METHOD REGION HOP DIST MED_P_0 MED_P_0_ERR MED_V MED_V_ERR MED_W_L MED_W_L_ERR MED_P_L MED_P_L_ERR MED_PHI0 MED_PHI0_ERR OPP_ELV OPP_ELV_LOW OPP_ELV_HIGH OPP_VH OPP_VH_ERR OPP_VH_METHOD OPP_REGION OPP_HOP OPP_DIST

Re-creating one of the figures from my paper leads to similar results. This scan is from a different time that was created during various spot checks.

fov_fig4_test_third_scan

Before moving this from draft status we need to:

  • merge Find frequency bands #444
  • Move general utility functions into appropriate libraries
  • Create a fitacf-esq non-ascii output file option
  • Finalize command-line options

aburrell and others added 30 commits June 14, 2019 08:58
Commented the exclude_outofscan subroutine.
Added more descriptions to the subroutine headers
First draft of a routine to read fit binaries to stdout
Renamed file to better reflect functionality
Copied the neccessary functions from make_grid.  Need these to be somewhere common in the future.
Created routine to write out fit data in a human readable format.
Updated file to not contain the duplicate grid routines and use the new write function
Created a makefile for the new procedure
Changed tabbing.
Renamed directory to match binary
Added the basis for a median fitting routine that uses fit data.
Removed unnecessary routine created for summer student.
Added writing routine that cycles through new structure.
Created a new structure that will be needed for FoV calculation.
Added routines to support new fit structure
Fixed bug in end of line.
Moved multscan.h to the inlcude folder to follow the typical
rst library structure and then modified the makefile to allow it
to be compiled correctly.
Updated the .gitignore to ignore the local user changes needed to install
RST.
Removed lines that weren't working.
Fixed errors and warnings present during compilation.  Added comments to
help ensure proper local libararies are included in the future.
Removed trailing whitespace from makefiles.
Added some whitespace to separate testing and reading parts of routine.
Moved read_fit from library to tools directory, a more appropriate location.
Generalized the loading routine, to either write the data line-by-line or store
the data in a multi-scan linked list.
Removed unused libraries.
Removed trailing whitespace in makefile
Added the outline of the FoV routine.  To avoid adding more data structures
than necessary, included all of the data and flags that will be needed for the
backscatter ID code as well.
Added make/free routines for the multscan structures.
Adapted Simon's elevation routine to incude a field of view flag.
@pasha-ponomarenko
Copy link
Contributor

Yep, I believe that what you are supposed to do with XMLs.

@aburrell
Copy link
Contributor Author

aburrell commented Feb 3, 2022

So, opening any of the files in a web browser doesn't do anything but display the code. I guess the good news is that my xml behaves the same as the rest of the xml. None of them open.

Given that it seems none of us have the knowledge to resolve this issue, here's what you need to run the code:

Create ascii output: make_fov -ascii
Use your own tdiff: make_fov -update-tdiff -tdiff <value_in_microseconds>
Specify your own frequency band: make_fov -tfmin -tfmax
Specify channel a: make_fov -cn a

Because you need to look at the range and azimuthal variations, you need to separate frequency bands for the analysis. For best results, use the TDIFF suited to the choses radar channel and frequency!

To convert the dmap binary to ascii, use fitmultbsidtoascii

@pasha-ponomarenko
Copy link
Contributor

Yep, the XML business looks like a separate issue.
:-(
Thanks for the clarification. Do you have a description of the extra fields in a more robust format like DOC or PDF? Or simple ASCII?

@ecbland
Copy link

ecbland commented Feb 3, 2022

To view the documentation in a web browser, I think you first need to run make.doc from the root directory of RST.

cd $RSTPATH
make.doc
cd doc/html
google-chrome index.html &   # or another web browser of your choice

You can then navigate through the html documentation in the web browser.
The documentation for the make_fov binary is here: $RSTPATH/doc/html/superdarn/src.bin/tk/tool/make_fov/index.html

@aburrell
Copy link
Contributor Author

aburrell commented Feb 3, 2022 via email

@ecbland
Copy link

ecbland commented Feb 3, 2022

@aburrell Now I am confused. Here's a screenshot from following the instructions I gave in my previous comment. Is this what @pasha-ponomarenko was looking for?

image

@aburrell
Copy link
Contributor Author

aburrell commented Feb 3, 2022

Is this what @pasha-ponomarenko was looking for?

Yes, @ecbland that is it!

@pasha-ponomarenko
Copy link
Contributor

@ecbland, thanks a lot, I can see this page after I run make.doc
;-)
I was able to process a FIATCF file last night, and I wanted to figure out what is *.fov data format so that I would be able to visualise it one way or another. From dmapdump output I saw that it is similar to FTITACF but with extra fields.

@pasha-ponomarenko
Copy link
Contributor

@aburrell, I might miss something, but there seem to be fewer fields in a binary output of make_fov than in its ASCII one, e.g., MED_* and GRPFLG, GRPNUM, GRPID are all missing.

@aburrell
Copy link
Contributor Author

aburrell commented Feb 10, 2022 via email

@pasha-ponomarenko
Copy link
Contributor

pasha-ponomarenko commented Feb 10, 2022

OK, here are some preliminaries.
I used 20150223.2001.00.sas.fitacf (v3, twofrequency mode) which contains both ionospheric and ground scatter. I also selected SAS data because the ground scatter here is normally dominated by signals received through the grating lobe. Here are plots for the 10 MHz band of the twofrequency mode. Importantly, for this interval there was no detectable tdiff offset for this frequency band. The data show both ionospheric and ground scatter, and ground scatter coming from rather high elevation angles which is indicative of the "back" lobe scatter.

image

image

I processed this data for frequencies between 10 and 11 MHZ using following command line options:
make_fov -vb -tfmin 10000 -tfmax 11000 20150223.2001.00.sas.fitacf > 20150223.2001.00.10.sas.fov
Then I converted DMAP *.fov file into NetCDF with following steps:

  • Analysed DMAP structure and produced supplementary files:
    dmaptocdl 20150223.2001.00.10.sas.fov fov.cdl fov.cdlmap
  • Adjusted the field lengths in the "mask" file fov.cdl to 75.
  • Generated an empty NetCDF file
    ncgen fov.cdl -o fov_10.netcdf
  • Filled the empty file with data from DMAP file:
    dmaptoncdf 20150223.2001.00.10.sas.fov fov.cdlmap fov_10.netcdf

After that used IDL to read and plot NetCDF FOV data:

Velocity
image
Elevation
image
FOV
image
GSCT
image
SCT
image

@pasha-ponomarenko
Copy link
Contributor

pasha-ponomarenko commented Feb 10, 2022

My understanding is that the majority of the 1.5 ionospheric scatter above rage gate 40 is removed due to one of the selection criteria (exceeding the maximum virtual height?).
The FOV determination seems to be a bit shaky for the majority of the ground scatter before 21:30 UT. I guess it might be related to an apparent partial overlap between the front- and "back"-lobe ground scatter. Indeed, the closer edge of the GS band between gates 18 and 25 shows very high elevation angle values in excess of 35 deg followed by a sharp boundary with elevation angle values dropping to 20-25 deg at gates ~25-32.

@aburrell, please correct me if I am wrong here.
BTW, I am not sure what SCT flag means.

@ecbland ecbland modified the milestones: RST 4.7, RST 5.0 Mar 14, 2022
aburrell and others added 4 commits April 7, 2022 16:05
Made changes needed to use the new channel-dependent tdiff.
Update the parameter description for the lobe dependent elevation calculation.
Added a propagation path warning for the users.
@aburrell
Copy link
Contributor Author

aburrell commented Apr 8, 2022

Ok, I played around with things and the main reason some things aren't working well is the propagation grouping. I applied the HAN defaults (basically the same, but -fhmax 900) and got much more of the ionospheric scatter back. Fine tuning the propagation virtual height ranges would further improve the results. I added a note to the docs to make sure users know how important this is, and know where the defaults come from.

Looking at the azimuthal variations, it makes sense that much of the groundscatter returns from the front FoV, since the variations are consistent across the beams.

sas_bm15_fhmax900_test

@aburrell
Copy link
Contributor Author

aburrell commented Apr 8, 2022

BTW, I am not sure what SCT flag means.

It's the sct parameter from the fitacf structure, I believe it just flags whether the scatter is valid or not.

@pasha-ponomarenko
Copy link
Contributor

It's the sct parameter from the fitacf structure, I believe it just flags whether the scatter is valid or not.

I don think there is such a parameter in the fitacf structure
https://radar-software-toolkit-rst.readthedocs.io/en/latest/references/general/fitacf/,
that is why I asked
:-P
It does not look like a quality flag...

@pasha-ponomarenko
Copy link
Contributor

@aburrell , What language did you make your plots in? Python? IDL?
Is there a way to integrate some plotting into this package (or utilise some existing RST plotting routines) in order to make tester's/user's life a bit easier?
:-{

@aburrell
Copy link
Contributor Author

I don think there is such a parameter in the fitacf structure

There is! Just grep "sct" in $RSTPATH/include/superdarn/*.h and you'll see which structures it's in. You can also examine the code to see how it's used. It may not be included in the user-description of the output, but it's a code structure parameter.

What language did you make your plots in

Python. I don't know how to plot in C, so I won't be adding this to RST. They're just RTI plots like the ones you made, but further separated by hop or propagation path. The plots you made were fine, they confirmed the code is working as it should be!

Added the IDL tarball created when computers without IDL compile RST to the list of ignored files.
@ksterne ksterne modified the milestones: RST 5.0, RST 5.1 Jul 6, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

Adding FoV detection, improved groundscatter flag
7 participants