Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1493 1716 (N54W13 L=249)   533     71        1 (1466/1685-1520/1739) -59.0 W09-W15 N50-N56
      2  535 1646 (N17E51 L=186)   502    166        1 ( 504/1600- 558/1688) -29.1 E55-E46 N14-N19
      3  472 1626 (N15E55 L=182)   444     60        1 ( 445/1603- 505/1645)  45.0 E57-E52 N13-N16
      4  844 1613 (N22E28 L=208)   239     81        0 ( 795/1599- 883/1624)  36.9 E30-E25 N20-N23
      5  411 1524 (N07E55 L=182)  1557    216        3 ( 328/1483- 503/1578)  28.3 E62-E48 N03-N11
      6 1011 1509 (N20E14 L=222)   309     73        1 ( 981/1485-1060/1533)  36.9 E16-E11 N18-N21
      7 1329 1362 (N20W08 L=245)   315    107        0 (1296/1303-1350/1413) -45.0 W05-W10 N17-N22
      8  194 1242 (S13E64 L=172)  2698    395        8 ( 119/1158- 243/1351)  90.0 E74-E58 S18-S05
      9  725 1174 (S06E21 L=216)   470    132        2 ( 655/1147- 781/1192)   7.1 E24-E16 S07-S04
     10  488  832 (S31E29 L=207)  1009    204        2 ( 422/ 801- 574/ 897)  -4.8 E35-E22 S33-S28
     11  861  749 (S27E02 L=234)   399    112        0 ( 797/ 726- 913/ 766)   8.1 E06-W01 S28-S24
     12 1252  730 (S17W21 L=258)   825    275        0 (1120/ 691-1384/ 758)  16.4 W13-W30 S19-S14
     13 1606  619 (S12W48 L=285)   334     93        1 (1570/ 588-1645/ 660)  71.6 W46-W49 S14-S08
     14 1696  614 (S09W56 L=293)   432     88        0 (1674/ 573-1715/ 656) -68.2 W54-W57 S11-S06
     15 1342  560 (S23W34 L=270)  1695    618        0 (1060/ 474-1548/ 617)  14.0 W15-W51 S30-S16
     16 1679  481 (S16W63 L=299)   494     84        1 (1650/ 450-1715/ 506)  51.3 W60-W65 S18-S13
     17 1217  411 (S36W34 L=270)   737    155        2 (1132/ 390-1290/ 435)  20.0 W27-W38 S37-S32
     18 1659  383 (S21W68 L=305)   628    148        1 (1639/ 356-1673/ 415) -15.9 W64-W71 S22-S19
  

    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