Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  544 1759 (N40E51 L=118)   593     75        1 ( 521/1721- 564/1795) -28.3 E59-E46 N36-N43
      2  572 1600 (N28E38 L=131)   737    233        3 ( 495/1568- 682/1630)   7.6 E45-E29 N25-N29
      3  295 1480 (N19E59 L=110)  1055    130        4 ( 262/1440- 350/1519)  39.8 E63-E54 N16-N21
      4 1531 1440 (N25W31 L=200)   836    157        0 (1471/1388-1589/1493) -32.0 W26-W34 N22-N28
      5  682 1473 (N20E26 L=143)   217     45        1 ( 657/1449- 708/1492) -14.0 E28-E24 N18-N20
      6 1080 1353 (N15W01 L=170)   651    147        4 (1008/1298-1126/1399)  26.6 E03-W03 N11-N18
      7  846 1286 (N09E13 L=156)   902    302        2 ( 719/1210- 935/1390) -32.7 E21-E07 N04-N14
      8 1475 1360 (N19W26 L=195)   235     53        2 (1453/1334-1497/1390) -76.0 W24-W27 N17-N21
      9 1915 1168 (N14W67 L=236)   726     68        3 (1883/1129-1941/1197) -31.0 W62-W71 N11-N15
     10  626  987 (S10E24 L=145)   687    229        0 ( 546/ 913- 711/1050)  39.3 E29-E18 S14-S05
     11  178  898 (S17E63 L=106)  1514    226        6 ( 126/ 810- 221/ 987)  48.0 E70-E56 S22-S11
     12  612  852 (S18E25 L=144)   303    135        0 ( 581/ 802- 634/ 913) -81.9 E26-E23 S21-S14
     13  938  818 (S18E04 L=165)  1632    753        3 ( 631/ 753-1234/ 894)  13.7 E23-W14 S23-S13
     14 1777  729 (S14W57 L=226)   508    159        0 (1739/ 685-1798/ 808)  82.9 W53-W58 S17-S09
     15  601  661 (S30E26 L=143)   359     71        2 ( 585/ 625- 623/ 696)  78.7 E27-E24 S32-S27
     16  474  629 (S33E37 L=132)   328     74        1 ( 445/ 597- 507/ 663) -53.1 E39-E34 S34-S30
     17 1294  499 (S35W25 L=194)   236     82        0 (1255/ 484-1337/ 514)  21.8 W22-W28 S35-S33
     18 1293  398 (S42W29 L=198)   798    230        2 (1219/ 357-1407/ 456)  -6.3 W21-W39 S44-S38
  

    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