Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  907 1894 (N52E37 L=292)   911     98        2 ( 867/1871- 955/1916)  -8.1 E44-E30 N48-N54
      2  662 1797 (N38E48 L=281)   823    100        1 ( 639/1752- 688/1850) -39.3 E55-E42 N33-N42
      3  808 1726 (N34E31 L=299)   281     67        1 ( 790/1689- 829/1763) -45.0 E33-E28 N31-N36
      4 1245 1442 (N22W06 L=336)   326     82        2 (1213/1412-1280/1467) -21.8 W03-W08 N20-N23
      5 1046 1375 (N14E05 L=325)   194     68        0 (1014/1355-1078/1394)  45.0 E06-E03 N12-N16
      6 1865 1082 (N14W58 L= 28)   598     88        1 (1833/1056-1890/1125)  45.0 W54-W61 N12-N17
      7  457  960 (S19E35 L=294)   642    178        0 ( 412/ 897- 515/1006)  48.8 E38-E31 S23-S15
      8 1700  642 (S14W54 L= 24)   454     66        1 (1681/ 604-1722/ 674) -63.4 W52-W56 S16-S12
      9  342  581 (S45E51 L=279)   727     77        1 ( 317/ 548- 368/ 623) -90.0 E54-E46 S46-S41
     10 1680  522 (S21W60 L= 29)  1974    264        8 (1610/ 440-1768/ 598)  69.0 W53-W65 S27-S14
     11 1131  533 (S34W18 L=347)   471    164        0 (1048/ 516-1207/ 552)  26.6 W12-W22 S35-S30
     12  991  257 (S57W21 L=351)   887    210        1 ( 883/ 221-1075/ 303)   0.0 W07-W31 S58-S54
  

    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