Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  590 1677 (N25E49 L=247)  1560    450        2 ( 532/1553- 629/1788) -35.0 E61-E40 N18-N32
      2  813 1737 (N35E37 L=259)   568    132        1 ( 769/1707- 871/1763)   0.0 E41-E31 N33-N36
      3  265 1487 (N07E71 L=226)   598    122        1 ( 253/1446- 282/1521) -45.0 E74-E67 N04-N09
      4  527 1413 (N08E41 L=255)   541     97        2 ( 487/1384- 561/1443)  -9.5 E44-E38 N05-N09
      5  230 1189 (S11E59 L=237)  1204    172        7 ( 205/1117- 264/1271)  90.0 E61-E55 S15-S05
      6  170 1122 (S16E64 L=232)  1155    135        4 ( 141/1066- 193/1168) -71.6 E68-E61 S19-S13
      7  635 1001 (S14E23 L=274)   465     90        2 ( 595/ 959- 670/1048)  90.0 E24-E20 S16-S10
      8  185  836 (S34E58 L=238)   570     63        1 ( 164/ 811- 208/ 859) -11.3 E60-E55 S35-S32
      9  606  761 (S29E20 L=276)   356     67        0 ( 580/ 727- 635/ 787)  59.0 E21-E18 S31-S26
     10  458  625 (S42E31 L=265)   654    137        1 ( 413/ 592- 505/ 665)  37.9 E35-E26 S45-S38
     11 1550  371 (S27W58 L=354)   480    107        0 (1517/ 355-1581/ 391)   9.5 W54-W61 S27-S25
  

    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