Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1413 1494 (N22W38 L= 14)   478     71        1 (1389/1445-1434/1534)  33.7 W35-W41 N18-N23
      2  891 1486 (N37W05 L=342)   364    113        0 ( 846/1458- 934/1527)  14.0 W01-W09 N35-N38
      3 1162 1473 (N28W22 L=358)   348    174        0 (1080/1456-1243/1492) -20.0 W16-W27 N26-N30
      4  693 1445 (N40E10 L=326)   385     73        1 ( 660/1418- 736/1464) -26.6 E12-E06 N38-N41
      5 1639 1283 (N04W44 L= 21)   478     78        0 (1596/1259-1672/1314)  -8.1 W40-W47 N02-N04
      6  163 1264 (N40E70 L=267)  1017     72        1 ( 120/1239- 198/1283)  -8.1 E83-E62 N38-N41
      7  608 1181 (N26E21 L=315)  1467    322        3 ( 453/1144- 749/1207) -23.2 E33-E12 N20-N30
      8  236  928 (N19E57 L=279)  2748    508        5 ( 127/ 777- 453/1011)   7.9 E75-E36 N10-N23
      9 1654  975 (S13W36 L= 13)   313     68        1 (1630/ 938-1682/1007) -90.0 W35-W39 S15-S10
     10 1402  895 (S10W18 L=354)   872    217        0 (1331/ 810-1449/ 980)  28.6 W11-W22 S13-S06
     11 1633  888 (S18W33 L=  9)   542    101        3 (1599/ 849-1671/ 918)   9.5 W30-W36 S19-S15
     12 1148  657 (S17E03 L=333)   351    109        0 (1117/ 613-1182/ 712)  38.7 E05-W00 S18-S14
     13  929  648 (S12E15 L=321)   340     79        2 ( 901/ 618- 961/ 672)  14.0 E17-E13 S12-S10
     14  964  609 (S15E14 L=322)   254     46        0 ( 940/ 586- 990/ 630) -18.4 E15-E12 S15-S13
     15  598  363 (S22E49 L=287)  1012    185        1 ( 576/ 281- 630/ 453)  38.7 E57-E42 S27-S15
     16  379  399 (S16E68 L=268)   825     95        3 ( 343/ 357- 408/ 430)  78.7 E70-E64 S18-S13
     17 1175  300 (S42E14 L=322)   407     52        0 (1158/ 266-1193/ 335)  56.3 E15-E11 S45-S39
     18 1401  261 (S56W02 L=338)  1491    190        3 (1318/ 235-1481/ 300) -23.6 E06-W11 S61-S50
  

    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