Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1266 1423 (N21W26 L= 99)   518    127        2 (1234/1368-1297/1487)  35.5 W22-W30 N18-N24
      2 1836 1252 (S06W61 L=134)  2876    390        2 (1778/1133-1917/1374) -76.0 W57-W66 S15-N01
      3 1071 1200 (N15W07 L= 80)   420    170        1 ( 987/1181-1158/1231) -16.7 W02-W12 N12-N17
      4  689 1060 (N17E18 L= 55)   852    292        2 ( 571/1015- 828/1101) -33.7 E25-E10 N11-N21
      5 1698 1037 (S13W41 L=114)   255     64        1 (1684/1012-1713/1081)  45.0 W39-W42 S13-S10
      6 1599  860 (S19W30 L=102)  1623    290        7 (1552/ 688-1630/1080)  56.3 W22-W36 S28-S07
      7  181  875 (N17E63 L= 10)   587    111        0 ( 167/ 831- 193/ 926)  45.0 E66-E60 N14-N19
      8  647  709 (S02E30 L= 43)   225     49        0 ( 633/ 680- 662/ 740)  45.0 E31-E28 S03-N00
      9 1343  572 (S29W06 L= 79)  1219    281        0 (1226/ 506-1418/ 681)  17.1 E01-W13 S32-S24
     10  808  618 (S11E23 L= 49)   347    150        0 ( 739/ 588- 877/ 656)   0.0 E28-E18 S11-S09
     11  888  509 (S19E23 L= 50)   593    198        0 ( 867/ 435- 918/ 605)  45.0 E26-E17 S23-S14
  

    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