Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1419 1447 (N18W34 L=286)  1395    220        5 (1314/1392-1487/1505)   7.6 W25-W40 N13-N20
      2 1718 1434 (N10W57 L=309)  1529    411        2 (1596/1381-1822/1478)  -9.1 W44-W69 N07-N12
      3 1029 1331 (N20W07 L=259)   577    202        1 ( 941/1306-1139/1356) -14.0 W01-W13 N18-N22
      4 1777 1096 (S10W51 L=303)   860    171        1 (1750/1026-1798/1191)  56.3 W47-W54 S13-S04
      5  496 1147 (N20E30 L=222)   349     63        1 ( 467/1111- 518/1177) -53.1 E32-E28 N17-N21
      6  365 1148 (N23E40 L=212)   701    195        2 ( 288/1126- 437/1173)  -8.7 E47-E34 N20-N24
      7  733 1123 (N14E15 L=237)   206     66        0 ( 698/1099- 774/1144) -45.0 E17-E12 N11-N15
      8 1737  810 (S26W43 L=295)  1753    301        3 (1637/ 692-1786/ 942)  42.7 W33-W47 S31-S18
      9  985  789 (S10E07 L=245)   234     63        1 ( 954/ 763-1023/ 804) -14.0 E08-E04 S11-S09
     10 1184  776 (S15W04 L=256)   197     44        1 (1162/ 756-1202/ 802)  33.7 W02-W05 S16-S13
     11 1090  698 (S18E03 L=249)   402    159        0 (1006/ 678-1175/ 722) -11.3 E08-W02 S19-S16
     12  758  682 (S12E23 L=229)   336    122        0 ( 696/ 638- 796/ 719) -50.2 E26-E21 S15-S08
     13 1544  481 (S44W26 L=278)   558     67        2 (1519/ 451-1567/ 523) -71.6 W23-W28 S46-S40
     14  604  319 (S32E49 L=204)  1662    307        3 ( 560/ 241- 645/ 414)  38.2 E57-E43 S36-S24
     15 1348  341 (S50W07 L=259)   838    135        2 (1324/ 294-1369/ 396)  63.4 W03-W09 S54-S45
     16 1138  320 (S45E11 L=241)   547    112        1 (1089/ 295-1176/ 357)   7.1 E15-E07 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