kristinemlarson / gnssrefl

GNSS Interferometric Reflectometry Software (GNSS-IR)
GNU General Public License v3.0
154 stars 76 forks source link

Unable to generate SNR #279

Closed zihelyu closed 5 months ago

zihelyu commented 5 months ago

The program I installed using docker cannot generate SNR after running it. Is there something wrong with the command I ran? Thank you for answering this for me.

iShot_2024-05-22_22 13 38 iShot_2024-05-22_22 13 49
kristinemlarson commented 5 months ago

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

zihelyu commented 5 months ago

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

Already uploaded to this link http://kh.dawnwl.com/zjgh2370.20o

kristinemlarson commented 5 months ago

i know for sure that my code will not allow 30 observable types. there is no reason to have doppler (all the observations starting with D). i will try to change the code so that it gives you a better warning about this.

in the mean time, try this command.

smoosh zjgh2370.20o 15

but run it where the file exists. i.e. ls zjgh* should be where you run the command.

this will remove all the extra observables and only keep SNR. and the 15 decimates so current sampling, which is 15. and then try rinex2snr.

On Wed, May 22, 2024 at 8:32 AM zihelyu @.***> wrote:

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

Already uploaded to this link http://kh.dawnwl.com/zjgh2370.20o

— Reply to this email directly, view it on GitHub https://github.com/kristinemlarson/gnssrefl/issues/279#issuecomment-2125094749, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFSDAPCRWV5P3ZX4KC5HFNLZDS27BAVCNFSM6AAAAABIDX4QJCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCMRVGA4TINZUHE . You are receiving this because you commented.Message ID: @.***>

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

kristinemlarson commented 5 months ago

I investigated a bit further and there is something about your file that does not like the line where I tested for RINEX version number:

rinext =np.loadtxt(r,usecols=0,dtype='str',max_rows=1)

it works for RINEX files from typical GNSS archives.

I am not sure why it does this - and I could provide a flag as a workaround

I will see if I can figure out why your file triggers this issue. Thanks for letting me know about it.

On Wed, May 22, 2024 at 8:39 AM Kristine Larson @.***> wrote:

i know for sure that my code will not allow 30 observable types. there is no reason to have doppler (all the observations starting with D). i will try to change the code so that it gives you a better warning about this.

in the mean time, try this command.

smoosh zjgh2370.20o 15

but run it where the file exists. i.e. ls zjgh* should be where you run the command.

this will remove all the extra observables and only keep SNR. and the 15 decimates so current sampling, which is 15. and then try rinex2snr.

On Wed, May 22, 2024 at 8:32 AM zihelyu @.***> wrote:

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

Already uploaded to this link http://kh.dawnwl.com/zjgh2370.20o

— Reply to this email directly, view it on GitHub https://github.com/kristinemlarson/gnssrefl/issues/279#issuecomment-2125094749, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFSDAPCRWV5P3ZX4KC5HFNLZDS27BAVCNFSM6AAAAABIDX4QJCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCMRVGA4TINZUHE . You are receiving this because you commented.Message ID: @.***>

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

kristinemlarson commented 5 months ago

from what I read you are creating RINEX files with incorrect strings, i.e. they are using Latin instead of whatever python expects (utf-8). i think smoosh uses gfzrnx and maybe it is friendlier to different strings.

On Wed, May 22, 2024 at 8:51 AM Kristine Larson @.***> wrote:

I investigated a bit further and there is something about your file that does not like the line where I tested for RINEX version number:

rinext =np.loadtxt(r,usecols=0,dtype='str',max_rows=1)

it works for RINEX files from typical GNSS archives.

I am not sure why it does this - and I could provide a flag as a workaround - but I don't think that will help you much because my code does not allow RINEX 2.11 files with more than 25 observables. So you will need to fix them regardless. And since most of those observables are zero, it is best that you stop writing them out.

I will see if I can figure out why your file triggers this issue. Thanks for letting me know about it.

On Wed, May 22, 2024 at 8:39 AM Kristine Larson < @.***> wrote:

i know for sure that my code will not allow 30 observable types. there is no reason to have doppler (all the observations starting with D). i will try to change the code so that it gives you a better warning about this.

