Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  971 1631 (N35E19 L=179)  3039    595        6 ( 713/1513-1204/1775)  -6.5 E44-W00 N30-N39
      2  389 1593 (N20E61 L=137)  1030    315        1 ( 371/1506- 409/1680) -47.7 E67-E56 N14-N25
      3 1277 1500 (N33W06 L=204)   253     60        0 (1253/1475-1308/1516)  53.1 W04-W08 N30-N34
      4  608 1455 (N15E36 L=162)   409     70        0 ( 584/1419- 644/1500)  78.7 E37-E34 N12-N18
      5  369 1277 (N00E48 L=149)   558    109        2 ( 329/1244- 409/1312)  45.0 E51-E45 S02-N02
      6  726 1050 (S06E18 L=180)   205     57        0 ( 710/1017- 751/1079)  90.0 E19-E17 S08-S04
      7 1867  965 (N13W62 L=260)  1447    261        2 (1810/ 859-1902/1037) -39.3 W54-W68 N06-N16
      8  747  965 (S11E15 L=182)   235     56        0 ( 720/ 938- 778/ 997) -21.8 E17-E12 S11-S09
      9  911  916 (S10E04 L=193)   376     84        1 ( 864/ 891- 961/ 946)  38.7 E07-E02 S11-S07
     10 1663  683 (S08W50 L=248)   696    169        0 (1648/ 619-1688/ 752) -49.4 W47-W54 S11-S04
     11 1319  628 (S19W27 L=225)   249     33        1 (1292/ 612-1346/ 644)   0.0 W25-W28 S20-S17
     12 1000  573 (S29W09 L=206)   332     65        0 ( 962/ 554-1039/ 591)  11.3 W06-W11 S30-S28
  

    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