Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  539 1641 (N17E51 L= 55)   747    225        3 ( 513/1568- 570/1714) -45.0 E56-E46 N11-N21
      2  871 1536 (N18E24 L= 82)   733    225        0 ( 845/1444- 896/1631) -63.4 E26-E20 N12-N24
      3 1248 1531 (N29E02 L=104)   558    103        0 (1213/1472-1272/1571) -51.3 E04-W01 N25-N31
      4 1057 1406 (N15E09 L= 97)   572    171        2 (1029/1334-1077/1488) -58.0 E13-E06 N10-N19
      5  477 1377 (S01E43 L= 63)   465    141        2 ( 430/1336- 515/1416)  -6.3 E48-E39 S02-N01
      6 1610 1173 (N19W31 L=136)   412    103        0 (1577/1133-1650/1202)  63.4 W28-W33 N15-N21
      7  847 1134 (S05E13 L= 93)   492    125        0 ( 815/1072- 870/1194) -82.9 E14-E11 S09-S01
      8  291  944 (S29E47 L= 59)  3862    702        6 ( 133/ 775- 407/1151)  52.1 E62-E40 S42-S15
      9 1605  992 (N09W35 L=140)   371     77        2 (1568/ 966-1644/1014)  31.0 W32-W37 N06-N10
     10  600  796 (S31E20 L= 86)  1063    332        2 ( 499/ 693- 695/ 906)  63.4 E25-E15 S39-S21
     11 1724  862 (N06W48 L=153)   739    221        1 (1632/ 826-1784/ 906)  42.0 W41-W51 N00-N09
     12  917  806 (S21W00 L=106)   589    175        1 ( 833/ 756-1003/ 869)  61.4 E03-W03 S26-S15
     13 1159  743 (S18W16 L=121)   501    207        0 (1071/ 703-1274/ 770)  45.0 W12-W21 S22-S13
     14 1644  717 (S05W47 L=153)   402     89        2 (1598/ 683-1682/ 744)  51.3 W44-W49 S07-S01
     15  930  631 (S31W06 L=112)   208     54        0 ( 899/ 618- 964/ 643)  36.9 W04-W08 S32-S29
     16 1528  576 (S15W45 L=151)   363     72        1 (1495/ 543-1564/ 593)  36.9 W42-W47 S17-S13
     17 1267  458 (S30W34 L=140)  1256    190        6 (1164/ 419-1380/ 494)  42.3 W28-W39 S34-S24
  

    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