Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1224 1834 (N61W40 L=212)   532     71        1 (1209/1795-1239/1866)  29.7 W34-W48 N56-N64
      2 1828 1324 (N15W67 L=239)   486     44        1 (1799/1307-1855/1339) -14.0 W63-W71 N13-N15
      3  555 1038 (N09E30 L=142)   324     81        0 ( 527/1014- 594/1070) -36.9 E32-E27 N06-N10
      4 1817  986 (S05W57 L=229)  1340    143        2 (1751/ 944-1882/1019) -14.0 W50-W66 S07-S02
      5  557  964 (N04E30 L=141)   275     57        1 ( 542/ 933- 574/ 995)  56.3 E31-E29 N02-N05
      6 1304  686 (S18W16 L=187)   602    137        1 (1244/ 647-1360/ 723)  26.6 W11-W19 S20-S16
      7  521  686 (S13E36 L=135)   802    230        2 ( 419/ 662- 601/ 722)   7.6 E44-E29 S14-S11
      8 1657  614 (S27W44 L=216)   586    108        0 (1617/ 575-1698/ 662)  39.8 W41-W47 S29-S24
      9 1106  472 (S32W01 L=173)  1197    311        1 (1032/ 373-1207/ 583)  49.4 E05-W08 S38-S24
  

    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