Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1115 1853 (N54E08 L=341)   537     67        0 (1081/1829-1142/1880)  56.3 E10-E06 N51-N58
      2  470 1709 (N33E58 L=291)   679     71        3 ( 449/1672- 486/1747) -36.9 E63-E53 N29-N35
      3 1208 1507 (N26W05 L=354)   670    205        0 (1167/1430-1238/1586) -70.0 W01-W07 N20-N31
      4 1022 1411 (N17E05 L=344)   326     91        0 ( 992/1370-1051/1452) -45.0 E07-E02 N15-N19
      5 1466 1359 (N21W24 L= 13)   230     43        1 (1444/1339-1490/1375)   0.0 W22-W25 N19-N21
      6 1808 1231 (N19W53 L= 42)  1245    266        2 (1767/1119-1836/1350) -77.9 W49-W55 N11-N25
      7  659  793 (S24E21 L=328)   554    198        0 ( 614/ 708- 685/ 883)  79.7 E23-E18 S29-S18
      8  474  834 (S23E34 L=314)   254     79        0 ( 438/ 811- 494/ 869)  38.7 E37-E32 S24-S20
      9 1717  652 (S17W55 L= 44)   367     30        0 (1694/ 629-1736/ 672) -11.3 W52-W57 S17-S15
     10 1115  174 (S63W37 L= 26)   615     44        1 (1095/ 141-1131/ 207) -30.3 W32-W45 S66-S59
  

    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