Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1821 (N42E31 L=175)  1101    212        2 ( 933/1759-1022/1888) -33.7 E38-E23 N37-N47
      2  932 1496 (N17E18 L=188)   466     90        0 ( 884/1453- 977/1536)  66.8 E19-E16 N13-N20
      3  902 1327 (N07E14 L=191)   473    136        0 ( 871/1262- 922/1397) -82.9 E16-E13 N02-N10
      4 1069 1195 (N03E02 L=204)   714    235        1 ( 943/1164-1182/1216)  17.1 E08-W05 N00-N05
      5 1014  907 (S14W02 L=208)   883    225        1 ( 940/ 870-1113/ 944)  45.0 E01-W06 S17-S09
      6 1738  895 (N06W48 L=253)   399    102        0 (1704/ 875-1785/ 917)   7.1 W44-W52 N04-N06
      7  730  596 (S40E07 L=198)  1628    332        3 ( 576/ 571- 905/ 631)  14.0 E19-W05 S43-S35
      8 1349  551 (S24W34 L=239)  1308    378        0 (1181/ 522-1519/ 588)  23.5 W23-W46 S29-S19
  

    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