Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1189 1702 (N40W29 L=111)   633    118        0 (1151/1663-1233/1737)  24.0 W24-W33 N37-N41
      2  808 1649 (N45E02 L= 80)   283     37        0 ( 786/1628- 832/1670) -63.4 E03-W00 N42-N46
      3 1628 1579 (N22W55 L=137)   749    134        1 (1583/1543-1661/1605)   7.1 W50-W58 N19-N23
      4  906 1563 (N36W03 L= 85)   286     22        0 ( 880/1547- 931/1579) -33.7 W01-W05 N35-N37
      5  787 1342 (N24E09 L= 73)   254     65        0 ( 748/1325- 824/1358) -21.8 E11-E06 N23-N25
      6 1670 1036 (S10W41 L=123)  1940    424        0 (1569/ 866-1819/1162)  26.6 W31-W56 S18-S04
      7 1870 1078 (S12W60 L=142)   656    101        0 (1835/1033-1897/1112) -56.3 W56-W63 S14-S08
      8  788  289 (S40E36 L= 46)   417     85        0 ( 772/ 261- 806/ 320)  38.7 E38-E33 S42-S37
  

    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