Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  709 1783 (N39E50 L=140)  1439    163        0 ( 637/1733- 770/1840)  -8.4 E68-E40 N36-N42
      2  731 1635 (N29E36 L=154)  1047    236        0 ( 618/1589- 820/1686)  47.7 E42-E32 N23-N34
      3 1577 1508 (N42W30 L=221)  1651    195        2 (1516/1413-1624/1586) -61.4 W24-W35 N35-N47
      4  419 1493 (N13E53 L=137)   671     90        0 ( 386/1452- 455/1533)  68.2 E55-E51 N10-N15
      5 1528 1302 (N26W27 L=217)  2307    662        2 (1366/1033-1638/1513) -41.6 W11-W38 N11-N36
      6 1114 1415 (N23E03 L=187)   450    198        1 (1041/1361-1171/1476) -21.8 E09-W01 N20-N25
      7 1005 1312 (N14E08 L=183)   417    127        3 ( 960/1274-1066/1346) -15.9 E10-E03 N12-N15
      8 1815 1179 (N25W55 L=245)  1162    107        0 (1792/1117-1833/1236) -82.9 W52-W57 N21-N29
      9  264 1037 (S17E53 L=137)  2489    457        5 ( 172/ 884- 348/1137)  48.4 E65-E46 S28-S09
     10  526  912 (S19E29 L=161)   225     65        0 ( 505/ 882- 544/ 944)  63.4 E30-E27 S21-S17
     11 1278  844 (S07W19 L=209)   260     96        0 (1237/ 825-1325/ 863)  31.0 W16-W21 S08-S05
     12 1496  795 (S05W34 L=225)   334    111        1 (1457/ 764-1547/ 815)  31.0 W31-W37 S07-S02
     13 1447  725 (S10W33 L=223)   276     58        1 (1425/ 697-1467/ 756) -45.0 W31-W34 S11-S08
     14 1432  658 (S14W34 L=225)   502    109        2 (1391/ 618-1465/ 690) -15.9 W30-W37 S16-S12
     15  356  618 (S43E45 L=145)  1178    142        2 ( 286/ 593- 427/ 651)   9.0 E56-E36 S45-S40
     16 1326  626 (S19W28 L=219)   380     83        0 (1279/ 613-1363/ 649)  38.7 W25-W30 S20-S16
     17 1446  537 (S21W40 L=231)   309     59        1 (1424/ 505-1473/ 561)   0.0 W37-W42 S23-S19
     18 1462  472 (S25W45 L=236)   710    232        0 (1381/ 400-1557/ 530)  63.4 W42-W50 S31-S19
  

    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