Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  506 1618 (N23E48 L= 45)   639    112        2 ( 490/1569- 525/1676) -53.1 E52-E45 N19-N27
      2  622 1604 (N23E38 L= 55)   359    109        0 ( 599/1558- 647/1654) -82.9 E39-E35 N19-N27
      3  334 1606 (N21E67 L= 27)   722     59        2 ( 314/1571- 359/1639) -32.0 E71-E62 N18-N23
      4  952 1587 (N27E15 L= 79)   764    202        0 ( 854/1522-1018/1631) -21.0 E22-E09 N23-N29
      5  366 1518 (N15E56 L= 38)   770    118        0 ( 316/1494- 431/1544)  -8.1 E63-E49 N13-N16
      6  304 1452 (N10E59 L= 34)   399     83        0 ( 286/1413- 318/1493) -71.6 E60-E57 N07-N13
      7  328 1148 (S08E48 L= 45)  1230    223        1 ( 277/1071- 387/1253)  65.6 E52-E44 S12-S01
      8  328 1048 (S14E47 L= 46)   412     72        2 ( 302/1017- 351/1083)  68.2 E49-E44 S16-S11
      9 1834  947 (N05W57 L=151)  1121    232        2 (1811/ 829-1876/1035) -63.4 W53-W61 S01-N09
     10  481  976 (S17E34 L= 59)   304     41        0 ( 458/ 955- 508/1002) -33.7 E35-E31 S17-S15
     11  154  949 (S21E67 L= 27)   700     84        1 ( 135/ 898- 177/ 990) -71.6 E69-E63 S24-S18
     12  394  878 (S24E40 L= 53)   906    205        3 ( 349/ 779- 431/ 978)  80.5 E44-E37 S29-S17
     13  568  776 (S28E26 L= 67)  1355    479        3 ( 527/ 589- 630/ 977)  90.0 E29-E23 S39-S15
     14  260  745 (S33E56 L= 37)   398     29        1 ( 237/ 725- 287/ 763)   9.5 E58-E52 S34-S32
     15 1747  658 (S13W58 L=152)   475    130        0 (1715/ 609-1769/ 717) -81.9 W57-W60 S16-S09
     16  843  599 (S35E05 L= 89)   250     74        1 ( 828/ 557- 865/ 632)  78.7 E05-E03 S37-S32
     17 1349  404 (S38W39 L=132)   286     47        1 (1326/ 374-1368/ 431) -63.4 W37-W40 S39-S35
  

    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