Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  625 1484 (N31E24 L=193)   443     90        0 ( 585/1460- 666/1508) -18.4 E26-E20 N28-N32
      2 1601 1486 (N22W44 L=261)   311     85        0 (1568/1469-1633/1508)  18.4 W41-W47 N21-N23
      3 1261 1375 (N18W18 L=234)   324    101        2 (1209/1350-1322/1399) -15.9 W14-W21 N15-N19
      4  474 1338 (N22E34 L=182)   443    130        1 ( 425/1290- 502/1381) -35.5 E38-E31 N19-N25
      5 1489 1168 (N03W30 L=246)   415     96        1 (1472/1113-1504/1219) -90.0 W28-W31 N00-N06
      6 1746 1116 (S01W49 L=265)  1004    210        2 (1702/1071-1812/1194) -41.2 W45-W55 S04-N03
      7  889 1158 (N08E07 L=209)   214     96        1 ( 848/1126- 928/1186)  36.9 E09-E04 N05-N09
      8  525 1081 (N06E30 L=186)   350     44        1 ( 502/1046- 548/1108)  71.6 E32-E28 N04-N07
      9 1672 1027 (S06W42 L=258)   320     38        1 (1647/1007-1694/1054) -56.3 W40-W43 S07-S04
     10 1763  960 (S11W49 L=266)   778    135        2 (1718/ 912-1800/1029)  45.0 W45-W53 S13-S06
     11 1599  904 (S13W36 L=252)   539    102        1 (1567/ 854-1630/ 941)  78.7 W33-W38 S15-S10
     12  298  791 (S09E52 L=164)  2394    436        3 ( 217/ 640- 392/ 894) -69.4 E59-E45 S18-S02
     13  682  848 (S09E22 L=194)   217     67        0 ( 667/ 810- 696/ 886)  78.7 E23-E21 S11-S06
     14 1552  722 (S24W33 L=249)   987    255        0 (1455/ 601-1648/ 836)  47.1 W26-W39 S31-S17
     15  404  620 (S20E47 L=169)   804    214        3 ( 375/ 533- 432/ 685)  45.0 E52-E42 S25-S16
     16  667  563 (S26E29 L=188)   402     34        1 ( 640/ 539- 694/ 588)  26.6 E30-E26 S27-S24
     17 1585  422 (S45W42 L=258)   858    103        1 (1555/ 385-1609/ 470) -74.1 W39-W45 S48-S41
     18 1498  399 (S46W34 L=250)   460    119        1 (1456/ 385-1544/ 413) -11.3 W29-W39 S47-S44
  

    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