Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1517 1606 (N32W49 L= 32)   577    191        2 (1439/1580-1563/1644) -10.3 W42-W54 N30-N33
      2 1258 1489 (N31W25 L=  8)   795    237        1 (1164/1438-1367/1535)  -3.6 W18-W34 N27-N32
      3 1746 1419 (N15W61 L= 43)  1438    197        2 (1683/1390-1809/1447) -26.6 W54-W67 N12-N18
      4  687 1102 (N16E20 L=323)   848    292        2 ( 546/1044- 799/1189) -34.5 E29-E13 N11-N22
      5 1328  951 (S03W17 L=360)  2145    440        4 (1271/ 741-1372/1160)  73.7 W10-W20 S14-N10
      6 1597  959 (S07W35 L= 18)   658    210        3 (1540/ 901-1677/1007) -45.0 W31-W40 S12-S03
      7 1506  667 (S24W26 L=  9)   984    179        0 (1474/ 595-1530/ 767)  65.6 W23-W28 S28-S17
      8 1112  385 (S37E06 L=337)  1438    316        3 (1022/ 304-1182/ 501)  38.7 E14-W01 S42-S28
  

    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