Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  538 1712 (N47E50 L=129)  3885    504        6 ( 428/1632- 673/1771)  16.7 E59-E35 N39-N51
      2 1125 1647 (N45W06 L=184)   310     48        1 (1101/1617-1147/1668)  63.4 W03-W07 N42-N46
      3  835 1631 (N42E19 L=160)   417     59        0 ( 804/1604- 867/1660)  45.0 E20-E16 N40-N44
      4 1236 1550 (N38W14 L=193)   265     53        0 (1213/1535-1266/1574) -26.6 W12-W16 N36-N39
      5  427 1365 (N21E44 L=134)   548    115        2 ( 400/1330- 461/1409) -45.0 E47-E40 N19-N24
      6  362 1267 (N15E48 L=131)  1058    232        2 ( 286/1229- 442/1292)   4.1 E55-E41 N12-N16
      7  344 1176 (N09E48 L=130)  1146    375        6 ( 257/1112- 450/1241) -17.5 E57-E38 N05-N12
      8  199 1156 (N07E64 L=115)   531    121        1 ( 180/1108- 216/1209) -56.3 E66-E60 N03-N09
      9  947 1056 (N04E05 L=174)   211     86        0 ( 915/1015- 985/1093) -45.0 E06-E02 N02-N06
     10  616  938 (S04E26 L=153)   250     78        0 ( 598/ 897- 640/ 979)  68.2 E26-E24 S06-S01
     11  384  791 (S15E44 L=135)  1575    260        4 ( 356/ 705- 427/ 898)  81.3 E45-E40 S20-S07
     12 1862  845 (S07W66 L=244)   784    105        1 (1839/ 805-1881/ 886) -45.0 W62-W68 S09-S04
     13  964  778 (S13E03 L=176)   245     73        2 ( 940/ 741- 990/ 812)  53.1 E04-E01 S15-S10
     14 1424  581 (S24W30 L=208)   305     69        0 (1387/ 567-1459/ 598)  -9.5 W26-W32 S25-S23
     15 1231  431 (S36W18 L=197)   857    145        1 (1164/ 408-1306/ 462)  -4.8 W12-W24 S37-S33
     16  826  401 (S40E13 L=165)  1335    290        2 ( 688/ 356- 958/ 436)  12.3 E25-E02 S44-S36
  

    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