Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1613 (N17E57 L=278)  1137    252        4 ( 405/1533- 475/1706) -54.0 E63-E53 N12-N23
      2  243 1456 (N04E71 L=263)  2368    316        7 ( 185/1410- 289/1525)  45.0 E82-E64 N00-N09
      3 1727 1116 (N20W43 L= 18)   339    111        0 (1707/1058-1744/1172) -50.2 W40-W45 N17-N23
      4  206  907 (S29E56 L=279)   593     79        1 ( 188/ 872- 229/ 946)  78.7 E57-E53 S31-S26
      5  528  884 (S23E27 L=308)   295     79        0 ( 491/ 851- 562/ 918) -18.4 E30-E24 S24-S22
      6 1826  815 (N06W61 L= 36)   617    112        0 (1792/ 758-1853/ 868)  90.0 W59-W62 N02-N09
      7 1581  616 (S12W47 L= 22)   309     88        2 (1542/ 584-1623/ 641)  51.3 W45-W49 S14-S09
      8 1603  553 (S15W51 L= 26)   711    206        0 (1504/ 515-1675/ 603)  39.3 W44-W55 S18-S09
      9 1712  497 (S14W66 L= 41)  1131    222        3 (1671/ 460-1743/ 538)  59.0 W61-W70 S16-S11
     10 1139  383 (S38W29 L=  4)   892    244        4 (1110/ 307-1170/ 488) -37.6 W21-W35 S43-S32
     11 1305  307 (S38W45 L= 20)  1024    230        2 (1260/ 284-1354/ 340)   0.0 W40-W49 S40-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