Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1251 1624 (N35E02 L= 34)   258     63        0 (1228/1591-1280/1657) -36.9 E04-W00 N33-N36
      2  190 1149 (S17E61 L=335)  1018    186        2 ( 165/1064- 213/1250) -84.8 E64-E59 S21-S10
      3  415 1178 (S12E41 L=355)   763    201        2 ( 381/1101- 449/1250) -58.0 E44-E37 S15-S07
      4  490 1171 (S11E36 L=  0)   353    120        0 ( 459/1135- 557/1197)  50.2 E37-E31 S13-S07
      5 1888 1000 (N16W60 L= 96)  2201    375        3 (1831/ 923-1929/1065) -53.1 W52-W66 N12-N20
      6  853  939 (S16E08 L= 28)   387    150        0 ( 811/ 860- 885/1008) -40.6 E11-E03 S19-S12
      7 1499  889 (S02W31 L= 67)  1054    225        2 (1408/ 836-1597/ 967)   0.0 W23-W38 S04-N00
      8 1080  858 (S15W07 L= 43)   182     70        0 (1042/ 828-1103/ 880)  18.4 W04-W08 S16-S14
      9 1018  824 (S18W05 L= 41)   196     34        0 (1006/ 796-1030/ 855) -56.3 W03-W05 S20-S16
     10  522  792 (S32E27 L=  9)   912    292        0 ( 396/ 767- 652/ 819)  10.8 E37-E16 S35-S29
     11 1859  778 (N04W65 L=101)   591     73        0 (1840/ 738-1878/ 816)  90.0 W63-W67 N00-N06
     12 1675  723 (S06W49 L= 85)   640    176        0 (1617/ 665-1713/ 761) -11.3 W43-W53 S08-S04
     13 1358  634 (S20W31 L= 67)   294     95        0 (1316/ 611-1410/ 652)   8.1 W27-W34 S20-S18
     14 1572  601 (S15W46 L= 82)   842    209        0 (1477/ 577-1662/ 643)  30.3 W39-W51 S18-S10
     15 1735  514 (S14W66 L=102)  4100    470        9 (1673/ 436-1797/ 618) -51.3 W58-W75 S19-S08
     16 1636  491 (S19W57 L= 93)   409     75        1 (1616/ 462-1661/ 521) -21.8 W54-W59 S20-S17
     17 1256  483 (S31W31 L= 67)   777    149        0 (1183/ 458-1341/ 513)  41.2 W27-W35 S34-S27
     18 1522  440 (S25W51 L= 87)   603    151        1 (1489/ 386-1543/ 509) -63.4 W48-W54 S28-S20
  

    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