Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1237 1303 (N14W18 L=137)   730    232        1 (1119/1247-1353/1354) -39.8 W11-W23 N08-N18
      2  727 1281 (N21E14 L=106)   216     98        0 ( 691/1248- 776/1311) -51.3 E15-E11 N18-N23
      3 1145 1197 (N09W10 L=129)   205     57        0 (1129/1172-1173/1229) -76.0 W09-W11 N07-N11
      4 1459  993 (S08W25 L=145)   525    157        0 (1408/ 918-1493/1079)  48.8 W21-W28 S11-S03
      5 1760  995 (S13W48 L=168)  1259    289        5 (1680/ 880-1799/1073) -77.0 W42-W52 S20-S07
      6  336  983 (N10E45 L= 74)  1034    314        2 ( 221/ 914- 430/1023) -22.4 E55-E38 N04-N13
      7 1186  980 (S04W09 L=128)   539    149        1 (1105/ 944-1252/1017) -26.6 W03-W11 S07-S01
      8 1724  955 (S15W44 L=164)   321     91        0 (1702/ 920-1746/ 991)  63.4 W42-W46 S17-S13
      9  675  838 (S04E24 L= 96)   220     65        0 ( 650/ 808- 709/ 872) -63.4 E25-E22 S05-S01
     10 1615  769 (S25W34 L=154)  1235    385        1 (1445/ 707-1752/ 838)   4.6 W21-W46 S26-S22
     11 1686  674 (S32W40 L=160)  1646    307        3 (1564/ 587-1792/ 765)  20.2 W30-W49 S36-S27
  

    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