Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1390 1439 (N18W34 L= 53)   684    186        0 (1309/1415-1488/1460) -31.0 W29-W39 N14-N20
      2 1003 1385 (N26W09 L= 28)   249     98        1 ( 979/1333-1018/1434)  45.0 W06-W11 N24-N29
      3 1743 1356 (N02W56 L= 76)   927    249        2 (1690/1298-1811/1408)   7.1 W50-W66 N00-N03
      4 1513 1050 (S06W28 L= 48)   895    345        2 (1425/ 979-1652/1164) -37.6 W24-W38 S12-N01
      5 1794  982 (S19W48 L= 67)   507    145        0 (1764/ 936-1841/1031)  12.5 W44-W53 S20-S17
      6 1620  826 (S22W31 L= 50)  1054    322        3 (1585/ 743-1661/ 969)  56.3 W27-W36 S28-S14
      7 1875  857 (S30W54 L= 73)   934    264        3 (1828/ 822-1930/ 903)  11.3 W49-W61 S31-S26
      8  669  806 (N03E26 L=353)   345    102        0 ( 639/ 756- 696/ 855)  45.0 E29-E23 N00-N05
      9  559  750 (N02E34 L=345)   226     60        2 ( 527/ 734- 597/ 773) -53.1 E36-E32 N00-N04
     10  385  447 (S12E62 L=317)   425    110        2 ( 369/ 404- 403/ 482)  63.4 E63-E59 S15-S09
  

    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