Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1508 (N53E48 L=  0)   560     74        0 ( 336/1484- 413/1528)   0.0 E53-E41 N51-N54
      2 1239 1477 (N27W26 L= 75)   271     90        1 (1212/1456-1275/1498) -11.3 W23-W28 N25-N28
      3 1033 1246 (N20W06 L= 54)   333    201        0 ( 956/1187-1090/1284) -49.4 W01-W08 N14-N22
      4  802 1079 (N15E12 L= 37)   594    129        3 ( 744/1031- 852/1134) -50.2 E14-E08 N12-N19
      5 1801  882 (S23W50 L= 98)  3234    646        2 (1606/ 679-1930/1117)  18.0 W29-W73 S30-S13
      6  512  962 (N14E33 L= 15)   267     98        2 ( 502/ 920- 540/1003)  63.4 E34-E30 N12-N16
      7  800  848 (N02E17 L= 31)  1032    379        2 ( 737/ 690- 867/1003)  49.1 E23-E09 S06-N09
      8 1644  849 (S19W35 L= 83)  1116    237        1 (1587/ 734-1704/ 936)  39.3 W28-W41 S24-S15
      9  977  550 (S19E15 L= 34)   280     75        0 ( 939/ 526-1013/ 570) -36.9 E16-E12 S21-S17
     10 1176  348 (S39E09 L= 39)  1387    233        2 (1087/ 270-1256/ 438)  22.8 E19-W00 S43-S34
     11  896  385 (S29E26 L= 22)   397    116        0 ( 845/ 371- 951/ 408) -33.7 E28-E22 S30-S25
  

    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