Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  898 1688 (N35E20 L=305)   357    125        0 ( 844/1662- 950/1726) -16.7 E25-E15 N33-N37
      2  795 1525 (N22E22 L=303)   724    206        1 ( 745/1435- 831/1609) -45.0 E27-E17 N16-N27
      3 1478 1469 (N28W24 L=349)   240     48        0 (1459/1440-1496/1496) -63.4 W22-W25 N26-N30
      4 1419 1314 (N17W21 L=346)   777    198        0 (1337/1246-1481/1389) -28.6 W14-W25 N13-N20
      5 1816 1282 (N23W54 L= 20)   534     48        1 (1798/1254-1838/1311) -56.3 W52-W56 N20-N24
      6  148 1047 (S12E67 L=258)   643     55        2 ( 130/1015- 170/1079) -90.0 E69-E63 S14-S10
      7  303  953 (S18E49 L=276)   388     88        0 ( 278/ 917- 347/ 973)  23.2 E51-E44 S19-S16
      8  147  867 (S23E69 L=256)   919     91        1 ( 118/ 831- 173/ 919)  26.6 E76-E64 S25-S20
      9  293  881 (S22E50 L=275)   300     53        1 ( 267/ 851- 320/ 905) -36.9 E52-E47 S23-S20
     10  978  704 (S26W02 L=327)   746    276        0 ( 871/ 616-1033/ 818)  63.4 E04-W04 S32-S18
     11 1483  426 (S35W47 L= 12)   599     88        0 (1439/ 393-1531/ 462)  56.3 W44-W49 S37-S31
     12  499  325 (S60E50 L=275)   706     64        1 ( 472/ 300- 535/ 348)  14.7 E60-E41 S63-S58
  

    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