Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1043 1752 (N46W03 L= 75)   401     32        0 (1012/1734-1073/1772)   0.0 W00-W05 N43-N47
      2 1077 1678 (N39W05 L= 77)   373     62        1 (1059/1648-1105/1716)  71.6 W03-W07 N36-N42
      3 1648 1562 (N31W50 L=122)  1508    202        1 (1566/1517-1738/1618) -28.3 W43-W58 N27-N34
      4  779 1573 (N32E16 L= 56)   716    143        0 ( 712/1525- 834/1601) -29.1 E21-E12 N28-N34
      5 1763 1501 (N27W60 L=132)   637     54        1 (1736/1475-1791/1530)  24.0 W55-W65 N24-N28
      6 1419 1471 (N24W27 L= 99)   703    159        1 (1348/1434-1511/1495)   0.0 W22-W34 N21-N25
      7 1595 1429 (N21W40 L=112)   660    122        0 (1565/1374-1620/1475)  74.1 W37-W42 N17-N24
      8 1846 1276 (N12W61 L=133)  3112    254        6 (1793/1167-1926/1381) -47.7 W53-W72 N05-N18
      9 1680 1245 (N10W44 L=116)   658    182        1 (1655/1153-1708/1330) -79.7 W42-W45 N04-N15
     10  407  725 (S19E43 L= 29)  4946    875        8 ( 303/ 375- 464/1087)  42.4 E89-E36 S40-N02
     11  230 1025 (S00E55 L= 17)   557    188        0 ( 182/ 961- 267/1074) -45.0 E60-E51 S04-N02
     12  190  852 (S10E62 L= 10)  5178    311       11 ( 123/ 707- 280/ 999)  73.6 E71-E52 S18-S01
     13  612  891 (S10E26 L= 46)   319     48        1 ( 573/ 870- 645/ 917)   0.0 E28-E23 S11-S08
     14  917  720 (S21E07 L= 65)  1485    487        1 ( 733/ 653-1016/ 873)  21.3 E19-W00 S25-S11
     15 1182  822 (S15W09 L= 81)   585    261        2 (1124/ 739-1261/ 862)  45.0 W05-W14 S20-S12
     16  289  753 (S17E53 L= 19)   659     91        0 ( 257/ 709- 318/ 796) -51.3 E55-E51 S19-S14
     17  233  762 (S16E59 L= 13)   593     57        2 ( 200/ 731- 261/ 783) -36.9 E62-E55 S17-S14
     18 1329  481 (S38W22 L= 94)  2528    594        6 (1088/ 372-1605/ 570)  16.7 W03-W43 S45-S31
     19 1174  250 (S56W14 L= 86)   394     73        1 (1131/ 233-1199/ 278)  33.7 W09-W17 S58-S53
  

    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