Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1328 1622 (N35W38 L=169)   607    147        2 (1251/1607-1394/1637) -26.6 W32-W42 N32-N37
      2  711 1434 (N38E12 L=119)   750    163        4 ( 656/1361- 760/1474) -56.3 E16-E08 N32-N42
      3  453 1235 (N30E35 L= 97)  1269    262        0 ( 382/1146- 518/1331) -54.5 E40-E30 N23-N37
      4  979 1240 (N20W02 L=133)   930    300        4 ( 847/1187-1142/1309) -38.7 E05-W10 N14-N27
      5 1263 1039 (N03W14 L=146)  1124    235        2 (1179/ 933-1329/1144) -69.0 W11-W16 S03-N10
      6 1815 1029 (S12W54 L=185)  1387    266        5 (1771/ 930-1845/1138)  71.6 W50-W57 S17-S05
      7  483 1016 (N17E35 L= 97)   372    127        0 ( 455/ 983- 535/1070) -56.3 E36-E31 N14-N20
      8  245  980 (N18E56 L= 76)   460     67        1 ( 227/ 949- 262/1020) -78.7 E57-E53 N15-N20
      9 1265  705 (S17W08 L=139)   874    174        3 (1196/ 631-1302/ 762) -58.0 W03-W10 S22-S13
     10  289  502 (S13E70 L= 62)   884    196        1 ( 270/ 473- 323/ 529)  29.7 E74-E64 S15-S11
     11  996  340 (S36E20 L=112)  1165    311        1 ( 904/ 266-1043/ 429)  29.4 E29-E12 S39-S30
     12  584  328 (S29E50 L= 82)   993    183        2 ( 527/ 295- 620/ 365)   6.3 E54-E45 S31-S26
  

    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