Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1607 1525 (N27W45 L= 20)   473     79        2 (1571/1492-1632/1557) -36.9 W42-W46 N24-N28
      2 1773 1389 (N17W56 L= 31)   422     89        0 (1737/1360-1804/1420) -38.7 W53-W58 N15-N19
      3 1563 1384 (N17W38 L= 13)   240     71        2 (1531/1365-1596/1401) -11.3 W35-W40 N16-N18
      4 1636 1297 (N12W42 L= 17)   250     46        0 (1624/1270-1653/1328) -90.0 W40-W42 N10-N13
      5 1253 1155 (N04W14 L=350)  1554    388        4 (1121/ 972-1325/1272) -59.5 W06-W18 S06-N12
      6  826 1107 (N03E11 L=324)   547    213        3 ( 790/1053- 884/1164) -50.2 E13-E08 N00-N06
      7  335  914 (S05E46 L=289)   367     79        3 ( 312/ 888- 359/ 945) -45.0 E48-E44 S06-S03
      8 1777  817 (S17W52 L= 28)   638    121        1 (1740/ 752-1818/ 862)  45.0 W49-W55 S21-S14
      9  265  770 (S13E55 L=280)  1172    193        2 ( 223/ 713- 321/ 820)  20.0 E60-E48 S16-S10
     10  942  619 (S27E07 L=328)   289     85        0 ( 914/ 579- 962/ 663)  63.4 E09-E05 S29-S23
     11 1511  436 (S42W37 L= 12)  1844    366        5 (1355/ 378-1670/ 520)  15.6 W23-W49 S46-S36
  

    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