Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  816 1428 (N36E05 L=167)   257     78        1 ( 769/1411- 850/1446)  -9.5 E08-E02 N33-N36
      2 1050 1391 (N28W10 L=182)   389    109        2 ( 998/1349-1093/1420) -26.6 W06-W12 N26-N30
      3  303 1240 (N32E48 L=124)  2119    401        1 ( 234/1084- 403/1364) -49.9 E56-E39 N21-N41
      4  483 1016 (N16E35 L=138)   356    101        0 ( 443/1000- 526/1040) -26.6 E37-E31 N14-N17
      5 1454 1005 (S03W25 L=197)   325     98        1 (1397/ 990-1493/1022) -26.6 W21-W27 S04-S01
      6 1516  928 (S09W28 L=200)   228     98        1 (1471/ 903-1551/ 946)  -9.5 W24-W30 S11-S07
      7  599  876 (N06E29 L=143)   478    104        3 ( 559/ 833- 641/ 913) -38.7 E31-E26 N02-N08
      8  652  812 (N01E27 L=145)   392     85        0 ( 619/ 788- 675/ 850)  36.9 E28-E24 N00-N03
      9  629  506 (S17E37 L=135)   992    151        3 ( 578/ 440- 673/ 554) -77.5 E40-E34 S22-S13
     10  546  271 (S34E58 L=115)   854    224        2 ( 528/ 234- 564/ 308)  32.5 E63-E52 S37-S30
  

    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