Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  684 1687 (N58E17 L=  2)  1969    232        6 ( 600/1617- 791/1729)   0.0 E28-E04 N52-N61
      2 1236 1646 (N41W32 L= 51)  2002    455        2 (1034/1615-1439/1670) -21.2 W16-W47 N35-N47
      3  259 1337 (N36E60 L=319)   632    116        0 ( 221/1289- 301/1367)  32.0 E65-E54 N33-N38
      4  480 1313 (N33E35 L=345)   347     67        2 ( 452/1290- 515/1334)  11.3 E36-E31 N31-N34
      5 1176 1232 (N16W13 L= 32)   761    205        2 (1103/1149-1222/1305) -74.7 W10-W14 N10-N21
      6 1811 1043 (S08W55 L= 75)   462     95        1 (1792/1004-1837/1084)  45.0 W53-W59 S10-S06
      7  475 1052 (N16E36 L=344)   232     52        0 ( 446/1038- 501/1070)  14.0 E37-E33 N15-N17
      8 1594  697 (S24W33 L= 52)  1217    172        3 (1559/ 618-1621/ 783) -84.8 W30-W35 S30-S18
      9 1299  691 (S18W11 L= 30)   723    171        2 (1249/ 615-1356/ 748)  45.0 W07-W15 S22-S15
  

    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