Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  737 1741 (N62W04 L=244)   596     97        1 ( 701/1702- 781/1778)   9.5 W00-W12 N58-N64
      2 1424 1688 (N32W53 L=292)  1258    192        3 (1337/1649-1525/1723) -45.0 W48-W59 N26-N37
      3 1100 1642 (N41W27 L=266)   346     62        0 (1083/1617-1128/1677)  78.7 W25-W28 N38-N43
      4 1141 1476 (N29W21 L=260)   209     59        0 (1112/1458-1175/1489) -11.3 W18-W23 N28-N30
      5  450 1196 (N31E34 L=205)   409    111        1 ( 393/1173- 502/1217)  -6.3 E38-E29 N29-N31
      6 1836  917 (S23W53 L=293)   790     67        1 (1816/ 874-1858/ 950)  53.1 W50-W56 S25-S20
      7 1798  802 (S29W48 L=287)   814    129        3 (1755/ 773-1854/ 842) -26.6 W43-W55 S32-S25
      8 1109  655 (S16E05 L=234)   325    121        0 (1073/ 636-1168/ 687) -38.7 E06-E01 S17-S13
      9 1713  645 (S37W38 L=278)   336     39        1 (1694/ 621-1728/ 668)  45.0 W36-W40 S38-S35
     10 1067  503 (S25E12 L=227)   748    273        0 ( 998/ 453-1125/ 569) -29.7 E16-E08 S27-S20
     11  683  484 (S16E36 L=203)   405    135        0 ( 621/ 463- 742/ 511) -35.5 E39-E32 S18-S13
     12  736  443 (S20E35 L=204)   487    132        2 ( 684/ 426- 813/ 466) -29.7 E37-E30 S22-S17
     13  641  386 (S22E44 L=195)  1387    189        2 ( 567/ 352- 729/ 432) -33.7 E49-E37 S25-S17
  

    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