Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1702 1427 (N25W56 L=346)   801    141        3 (1659/1368-1738/1493) -90.0 W53-W58 N21-N29
      2 1127 1431 (N30W09 L=299)   556    140        2 (1095/1380-1159/1492)  50.2 W06-W12 N26-N33
      3  293 1414 (N30E61 L=229)   480    129        0 ( 259/1364- 340/1449)  40.6 E64-E56 N26-N32
      4  301 1335 (N25E56 L=233)   422     57        2 ( 283/1304- 319/1371)  78.7 E58-E54 N22-N27
      5 1211 1332 (N23W14 L=304)   432    106        0 (1155/1313-1261/1356) -23.2 W10-W17 N21-N24
      6 1751  855 (S10W51 L=341)   587     58        1 (1723/ 827-1775/ 882) -38.7 W48-W54 S12-S08
      7 1087  807 (S09W03 L=293)   338     82        1 (1057/ 772-1110/ 840) -53.1 W01-W04 S11-S07
      8  900  775 (S11E09 L=281)   277     69        1 ( 864/ 747- 932/ 799) -14.0 E10-E06 S12-S09
      9 1353  743 (S15W20 L=310)   219     68        0 (1339/ 708-1363/ 783)  76.0 W19-W20 S16-S12
     10 1704  589 (S28W53 L=343)  1713    275        4 (1674/ 462-1748/ 693) -58.6 W48-W61 S38-S20
     11 1624  603 (S26W43 L=333)   466    127        1 (1600/ 559-1654/ 653)  66.8 W41-W45 S29-S22
  

    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