Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  910 1390 (N30E02 L=  5)  1276    326        2 ( 771/1339-1061/1433)  -5.7 E11-W09 N27-N32
      2 1550 1297 (N13W39 L= 45)   654    149        1 (1477/1258-1605/1323) -32.0 W33-W42 N10-N15
      3 1843 1053 (S07W59 L= 66)   943    239        0 (1807/ 940-1872/1155)  47.5 W53-W65 S13-S01
      4  869 1113 (N14E08 L=358)   188     68        0 ( 830/1095- 896/1126) -36.9 E10-E06 N12-N15
      5  361 1056 (N16E45 L=322)   598    196        1 ( 303/1019- 417/1120) -37.9 E49-E40 N13-N20
      6  207 1069 (N18E60 L=307)   401    130        0 ( 177/1024- 235/1100) -33.7 E63-E56 N15-N19
      7  954  864 (S02E07 L=360)   951    297        3 ( 911/ 739- 990/1028)  83.3 E08-E04 S09-N08
      8 1248  859 (S07W11 L= 17)   451    163        2 (1179/ 820-1334/ 890) -36.9 W07-W16 S10-S04
      9  944  711 (S11E10 L=357)   202     26        0 ( 922/ 693- 964/ 730) -45.0 E10-E08 S12-S10
     10 1135  546 (S25E00 L=  6)   393    100        0 (1112/ 502-1165/ 596)  51.3 E02-W02 S28-S22
     11  830  318 (S37E28 L=338)  1209    287        3 ( 766/ 255- 874/ 422)  49.8 E36-E24 S42-S29
  

    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