Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1131 1628 (N41W25 L=338)   293     52        1 (1112/1598-1147/1656)  53.1 W23-W27 N38-N42
      2 1754 1428 (N11W62 L= 15)  1479    202        4 (1702/1371-1803/1512) -90.0 W58-W68 N07-N17
      3 1183 1486 (N30W22 L=335)   296     57        1 (1162/1458-1217/1509)   0.0 W19-W24 N29-N32
      4 1490 1209 (N07W32 L=345)   845    231        0 (1469/1096-1509/1315)  72.9 W28-W35 N00-N13
      5  160  952 (N17E66 L=247)   844    192        2 ( 117/ 906- 204/ 998) -36.9 E70-E60 N13-N19
      6 1600  789 (S20W31 L=344)   741    157        2 (1549/ 724-1651/ 863)  31.0 W26-W36 S23-S17
      7  297  737 (N03E56 L=257)   838    164        1 ( 245/ 665- 340/ 815) -69.4 E59-E53 S01-N07
      8 1357  445 (S37W09 L=322)   321     49        2 (1335/ 423-1378/ 467) -63.4 W07-W10 S39-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