Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1726 1488 (N43W41 L=142)   984    159        2 (1646/1444-1800/1525)  10.8 W31-W53 N40-N46
      2  366 1251 (S08E48 L= 53)   504    124        0 ( 320/1207- 414/1298)  49.4 E50-E44 S11-S04
      3 1789 1065 (N16W49 L=150)   726    155        2 (1757/ 999-1820/1147) -48.8 W44-W52 N12-N20
      4 1867 1029 (N17W58 L=159)  1207    241        0 (1843/ 939-1921/1137) -28.1 W54-W69 N12-N22
      5  511  956 (S23E30 L= 71)   395    159        0 ( 440/ 934- 574/ 983)  11.3 E35-E25 S24-S21
      6  324  745 (S40E47 L= 54)   748    137        1 ( 282/ 704- 358/ 810)  41.6 E52-E42 S42-S34
      7 1633  585 (S14W53 L=154)   840    247        0 (1546/ 499-1702/ 680)  81.9 W50-W54 S21-S07
      8 1197  565 (S28W24 L=125)   339     83        1 (1167/ 533-1234/ 613)  80.5 W22-W25 S30-S24
      9 1137  513 (S33W23 L=124)   346     68        1 (1107/ 493-1172/ 537)   0.0 W19-W25 S34-S31
  

    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