Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1618 1528 (N15W55 L=192)   380     52        1 (1589/1507-1649/1547) -26.6 W53-W57 N14-N17
      2 1055 1266 (N19W08 L=146)   514    134        1 ( 992/1246-1130/1283) -29.7 W04-W12 N16-N20
      3 1877 1210 (S10W64 L=202)  2554    307        8 (1822/1122-1922/1274) -57.5 W58-W73 S15-S04
      4  874 1183 (N19E04 L=133)   299    141        0 ( 817/1136- 925/1218) -50.2 E07-E02 N15-N21
      5 1834 1101 (S14W55 L=192)  1407    123        2 (1799/1043-1877/1165) -71.6 W51-W59 S18-S09
      6  302  977 (N21E49 L= 89)  2056    381        2 ( 192/ 840- 388/1147) -55.7 E60-E43 N11-N33
      7  637 1094 (N20E21 L=117)   211     86        0 ( 605/1068- 670/1126)  18.4 E23-E17 N19-N21
      8 1781 1023 (S16W48 L=185)   700    101        1 (1747/ 975-1807/1075) -90.0 W44-W50 S19-S13
      9  519  990 (N17E31 L=107)   368    129        0 ( 468/ 955- 581/1022) -45.0 E33-E27 N13-N19
     10 1263  830 (S11W08 L=145)   748    338        4 (1225/ 680-1282/ 994)  57.1 W01-W12 S18-S01
     11 1616  903 (S17W32 L=169)   250     55        0 (1582/ 893-1650/ 914) -36.9 W29-W34 S18-S15
     12 1373  542 (S31W07 L=145)  2952    539        5 (1230/ 373-1539/ 746)  27.3 E07-W22 S38-S23
     13  590  601 (S06E37 L=100)   747    249        0 ( 544/ 482- 612/ 692)  59.7 E43-E34 S12-N00
     14  431  617 (S02E48 L= 89)   724    129        1 ( 365/ 578- 487/ 657) -45.0 E52-E44 S05-N01
     15 1194  561 (S24E03 L=134)   271     95        0 (1174/ 513-1219/ 592)  45.0 E05-E01 S27-S22
     16  787  426 (S22E32 L=105)  1182    339        3 ( 762/ 305- 840/ 587)  56.3 E39-E24 S30-S12
     17  970  498 (S22E18 L=119)   391    111        1 ( 929/ 461- 990/ 555)  63.4 E21-E15 S24-S18
  

    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