Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  539 1565 (N32E44 L=179)  1021    146        1 ( 475/1527- 600/1607)  -8.1 E51-E37 N29-N34
      2  296 1525 (N26E70 L=153)  1109     75        2 ( 269/1488- 332/1561)  -7.8 E84-E62 N24-N28
      3  203  822 (S18E63 L=160)  1232    165        1 ( 181/ 773- 233/ 881)  49.4 E67-E58 S21-S14
      4  434  815 (S16E39 L=185)   326     80        1 ( 410/ 788- 467/ 847) -56.3 E40-E35 S17-S14
      5  811  788 (S15E12 L=212)   538    197        1 ( 722/ 751- 895/ 811)  -9.5 E17-E05 S16-S12
      6 1336  666 (S18W24 L=247)   326    104        0 (1315/ 620-1364/ 710)  90.0 W22-W24 S21-S15
      7 1461  660 (S18W33 L=256)   266     63        1 (1426/ 642-1488/ 681) -11.3 W30-W35 S18-S16
      8 1424  537 (S26W33 L=256)   264     72        0 (1391/ 516-1453/ 564)  45.0 W31-W34 S27-S24
      9 1547  488 (S29W46 L=269)   361     36        1 (1528/ 466-1568/ 509) -21.8 W43-W48 S30-S27
     10  515  417 (S44E40 L=183)   565     66        2 ( 488/ 385- 547/ 452)  40.6 E44-E35 S47-S41
  

    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