Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1005 1761 (N38E24 L=271)   768    128        1 ( 960/1735-1062/1795)  -5.7 E28-E18 N36-N40
      2  395 1623 (N15E61 L=234)  2249    313        5 ( 354/1541- 446/1695) -65.6 E66-E55 N09-N20
      3  488 1453 (N06E44 L=250)   802    174        2 ( 428/1415- 556/1493)  -4.8 E50-E38 N04-N07
      4  917 1331 (N07E13 L=281)  1297    344        5 ( 881/1198- 949/1469) -59.0 E17-E08 N00-N15
      5 1710 1104 (N15W40 L=335)  1184    217        0 (1663/1004-1751/1220) -45.0 W34-W45 N10-N21
      6  493 1085 (S15E34 L=261)   410    165        0 ( 425/1065- 568/1109)   5.2 E39-E28 S15-S14
      7  247 1016 (S23E53 L=242)  1756    335        0 ( 161/ 938- 353/1080)  29.4 E60-E43 S28-S18
      8  733 1058 (S12E17 L=277)   216     74        0 ( 692/1036- 769/1078)  -9.5 E20-E14 S12-S11
      9 1753  549 (S12W65 L=360)  1660    203        1 (1677/ 506-1805/ 609)  48.4 W59-W70 S16-S07
  

    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