Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1608 1530 (N17W52 L= 35)   334     64        1 (1579/1508-1639/1549)   0.0 W48-W55 N16-N18
      2  780 1356 (N27E07 L=335)   220     81        0 ( 746/1331- 819/1383) -59.0 E09-E05 N24-N29
      3 1842 1212 (S06W60 L= 42)  2972    357        8 (1794/1063-1896/1341)  47.3 W54-W67 S13-N01
      4 1748 1035 (S13W46 L= 29)  4121    877       10 (1588/ 650-1864/1307) -67.1 W40-W60 S40-N05
      5  717 1244 (N22E13 L=329)   195     91        0 ( 673/1235- 761/1254) -18.4 E16-E10 N20-N22
      6 1125 1210 (N11W10 L=352)   241     89        0 (1085/1192-1174/1229) -31.0 W07-W12 N09-N12
      7 1325  872 (S12W14 L=356)   634    163        3 (1282/ 832-1379/ 914)   9.5 W11-W17 S14-S09
      8 1030  535 (S25E11 L=331)  4209   1117        7 ( 727/ 178-1424/ 795)  18.2 E55-W18 S43-S18
      9  434  485 (S16E53 L=289)   404     75        1 ( 401/ 458- 457/ 510) -53.1 E55-E51 S17-S13
     10  676  328 (S31E44 L=299)   356     93        0 ( 630/ 299- 708/ 348) -59.0 E46-E41 S33-S28
     11  480  292 (S29E66 L=276)   573     44        2 ( 455/ 271- 508/ 311)   5.2 E72-E61 S30-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