Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  722 1765 (N30E45 L=140)  1899    348        0 ( 641/1665- 811/1825) -81.9 E51-E35 N21-N36
      2 1142 1727 (N39E14 L=171)   647    127        1 (1096/1670-1173/1779) -32.5 E20-E09 N35-N42
      3  711 1637 (N21E37 L=148)   368     86        2 ( 671/1606- 733/1662) -18.4 E40-E34 N19-N23
      4  333 1512 (N06E62 L=123)   936    120        1 ( 285/1475- 377/1563)  50.2 E66-E58 N03-N10
      5 1702 1398 (N35W36 L=221)   812    276        1 (1658/1265-1734/1502) -68.2 W32-W40 N26-N41
      6  237 1129 (S18E57 L=128)   339     85        1 ( 217/1084- 256/1166)  59.0 E59-E55 S20-S15
      7 1822 1067 (N18W52 L=237)   752    153        0 (1800/1000-1844/1143)  90.0 W49-W55 N14-N23
      8  436  971 (S24E36 L=149)   458    118        2 ( 409/ 937- 482/1013) -18.4 E38-E32 S25-S21
      9  656  909 (S22E19 L=166)   538    200        1 ( 602/ 818- 686/ 987) -48.8 E24-E16 S27-S18
     10  182  867 (S34E64 L=121)  1022    153        1 ( 146/ 804- 206/ 934)  33.7 E72-E59 S38-S30
     11 1416  838 (S07W27 L=212)   298     61        0 (1398/ 799-1433/ 875) -78.7 W25-W28 S09-S04
     12  310  812 (S36E47 L=138)   384    110        0 ( 285/ 766- 346/ 864) -39.8 E49-E43 S38-S32
     13 1585  657 (S11W45 L=230)  1002    248        1 (1482/ 595-1691/ 705)  47.5 W40-W51 S17-S05
     14 1719  646 (S08W57 L=242)   374     39        2 (1700/ 623-1741/ 670) -45.0 W54-W59 S08-S06
     15 1292  521 (S28W32 L=217)   475    136        0 (1229/ 496-1351/ 554)  39.8 W28-W35 S30-S24
     16 1498  392 (S28W55 L=240)  1732    292        1 (1363/ 357-1625/ 425)  15.9 W44-W72 S32-S24
     17  596  345 (S64E15 L=170)  1421    183        5 ( 525/ 304- 687/ 391)   4.8 E26-E02 S66-S60
     18 1184  345 (S42W37 L=222)  1133    175        4 (1110/ 318-1262/ 378)  33.7 W29-W42 S44-S37
  

    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