Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1080 1390 (N28W07 L=156)   223     28        1 (1058/1370-1103/1406)  18.4 W05-W08 N26-N28
      2 1663 1331 (N17W49 L=198)   441    142        0 (1599/1312-1721/1368) -21.8 W44-W54 N15-N20
      3  872  938 (N01E10 L=139)   904    388        1 ( 820/ 800- 923/1089)  86.6 E13-E06 S07-N10
      4 1720 1055 (S00W48 L=197)   323     44        2 (1694/1027-1743/1076) -63.4 W45-W50 S01-N01
      5 1188  942 (S01W09 L=159)   252     95        0 (1146/ 909-1221/ 964) -36.9 W06-W11 S03-N00
      6 1861  857 (S15W64 L=213)   583    103        0 (1814/ 837-1907/ 883)   0.0 W57-W71 S15-S13
      7 1577  699 (S21W35 L=185)  1675    287        4 (1471/ 612-1660/ 797)  37.6 W27-W42 S25-S14
      8  915  716 (S13E09 L=140)   292    131        0 ( 858/ 694- 950/ 754)   9.5 E12-E06 S14-S10
      9 1379  488 (S33W21 L=170)   454     75        0 (1345/ 463-1409/ 520)  31.0 W18-W23 S35-S31
     10  724  341 (S39E32 L=117)   483    107        0 ( 698/ 295- 756/ 386)  37.9 E36-E27 S42-S35
  

    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