Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1134 1833 (N49E21 L=153)   598     87        1 (1088/1811-1174/1855)   0.0 E25-E16 N47-N50
      2  608 1685 (N23E47 L=126)  1223    207        1 ( 553/1595- 645/1779) -65.0 E53-E44 N15-N30
      3 1327 1362 (N20W09 L=183)   245     86        0 (1301/1329-1360/1396) -21.8 W06-W11 N18-N21
      4  911 1337 (N08E14 L=160)   261     50        0 ( 886/1319- 933/1359)  71.6 E15-E12 N06-N09
      5 1258 1209 (N09W09 L=182)   482    165        1 (1233/1139-1284/1297) -53.1 W05-W11 N05-N13
      6 1750 1012 (N12W46 L=220)   433     83        2 (1732/ 970-1765/1052) -45.0 W43-W47 N09-N14
      7  180  934 (S30E63 L=111)  1212    161        1 ( 113/ 903- 260/ 958)   9.5 E76-E52 S31-S27
      8  492  761 (S35E29 L=144)   648    221        0 ( 428/ 727- 554/ 822)  48.8 E35-E24 S37-S29
      9  238  793 (S38E57 L=117)   598    115        3 ( 191/ 767- 292/ 821)  18.4 E64-E49 S40-S35
     10  702  627 (S38E10 L=163)   212     63        0 ( 673/ 604- 732/ 651)  -9.5 E13-E07 S39-S37
     11  780  394 (S52W03 L=176)   518    104        1 ( 732/ 362- 826/ 437) -11.3 E02-W08 S54-S50
  

    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