Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1439 1705 (N31W53 L=247)   979    215        1 (1374/1683-1517/1723) -17.1 W46-W59 N28-N34
      2 1014 1426 (N29W11 L=206)   320     82        0 ( 966/1399-1056/1454)   0.0 W07-W14 N28-N30
      3 1484 1406 (N14W39 L=233)  2692    630        4 (1213/1364-1708/1433) -21.8 W21-W56 N06-N21
      4  714 1234 (N27E13 L=181)   694    157        1 ( 674/1160- 743/1310)  84.3 E14-E10 N21-N31
      5 1857 1194 (S10W62 L=256)  2510    365        9 (1801/1094-1915/1260)  38.7 W55-W72 S15-S05
      6 1121 1135 (N10W08 L=203)   374    133        1 (1108/1063-1142/1207)  54.5 W06-W11 N06-N13
      7 1810 1107 (S12W53 L=248)   671    104        3 (1786/1071-1837/1152) -78.7 W50-W56 S14-S09
      8 1753 1071 (S12W47 L=241)  1522    339        4 (1711/ 964-1818/1146) -77.0 W42-W50 S20-S07
      9 1672 1043 (S11W40 L=234)   620    102        3 (1649/ 990-1695/1096)  81.9 W36-W41 S14-S07
     10 1164 1007 (N02W07 L=202)   471     93        0 (1111/ 984-1215/1034) -38.7 W04-W09 N00-N04
     11  316  942 (N19E48 L=146)   420     84        0 ( 301/ 892- 330/ 989)  71.6 E50-E46 N15-N21
     12 1518  791 (S20W23 L=217)  2055    424        4 (1404/ 625-1568/ 940) -87.5 W18-W25 S31-S08
     13  688  488 (S16E36 L=159)   721    258        1 ( 645/ 392- 731/ 580)  45.0 E42-E30 S20-S09
     14  460  544 (S07E50 L=145)   581    124        0 ( 417/ 512- 523/ 577) -39.8 E52-E45 S10-S04
     15  278  526 (S06E71 L=124)   807    181        2 ( 244/ 501- 324/ 566) -26.6 E75-E64 S07-S02
     16  921  341 (S31E28 L=166)  1547    322        3 ( 794/ 286-1007/ 453)  15.9 E39-E18 S34-S25
     17  552  419 (S17E49 L=146)   369     90        2 ( 507/ 401- 586/ 434) -36.9 E51-E46 S18-S15
     18  701  223 (S35E51 L=143)   802    184        2 ( 680/ 186- 730/ 265)  33.7 E56-E47 S37-S31
  

    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