Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  709 1783 (N68E14 L=191)   429     72        2 ( 678/1755- 740/1804) -15.3 E19-E08 N64-N70
      2 1697 1544 (N22W67 L=271)  1457    246        1 (1579/1512-1800/1577) -23.7 W55-W81 N16-N27
      3  845 1261 (N24E07 L=198)   901    268        5 ( 766/1167- 908/1373) -75.1 E11-E04 N17-N32
      4 1619 1286 (N10W44 L=249)   276    117        1 (1578/1247-1670/1306)   6.3 W40-W49 N08-N11
      5 1833  936 (S16W57 L=262)   751    175        6 (1800/ 869-1881/1005)  36.9 W52-W64 S20-S13
      6 1341  805 (S12W15 L=220)   298    117        0 (1320/ 748-1372/ 847)  51.3 W13-W17 S15-S10
      7  993  791 (S06E06 L=199)   272    149        0 ( 967/ 731-1044/ 830)  33.7 E08-E02 S09-S04
      8  760  737 (S06E21 L=184)   467    171        2 ( 691/ 699- 837/ 777) -40.6 E24-E17 S09-S02
      9 1614  721 (S24W34 L=239)   280     74        1 (1572/ 699-1640/ 750)  18.4 W30-W36 S25-S22
     10 1754  710 (S29W48 L=253)   531     97        0 (1719/ 671-1793/ 746)  20.6 W44-W52 S30-S26
     11 1197  450 (S33W00 L=205)   327     99        0 (1164/ 417-1237/ 496)  38.7 E02-W03 S35-S30
     12  956  419 (S30E17 L=188)   328    160        0 ( 918/ 354-1002/ 466)  35.0 E22-E12 S34-S27
  

    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