Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1066 1632 (N39W23 L=161)   932    203        0 (1010/1547-1118/1697)  21.8 W15-W30 N34-N41
      2 1480 1600 (N24W48 L=185)  1092    161        2 (1414/1556-1552/1630) -49.4 W42-W51 N19-N26
      3  526 1247 (N31E25 L=112)   920    346        0 ( 387/1218- 680/1278) -22.2 E36-E14 N26-N35
      4 1532 1163 (S01W33 L=170)  1037    269        4 (1426/1123-1626/1200) -21.0 W25-W38 S05-N01
      5 1493 1091 (S04W28 L=166)   267     92        1 (1455/1058-1531/1109)   0.0 W25-W30 S04-S02
      6  796 1073 (N14E12 L=126)   184     62        1 ( 771/1057- 828/1101) -45.0 E12-E09 N12-N15
      7 1697  986 (S16W40 L=177)   641    119        1 (1651/ 949-1735/1024)  26.6 W36-W43 S17-S14
      8 1426  836 (S16W18 L=155)   618    161        2 (1405/ 764-1462/ 930)  78.7 W16-W20 S20-S10
      9 1099  689 (S15E05 L=133)   811    212        0 ( 972/ 648-1188/ 718) -39.3 E11-W00 S19-S10
     10 1195  565 (S25E03 L=134)   287     97        0 (1181/ 512-1209/ 600)  59.0 E04-E01 S28-S23
     11 1040  218 (S46E29 L=108)  1038    140        1 (1012/ 173-1065/ 264)  35.0 E34-E24 S49-S41
  

    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