Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  573 1759 (N27E61 L=150)   521     35        0 ( 550/1735- 596/1780) -18.4 E65-E56 N25-N28
      2  431 1693 (N19E69 L=142)   818     95        0 ( 392/1659- 470/1729)  71.6 E72-E67 N16-N22
      3  761 1639 (N23E36 L=175)   435    134        0 ( 710/1583- 787/1675) -18.4 E40-E31 N20-N25
      4  598 1639 (N19E47 L=164)   327     59        0 ( 575/1620- 627/1666) -18.4 E50-E44 N18-N20
      5  990 1442 (N17E14 L=197)  1314    542        1 ( 879/1261-1105/1664) -33.7 E29-E02 N09-N28
      6 1509 1419 (N32W19 L=230)   601    192        1 (1474/1330-1529/1496) -63.4 W15-W21 N26-N36
      7  406 1340 (S03E48 L=163)  1862    298        4 ( 311/1230- 467/1463)  82.9 E52-E43 S11-N05
      8  688 1448 (N10E32 L=179)   226     92        0 ( 650/1425- 737/1462)  38.7 E33-E28 N07-N11
      9 1340 1391 (N25W08 L=219)   542    188        1 (1293/1310-1365/1453) -45.0 W03-W11 N19-N27
     10  268 1360 (S04E62 L=149)  1056    150        3 ( 236/1284- 293/1416) -32.5 E69-E57 S08-S01
     11 1708 1374 (N37W38 L=249)   340     48        0 (1692/1339-1720/1399) -63.4 W35-W39 N34-N39
     12 1292 1256 (N15W09 L=220)   230     80        0 (1255/1234-1327/1277)   0.0 W06-W11 N14-N15
     13 1774 1111 (N22W47 L=258)   355     54        1 (1753/1076-1790/1137) -45.0 W44-W48 N20-N23
     14  417  983 (S23E37 L=174)   877    212        0 ( 342/ 924- 478/1065)  56.3 E43-E34 S28-S16
     15 1051  868 (S13W06 L=217)   781    268        3 ( 962/ 830-1155/ 906)  45.0 W00-W10 S17-S08
     16  947  888 (S15E01 L=210)   194     24        0 ( 922/ 871- 967/ 905)  63.4 E02-W00 S15-S13
     17 1192  833 (S11W15 L=226)   412    109        2 (1151/ 797-1248/ 870)  56.3 W12-W17 S13-S07
     18  771  724 (S29E07 L=204)  1041    386        1 ( 630/ 685- 949/ 794)  13.4 E17-W04 S31-S23
     19  658  427 (S53E07 L=204)   653    131        1 ( 595/ 401- 723/ 464)   4.4 E13-W00 S54-S51
  

    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