Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  497 1664 (N58E46 L= 28)   632     72        0 ( 470/1637- 530/1690)  68.2 E52-E41 N55-N60
      2 1703 1458 (N19W60 L=134)   544    131        0 (1644/1421-1746/1487) -32.0 W54-W62 N15-N21
      3  740 1436 (N36E13 L= 61)   453    125        1 ( 680/1419- 807/1447)  -5.7 E18-E08 N35-N37
      4  761 1185 (N20E14 L= 60)   354     99        1 ( 720/1161- 804/1210) -21.8 E16-E11 N17-N21
      5  380 1071 (N17E43 L= 31)   538    122        0 ( 348/1025- 411/1125) -56.3 E46-E40 N14-N20
      6 1350  594 (S26W14 L= 88)  1888    654        4 (1097/ 429-1547/ 785)  23.2 E05-W30 S33-S16
      7 1196  616 (S21W04 L= 78)   847    337        1 (1091/ 530-1306/ 743)  33.7 E03-W12 S25-S15
      8  699  585 (S16E28 L= 46)   491    128        0 ( 675/ 524- 731/ 637)  49.4 E30-E24 S19-S12
      9 1661  543 (S37W42 L=116)   478     85        0 (1626/ 522-1698/ 579)  15.9 W38-W45 S38-S35
     10 1434  454 (S38W20 L= 95)   620    171        0 (1368/ 404-1504/ 509)  24.4 W14-W26 S41-S35
     11  800  444 (S27E26 L= 48)   618    106        2 ( 779/ 400- 830/ 493)  81.9 E27-E23 S30-S23
  

    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