Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1055 1740 (N42E08 L=205)   287     30        0 (1031/1725-1080/1756)  26.6 E09-E05 N40-N43
      2 1124 1420 (N19W02 L=214)   230     45        0 (1093/1402-1152/1438)  33.7 W00-W03 N17-N20
      3  701 1252 (N05E22 L=190)   206     62        0 ( 689/1215- 712/1290) -78.7 E22-E21 N02-N07
      4 1922 1095 (N12W67 L=279)   795    187        0 (1898/1060-1941/1129)  59.0 W63-W69 N09-N14
      5  985 1100 (S02E03 L=209)   298    105        0 ( 931/1069-1044/1122)  -9.5 E06-W00 S02-N00
      6 1897  946 (N02W65 L=278)   705    145        0 (1875/ 906-1928/ 994) -24.0 W61-W71 N00-N04
      7  701  901 (S16E19 L=194)   225    108        0 ( 681/ 848- 732/ 965)  76.0 E19-E17 S19-S11
      8 1081  882 (S13W05 L=217)   389    129        0 (1011/ 867-1148/ 899)  12.5 W00-W09 S14-S11
      9  174  768 (S27E62 L=150)  2600    371        4 ( 146/ 660- 216/ 860)  59.7 E66-E56 S33-S21
     10  490  811 (S23E33 L=179)   344    116        0 ( 473/ 745- 509/ 860)  81.9 E34-E31 S27-S20
     11 1172  832 (S15W11 L=224)   210     68        0 (1129/ 816-1207/ 852)  21.8 W08-W13 S16-S13
     12 1609  758 (S14W42 L=254)  1616    296        2 (1480/ 657-1733/ 851)  47.0 W33-W49 S21-S06
     13  815  714 (S26E10 L=202)   308     86        0 ( 769/ 670- 857/ 754)  50.2 E13-E07 S29-S23
     14 1292  714 (S21W21 L=233)   462    160        1 (1215/ 696-1374/ 726)  10.3 W15-W26 S22-S19
     15 1039  666 (S27W05 L=218)   305     84        0 ( 991/ 635-1087/ 689)  38.7 W02-W08 S29-S24
     16 1799  637 (S17W63 L=275)   552     73        2 (1779/ 604-1815/ 668)  90.0 W60-W65 S19-S15
     17 1663  517 (S27W54 L=267)   418     48        1 (1643/ 492-1685/ 539)  71.6 W52-W56 S28-S25
     18 1122  313 (S51W21 L=233)   305     30        0 (1100/ 294-1148/ 330)  11.3 W18-W23 S52-S48
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y