Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  752 1690 (N38E25 L= 38)   343     80        0 ( 716/1670- 793/1717) -26.6 E29-E21 N36-N40
      2  354 1571 (N29E58 L=  6)  2094    321        4 ( 296/1407- 406/1681) -59.0 E74-E52 N17-N37
      3  454 1408 (N17E41 L= 23)   887    223        2 ( 396/1317- 521/1509)  72.9 E44-E37 N11-N24
      4 1490 1449 (N25W30 L= 93)   268    103        0 (1439/1424-1536/1471) -15.9 W26-W33 N23-N26
      5  627 1327 (N13E27 L= 37)   768    270        0 ( 540/1257- 731/1385) -23.2 E33-E19 N08-N16
      6  549 1265 (N08E31 L= 32)   235     73        0 ( 517/1232- 566/1290) -36.9 E33-E29 N06-N10
      7 1934 1098 (N07W71 L=135)   970     72        2 (1911/1072-1959/1129) -63.4 W67-W77 N05-N09
      8  653  924 (S12E23 L= 41)   214     57        0 ( 624/ 909- 681/ 943)  14.0 E24-E20 S12-S10
      9 1222  879 (S12W13 L= 76)   503    198        1 (1162/ 815-1271/ 936)  58.0 W09-W15 S16-S08
     10  355  856 (S16E45 L= 18)   882    214        3 ( 266/ 813- 459/ 884)   3.2 E54-E36 S18-S14
     11  229  829 (S18E59 L=  5)   754     90        3 ( 204/ 790- 258/ 874) -68.2 E61-E55 S20-S15
     12 1326  736 (S20W21 L= 85)   742    164        4 (1279/ 676-1367/ 789) -50.2 W17-W24 S23-S16
     13 1800  704 (S18W61 L=125)  1057    126        0 (1744/ 670-1844/ 753)  36.9 W55-W65 S20-S14
     14 1107  719 (S22W07 L= 70)   409    169        0 (1034/ 696-1185/ 752) -15.3 W01-W12 S23-S20
     15 1463  500 (S34W37 L=101)   607     57        2 (1426/ 473-1496/ 527) -23.2 W33-W40 S35-S32
     16  937  410 (S44E03 L= 60)   309     99        2 ( 890/ 390- 977/ 437) -15.9 E07-W00 S45-S42
  

    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