Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  913 1579 (N44W06 L=350)   767    168        2 ( 861/1513- 974/1644) -84.8 W03-W10 N38-N49
      2 1063 1554 (N38W16 L=360)   814    162        1 ( 993/1523-1144/1589) -36.9 W11-W21 N34-N41
      3 1389 1446 (N23W34 L= 17)   422    162        0 (1366/1379-1436/1487)  35.5 W31-W38 N18-N25
      4 1549 1319 (N12W40 L= 24)   383     89        1 (1509/1285-1586/1351)   8.1 W36-W43 N11-N13
      5  358 1308 (N36E45 L=299)  1063    205        2 ( 261/1263- 443/1339) -13.4 E57-E36 N32-N38
      6  582 1217 (N27E25 L=319)   692    156        3 ( 493/1177- 649/1244) -20.0 E31-E20 N23-N29
      7  521 1038 (N17E32 L=312)   290    147        0 ( 497/ 985- 568/1068)  38.7 E33-E27 N13-N18
      8  331  927 (N13E48 L=295)   452    129        2 ( 309/ 882- 362/ 983) -81.9 E50-E46 N09-N16
      9 1638  859 (S16W36 L= 19)   691    191        4 (1583/ 793-1690/ 948)  31.0 W30-W41 S19-S12
     10  960  672 (S12E11 L=333)   244    123        0 ( 913/ 639-1007/ 708)  15.9 E14-E07 S13-S11
  

    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