in the mean time, try this command.

smoosh zjgh2370.20o 15

but run it where the file exists. i.e. ls zjgh* should be where you run the command.

this will remove all the extra observables and only keep SNR. and the 15 decimates so current sampling, which is 15. and then try rinex2snr.

On Wed, May 22, 2024 at 8:32 AM zihelyu @.***> wrote:

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

Already uploaded to this link http://kh.dawnwl.com/zjgh2370.20o

— Reply to this email directly, view it on GitHub https://github.com/kristinemlarson/gnssrefl/issues/279#issuecomment-2125094749, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFSDAPCRWV5P3ZX4KC5HFNLZDS27BAVCNFSM6AAAAABIDX4QJCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCMRVGA4TINZUHE . You are receiving this because you commented.Message ID: @.***>

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

kristinemlarson commented 5 months ago

yeap, if i change the loadtxt command to say encoding = 'latin-1', it does not crash. i could put it in a try/except i guess....

On Wed, May 22, 2024 at 8:55 AM Kristine Larson @.***> wrote:

from what I read you are creating RINEX files with incorrect strings, i.e. they are using Latin instead of whatever python expects (utf-8). i think smoosh uses gfzrnx and maybe it is friendlier to different strings.

On Wed, May 22, 2024 at 8:51 AM Kristine Larson < @.***> wrote:

I investigated a bit further and there is something about your file that does not like the line where I tested for RINEX version number:

rinext =np.loadtxt(r,usecols=0,dtype='str',max_rows=1)

it works for RINEX files from typical GNSS archives.

I am not sure why it does this - and I could provide a flag as a workaround - but I don't think that will help you much because my code does not allow RINEX 2.11 files with more than 25 observables. So you will need to fix them regardless. And since most of those observables are zero, it is best that you stop writing them out.

I will see if I can figure out why your file triggers this issue. Thanks for letting me know about it.

On Wed, May 22, 2024 at 8:39 AM Kristine Larson < @.***> wrote:

i know for sure that my code will not allow 30 observable types. there is no reason to have doppler (all the observations starting with D). i will try to change the code so that it gives you a better warning about this.

in the mean time, try this command.

smoosh zjgh2370.20o 15

but run it where the file exists. i.e. ls zjgh* should be where you run the command.

this will remove all the extra observables and only keep SNR. and the 15 decimates so current sampling, which is 15. and then try rinex2snr.

On Wed, May 22, 2024 at 8:32 AM zihelyu @.***> wrote:

I don't know what the problem is. If you can post the rinex file somewhere public (but not on github and maybe gzipped), I will look at it.

Already uploaded to this link http://kh.dawnwl.com/zjgh2370.20o

— Reply to this email directly, view it on GitHub https://github.com/kristinemlarson/gnssrefl/issues/279#issuecomment-2125094749, or unsubscribe https://github.com/notifications/unsubscribe-auth/AFSDAPCRWV5P3ZX4KC5HFNLZDS27BAVCNFSM6AAAAABIDX4QJCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDCMRVGA4TINZUHE . You are receiving this because you commented.Message ID: @.***>

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

-- Kristine M. Larson @.*** https://www.kristinelarson.net http://www.kristinelarson.net https://gnss-reflections.org https://github.com/kristinemlarson

Are you waiting for your code to be perfect?

zihelyu commented 5 months ago

(utf-8

I am running it for my friend. I don’t know anything about GPS. Can you give me the running command of rinex2snr corresponding to this file? Thank you very much for answering my questions.

kristinemlarson commented 5 months ago

I have just pushed a new version of the gnssrefl code to github. It will take about 15 minutes to make a new docker. When that version is ready, please update your docker, and then test with your existing RINEX file that you were storing on your system with this command:

rinex2snr zjgh 2020 237 -orb gnss -strip T -nolook T

the extra boolean command strip will remove those extra observables I was previously talking about. and the command that was failing is now behind a try/except.

zihelyu commented 5 months ago

zjgh2370.20o

OK, I have successfully run and generated images, thank you very much for your help.

kristinemlarson commented 5 months ago

i will make a new version number for this later today. thanks again for helping me figure this out.