Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  701 1657 (N41E34 L=291)  4638    541       11 ( 474/1518- 893/1776) -10.3 E80-E14 N31-N47
      2 1087 1543 (N36W00 L=325)  1704    291        4 ( 961/1496-1227/1575)   9.5 E08-W10 N31-N38
      3  374 1493 (N26E56 L=269)   761    136        1 ( 326/1449- 419/1539) -14.9 E64-E49 N23-N28
      4 1351 1428 (N30W20 L=344)   531    101        1 (1312/1384-1388/1478) -51.3 W17-W22 N26-N32
      5  214 1412 (N19E73 L=251)  1779    111        4 ( 183/1367- 239/1467)  90.0 E79-E68 N15-N22
      6  316 1216 (N08E51 L=273)  2957    490        5 ( 272/1010- 401/1449) -87.9 E54-E46 S05-N22
      7 1169 1363 (N24W07 L=331)   534    121        2 (1148/1301-1190/1427)  83.7 W05-W08 N19-N28
      8 1588 1247 (N19W37 L=  1)   247     68        1 (1546/1235-1624/1256)   8.1 W33-W40 N18-N20
      9 1887 1215 (N19W70 L= 34)   706     71        2 (1870/1187-1907/1249)  26.6 W66-W74 N16-N20
     10 1753 1218 (N18W52 L= 16)   320     62        1 (1726/1199-1779/1236)   9.5 W49-W55 N17-N19
     11 1871 1160 (N15W66 L= 30)  1335    104        3 (1820/1132-1914/1192)   0.0 W58-W74 N13-N17
     12 1679 1050 (N08W44 L=  9)   259     79        0 (1648/1014-1698/1080) -45.0 W41-W45 N05-N09
     13 1878  959 (N03W66 L= 30)   740     78        1 (1845/ 928-1917/ 984)  15.3 W61-W72 N00-N04
     14 1543  920 (S01W34 L=359)   272     62        0 (1519/ 901-1570/ 948) -14.0 W32-W36 S02-N00
     15 1257  657 (S19W18 L=342)  1019    217        2 (1165/ 625-1359/ 694)  22.6 W12-W24 S21-S16
     16  931  513 (S32E02 L=322)  1072    245        5 ( 812/ 450-1034/ 552)  28.1 E11-W04 S37-S28
     17  638  388 (S45E28 L=297)   408     54        0 ( 613/ 371- 672/ 414) -26.6 E29-E24 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