Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1398 1380 (N27W16 L=332)   655    133        1 (1331/1337-1442/1412)   0.0 W10-W19 N25-N28
      2  826 1337 (N12E19 L=298)   257     85        1 ( 779/1313- 856/1362)  11.3 E21-E16 N10-N12
      3 1783 1119 (N21W51 L=  7)   429     90        0 (1764/1083-1814/1163) -26.6 W48-W54 N18-N23
      4 1668  875 (N04W44 L=  1)   283     38        0 (1640/ 857-1696/ 894)  26.6 W42-W46 N02-N05
      5  619  829 (S22E21 L=295)   399     76        0 ( 568/ 814- 664/ 846)  26.6 E24-E18 S23-S20
      6  980  619 (S27W07 L=323)   220     76        0 ( 943/ 604-1019/ 642)  -9.5 W03-W09 S27-S26
      7 1598  481 (S21W56 L= 13)   965    190        0 (1513/ 441-1679/ 517)  42.0 W50-W62 S25-S16
      8 1663  464 (S20W64 L= 21)   532     46        0 (1639/ 440-1692/ 489)  45.0 W62-W67 S21-S18
  

    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